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Abstract 

Experimental data are presented on particle correlations and fluctuations in various 
high-energy multiparticle collisions, with special emphasis on evidence for scaling-law 
evolution in small phase-space domains. The notions of intermittency and fractality as 
related to the above findings are described. Phenomenological and theoretical work on 
the subject is reviewed. 
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CONTENTS 



Of the achieved triumph pangs and tricks 
Are just tightly stretched bow-strings. 

B. Pasternak 



Chapter 1 
Introduction 



Recent years have witnessed a remarkably intense experimental and theoretical activity 
in search of scale-invariance and fractality in multihadron production processes, for short 
also called "intermittency" . These investigations cover all types of reactions ranging from 
e + e~ annihilation to nucleus-nucleus collisions, up to the highest attainable energies. 
The creation of soft hadrons in these processes, a major fraction of the total cross section, 
relates to the strong-coupling long-distance regime of Quantum Chromodynamics (QCD), 
at present one of the least explored sectors in the whole of high-energy particle physics. 

A primary motivation is the expectation that scale-invariance or self-similarity, anal- 
ogous to that often encountered in complex non-linear systems, might open new avenues 
ultimately leading towards deeper insight into long-distance properties of QCD and the 
unsolved problem of colour confinement. 

History shows that studies of fluctuations have often triggered significant advances in 
physics. In the present context, it was the observation of "unusually large" particle density 
fluctuations, reminiscent of intermittency spikes in spatio-temporal turbulence, which 
prompted the pioneering suggestion to investigate the pattern of multiplicity fluctuations 
in multihadron events for ever decreasing domains of phase-space. Scale-invariance or 
fractality would manifest itself in power-law behaviour for scaled factorial moments of the 
multiplicity distribution in such domains. It is important to stress here that, in practice, 
one deals with the problem of evolution of particle number distributions for ever smaller 
bins and intermittent behaviour implies that, for small phase space bins, the distributions 
become wider in a specific way. The same problem can be stated as an increasing role of 
correlations within a small phase space volume. 

Through a multitude of increasingly sophisticated experimental studies of factorial 
moments, much new information has been gathered in a surprisingly short time. This 
work indeed confirms approximate power behaviour down to the experimentally possible 
resolution, especially when carried out in two- and three-dimensional phase space. 

The proposal to look for intermittency also has triggered a thorough revival of in- 
terest in the old subject of particle correlations, by experimentalists and theorists alike. 
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CHAPTER 1. INTRODUCTION 



The need for greater sensitivity in measurements of correlation functions has directly 
inspired important work on refined analysis techniques. A promising and long overdue 
systematic approach to correlation phenomena of various sorts, including Bose-Einstein 
interferometry, is finally emerging. 

The large body of experimental observations now available is calling for satisfactory 
explanation and, indeed, theoretical ideas of all sorts abound. 

The level of theoretical understanding is quite different for the various types of colli- 
sion processes. In e + e~ annihilation, parton cascade models based on leading-log QCD 
have met considerable success and good overall description of multiplicity fluctuations is 
claimed. Closer inspection, nevertheless, reveals potentially serious deviations from the 
data, thus requiring further study. 

For other processes, in particular hadron initiated collisions, models are faced with 
large and partly unexpected obstacles. This may be a reflection of insufficient knowledge 
of the reaction dynamics, although present evidence points to hadronization as the main 
culprit. 

Within the framework of perturbative QCD, results of considerable interest on the 
emergence of power behaviour and multifractality have been obtained. However, these 
are asymptotic in nature and most likely quite unrelated to present-day experiment. Be- 
ing related to the mechanism of confinement, not surprisingly, the role of hadronization 
remains unclear. 

Random self-similar multiplicative branching models have inspired much of the original 
work on intermittency. Among many scale-invariant physical systems, the cascade process 
is a particularly natural candidate for the description of strong fluctuations self-similar 
over a wide range of scales. It finds support in the cascade nature, not only of perturbative 
QCD, but also of the subsequent hadronization. However, further work is needed to help 
understand the details of the process. 

Alternatively, "classic" and extensively studied possibilities are scale-invariant systems 
at the critical point of a high-order phase transition. This subject has attracted particular 
attention in view of potential application to quark-gluon plasma formation in heavy-ion 
collisions. 

This paper contains a review of the present status of work on intermittency and correla- 
tions as performed over the last years. In Chapter 2 we introduce the necessary formalism 
and collect useful results and relations widely scattered in the literature. Chapter 3 de- 
scribes experimental data on correlations in various experiments and discusses predictions 
of popular models. Chapter 4 is devoted to data and models on the subject of particle 
fluctuations and the search for power laws. Chapter 5 gives an overview of the many 
theoretical ideas related to the problem of multiplicity scaling and fractality. Conclusions 
are summarized in Chapter 6. 



Chapter 2 
Formalism 



2.1 Definitions and notation 

In this section, we compile and summarize definitions and various relations among the 
physical quantities used in the sequel of this paper. No originality is claimed in the 
presentation of this material. It merely serves the purpose of fixing the notation and 
assembling a number of results scattered throughout the literature. 

2.1.1 Exclusive and inclusive densities 

We start by considering a collision between particles a and b yielding exactly n particles 
in a sub- volume Q of the total phase space fitot- Let the single symbol y represent 
the kinematical variables needed to specify the position of each particle in this space 



(for example, y can be the cm. rapidityQ variable of each particle and fl an interval of 
length 8y). The distribution of points in Q can be characterized by continuous probability 
densities P n (yi, • • • , y n )', n — 1,2, . . .. For simplicity, we assume all final-state particles to 
be of the same type. In this case, the exclusive distributions P n {y\, ■ ■ ■ , y n ) can be taken 
fully symmetric in y±, . . . , y n \ they describe the distribution in Q when the multiplicity is 
exactly n. 

The corresponding inclusive distributions are given for n — 1,2,... by: 




Pn(Vl, ■ ■ ■ ,Un) 



P n (yi, ■ ■ ■ ,y n 




1 The rapidity y is defined as y = ~ ln[(E +pi,)/(E — pi,)], with E the energy and the longitudinal 
component of momentum vector p along a given direction (beam-particles, jet-axis, etc.); pseudo-rapidity 
is defined as r) = \ ln[(p + p L )/ (p - p L )] 
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The inverse formula is 

Pn(yi,---,y n ) = Pn(yi,---,Vn) 



oo i „ m 

+ E (-i) m ^ / Pn+m(yi, ■ ■ ■ , v», y'l, ■ ■ ■ , y'J U dy' t . (2.2) 



p n (yi, . . . , y n ) is the probability density for n points to be at y x , ... , y n , irrespective of 
the presence and location of any further points. The probability Pq of multiplicity zero is 
given by 

00 „ n 

This suggests to define po — 1 i n ( |2.1| ). It is often convenient to summarize the above 
results with the help of the generating functional!! 



oo i „ n 

nl Jn t=i 



dy t , (2.4) 



where z(y) is an arbitrary function of y in Q. The substitution 

z(y) = l + u(y) (2.5) 
gives through fl2.1| ) the alternative expansion 

oo i . n 

^ ncl Ky)] = l + ^- / Pn (y 1 ,...,y n )u(y 1 )...u(y n ) J[dy t (2.6) 

£i n] Jn t=i 

and the relation 

£ incl [z{y)\ = £ cxcl [z(y) + 1] . (2.7) 
From (|2.4|) and ( |2.7|) one recovers by functional differentiation: 

Qngexcl 



and 



P n (vi, ■ ■ ■ ,yn) 



pn{yi, ■ ■ -,y n ) 



dz(yi) . ..dz(y n ) 

Qngind 



du(yx) . . . du(y n ) 



z=0 



u=0 



(2.8) 



(2.9) 



To the set of inclusive number-densities p n corresponds a sequence of inclusive differential 
cross sections: 

— dcT = Pl (y)dy, (2.10) 

^inel 

1 



^incl 



d a = p 2 (yi, y 2 ) dy!dy 2 . (2.11) 



2 The technique of generating functions has been known since Euler's time and was used for functionals 
by N.N. Bogoliubov in statistical mechanics already in 1946 see also Q 
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Integration over an interval Q in y yields 

/ Pi(y)dy = (n) 
Jn 

/ / P2(Vi,y2)dy 1 dy 2 = (n(n - 1)) 
Jn Jn 

f d Vl ... [ dy qPq ( yi , ...,y q ) = (n{n - 1) . . . (n - q + 1)) , (2.12) 
Jn Jn 

where the angular brackets imply the average over the event ensemble. 
2.1.2 Cumulant correlation functions 



The inclusive g-particle densities p g (yi, ■ ■ ■ ,y q ) in general contain "trivial" contributions 
from lower-order densities. Under certain conditions, it is, therefore, advantageous to 
consider a new sequence of functions C q (yi, . . . ,y q ) as those statistical quantities which 
vanish whenever one of their arguments becomes statistically independent of the others. 
It is well known that the quantities with such properties are the correlation functions-also 
called (factorial) cumulant functions-or, in integrated form, Thiele's semi- invariants [|J. A 
formal proof of this property was given by Kubo ||] (see also Chang et al. ||). The cumu- 
lant correlation functions are defined as in the cluster expansion familiar from statistical 
mechanics via the sequence 



Pi(l) = QlCI), (2.13) 

p 2 (l,2) = C7 1 (l)Ci(2) + C 2 (l,2), (2.14) 
#,(1,2,3) = C 1 (1)C 1 (2)C 1 (3) + C 1 (1)C 2 (2,3) + C 1 (2)C 2 (1,3) + 

+ C 1 (3)C 2 (1,2)+C 3 (1,2,3); (2.15) 

and, in general, by 

p m (l,...,m) = E E lCiO---C 1 Q\p 2 (,)-yC 2 (,)\ -- 

^ m h factors i 2 factors 

•••[C m (,...,).--C m (,. ..,)]. (2.16) 

V v ' 

i m factors 

Here, U is either zero or a positive integer and the sets of integers {li} m satisfy the 
condition 

Y J ik = m. (2.17) 

i=l 

The arguments in the Cj functions are to be filled by the m possible momenta in any order. 
The sum over permutations is a sum over all distinct ways of filling these arguments. For 
any given factor product there are precisely [|7| 

^ r2 18^ 
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terms. The complete set of relations is contained in the functional identity: 

g™ l [u(y)]=e W {g[u(y)}} , (2.19) 

where 

Pi(y)u(y) dy + — Cg(y 1 ,...,y q )u{y 1 )...u(y q )Y[dy i . (2.20) 
2 q. Jn i=1 



It follows that 

C q (yi,...,y q ) - 



du{y x ) . . . du(y n ) 
The relations (|2.16 ) may be inverted with the result: 



(2.21) 



u=0 



(72(1,2) = p 2 (l,2)- Pl (l) Pl (2) , 
C 3 (l,2,3) = p 3 (l,2,3)-^p 1 (l)p 2 (2,3) + 2p 1 (l)p 1 (2)p 1 (3) , 

(3) 

C 4 (l,2,3,4) = p 4 (l,2,3,4)-^ Pl (l)p3(l,2,3)-^p 2 (l,2)p 2 (3,4) 

(4) (3) 

+ 2^p 1 (l)p 1 (2)p 2 (3,4)-6p 1 (l)p 1 (2)p 1 (3)p 1 (4). (2.22) 

(6) 

In the above relations we have abbreviated C q (yi, . . . ,y q ) to C q (l, 2, . . . , g); the summa- 
tions indicate that all possible permutations have to be taken (the number under the 
summation sign indicates the number of terms). Expressions for higher orders can be 
derived from the related formulae given in ||. 

It is often convenient to divide the functions p q and C q by the product of one-particle 
densities. This leads to the definition of the normalized inclusive densities and correla- 
tions: 

r q (yu---,y q ) = Pq(vi,---,Vq)/pi(yi)--Pi(v q ), (2-23) 
Kg(Vi,---,Vq) = C q (y u ...,y q )/p 1 (y 1 )...p 1 (y q ). (2.24) 



From expression ( 2.19|) it can be deduced that, at finite energy, an infinite number of 



C q will be non-vanishing: The densities p q vanish for q > N, where N is the maximal 
number of particles in Q allowed e.g. by energy-momentum conservation. As a conse- 
quence, the functional Q is a "polynomial" in u(y). This in turn requires the exponent 
in ( |2.19| ) to be an "infinite series" in u(y). In other words, the higher-order correlation 
functions must cancel the lower-order ones that contribute to a vanishing density func- 
tion. Phenomenologically, this implies that it is meaningful to use correlation functions 
C q only if the number of correlated particles in the considered phase-space domain Q is 
considerably smaller than the average multiplicity in that region . These conditions are 
not always fulfilled in present-day experiments for very small phase-space cells, with the 
exception of perhaps AA-collisions. 
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2.1.3 Correlations for particles of different species 



The generating functional technique of Sect. 2.1.1 can be extended to the general situation 
where several different species of particles are distinguished. This will not be pursued here 
and we refer to the literature for details P, |TT)|-|T2"|] . Considering two particle species a and 
b, the two-particle rapidity correlation function is of the form: 

Cf{ Vl) y 2 ) = pf( yi ,y 2 ) - /pl(yi)p b (y 2 ), (2.25) 

with 

Pt(yi) = f-] pf{VuV2) = f^— (2.26) 

0"inel d?/i CT ine l d?/id?/ 2 

Here, y\ and y 2 are the cm. rapidities, <r inel the inelastic cross section and a, b represent 
particle properties, e.g. charge. 

The normalization conditions are: 

ft(yi)d yi = (n a > ; / / pf (y u y 2 )d yi dy 2 = (n a (n h - 5 ab )> , (2.27) 



C% b (y 1 ,y 2 )dy 1 dy 2 = {n !i (n h -5 &b ))-f{n s ){n h ) , (2.28) 

where 5 ab = for the case when a and b are particles of different species and 5 ab = 1 for 
identical particles, and n a and n h are the corresponding particle multiplicities. 

Most experiments use 

f = l , (2.29) 
so that the integral over the correlation function (equal to the ratio n 2 /k of the negative 



binomial parameters [|1|) vanishes for the case of a Poissonian multiplicity distribution. 
Other experiments use 

j = M^)> (2 . 30) 

to obtain a vanishing integral also for a non- Poissonian multiplicity distribution. 

To be able to compare the various experiments, we use both definitions and denote 
the correlation function Cf'{yi,y 2 ) when following definition ( p. 29 ) and C 2 ab {yi,y 2 ) when 



following definition (|2.30|) . We, furthermore, use a reduced form of definition ( |2.30| ), 

C?{yi,y2) = C' 2 ah (yi,y2)/(n a (n h - 5 ab )>. (2.31) 



The corresponding normalized correlation functions 

Cf(y 

fpi(yi)pi{y. 



K?{y^) = f^ im t\ (2.32) 
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follow the relations 

K 2 ' = j(K 2 + l)-l , (2.33) 

and K 2 is denned as K 2 = K 2 . These are more appropriate than C 2 when comparisons 
have to be performed at different average multiplicity and are less sensitive to acceptance 
problems. 

The correlation functions defined by expressions ( |2.25| )-( p.33| ), contain a pseudo-correl- 
ation due to the summation of events with different charge multiplicity n and different 
semi-inclusive single-particle densities pi . 

The relation between inclusive and semi-inclusive correlation functions has been care- 



fully analyzed in |14|. Let a n be the topological cross section and 

P n = a n /Za n . (2.34) 

The semi-inclusive rapidity single- and two-particle densities for particles a and b are 
defined as 

p[ n \y) = and A n \y uy2 ) = • (2.35) 

cr n dy a n dy 1 dy 2 

The inclusive correlation function C 2 (yi,y 2 ) can then be written as 

C 2 ( yi ,y 2 ) = C s (y 1 ,y 2 )+C h (y 1 ,y 2 ) , (2.36) 

where 

C a { yi ,y 2 ) = EP n (fr\ yi ,y 2 ) (2.37) 

C L (y u y 2 ) = ZP n ApW( yi )ApW(y 2 ) (2.38) 



with Cf V.jfe) = P ( 2 \yi,y2) - P t L \yi)p ( r\y2) and ApW(y) = P f\y)- Pl {y). In O 



C§ is the average of the semi-inclusive correlation functions (often misleadingly denoted as 
"short-range" ) and is more sensitive to dynamical correlations. The term Cl (misleadingly 
called "long-range" ) arises from mixing different topological single-particle densities. 

A normalized form of Cg can be defined as 



K s (y h y 2 ) = C fy^ () = ^nPnf^y.) 

' En Pnpt\yi)p ( l\y2) En Pn/tfV^Gfe) 



C' s and C*s and their normalized forms K§ and Kg are defined accordingly, with the 
averages (n) and (n a (n h — <5 ab )) replaced by n and n a (n h — 5 ab ), respectively. 

Analogous expressions may be derived for three-particle correlations. They are dis- 
cussed in Sect. 3.4. 
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2.1.4 Factorial and cumulant moments 



When the parametric function z(y) is replaced by a constant z, the generating functionals 
reduce to the generating function for the multiplicity distribution. Indeed, the probability 
P n for producing n particles is given by 



P n = (CV^el (2.40) 



and we have 



G{z) = P„{1 + z) n = Q cxcl [z + 1} =g' mcl [z] (2.41) 

n=0 

00 z q r 

= 1 + E — / PqiVl, •••'%) d Vl • • ( 2 - 42 ) 
0=1 Jn 



oo q 



= l + E^^i- ( 2 - 43 ) 
=i y- 

The F are the unnormalized factorial (or binomial) moments 

F q = (n [q] ) = (n(n-l)...(n-q + l)) 

= dyi... dy q p q {yi,...,y q ) 
Jn Jn 

= Y / Pnn(n-l)...(n-q + l). (2.44) 

n 

This relation can (formally) be inverted. If P n = for n > N then an approximation for 
P n is given by: 

i N~n rp 

^ = ^E(- 1 ) J ^f i (n = 0,1, . . . N), (2.45) 

and P n is included between any two successive values obtained by terminating the sum 
at j = s and j = s + 1, respectively. 



In (|2.44j ) n denotes the multiplicity in f2 and the average is taken over the ensemble 
of events. All the integrals are taken over the same volume Q such that ?/; 6 H Vi 6 
{1, . . . , q}. Using the correlation-function cluster decomposition, one further has 



logG(*) = <«>* + i;-r/«- (2 " 46) 

q=2 1- 

The f q are the unnormalized factorial cumulants, also known as Mueller moments || 

fq= [ d Vl ... [ dy g C q ( yi ,...,y q ), (2.47) 
Jo, Jn 



16 



CHAPTER 2. FORMALISM 



the integrations being performed as in ( |2.44 ). The quantities F q and f q are easily found 
if G(z) is known: 



d q G(z) 



z=0 

dnogG(z 



and 



dzi 
1 d q G(z) 



2=0 



q\ dz q 



(2.48) 
(2.49) 

(2.50) 



z=-l 



Using Cauchy's theorem, this can also be written as 

1 r G{z) 



2m J (1 + z 



,n+l 



dz. 



(2.51) 



where the integral is on a circle enclosing z = —1. Equation ( [2.51| ) is sometimes useful in 
deriving asymptotic expressions for P n in terms of factorial moments or cumulants fT5|,|8]]. 

As a simple example, we consider the Poisson distribution 



P n = e" (n) 



(n) n 



for which 



G(z) = "£P n (l + z) n = exp{(n)z} 
o 



(2.52) 



showing that f q = for q > 1. In that case one has: 



F q = (n(n - 1) . . . (n - q + 1)) = (n) q . 



(2.53) 



The expressions of density functions in terms of cumulant correlation functions, and 
the reverse relations, are duplicated for their integrated counterparts. They follow directly 
from the equations: 



oo z q 



oo z q 



q =i y- q =2 y- 



or 



oo z q 



g=i y- 



q=2 y- 



(2.54) 



(2.55) 



by expanding either the exponential in (p.54| ) or the logarithm in ( [2.55|) and equating the 
coefficients of the same power of z. One finds |J: 
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h 

F, 
Fi 

h 

and in general: 



h , 

f2 + fl , 

h ■ • ./? , 

h + + 3/ 2 2 + Qf 2 fi + ft , 

h + 5/4/1 + 10/3/2 + 10/3A 2 + i5/ 2 2 /i + 10/2A 3 + A 5 



(2.56) 
(2.57) 



with the summation as in (|2.16|) and J2i=i ih — Q- 
The latter formula can also be written as: 

q-l 



E 

1=0 



q-l 
I 



fq-l Fl 



(2.58) 



(with F = 1, f Q = 0) and is well-suited for computer calculation. An equivalent relation 
was derived in [§]. The (ordinary) moments: 



n=0 

may be derived from the moment generating function 

00 

M(s) = Ee™P B , 



since 



n=0 



d 9 M(z) 



dz« 



2=0 



We note the useful relations 



M(z) 
G(z) 



G(e z -l) , 
M(log(l + z)) 



^2.59) 



^2.60) 



(2.61) 



(2.62) 
(2.63) 



Moments and factorial moments are related to each other by series expansions. From 



the identities 16 



n 



n(n-l)...(n-q + l) = £ S^ m) 

-- Y,S^ m) n{n-l)...{n-m + l) 



rr 



(2.64) 
(2.65) 



m=0 
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where and are Stirling numbers of the first and second kind, respectively, follows 
directly: 



Mr, 



m=0 
9 



Mg 



E ^ 



(2.66) 
(2.67) 



m=0 



Cumulants K q can be defined in terms of the moments p q in the standard way |L7|, 
||. They obey relations identical to (|2.56| ). The cumulants are integrals of the type 
( |2.47| ) of differential quantities known as density moments. These are discussed in [ |T8| 
[XQfl . Relations expressing central moments in terms of factorial moments via non-central 
Stirling numbers are derived in EO . 



2.1.5 Cell- averaged factorial moments and cumulants; 
generalized moments 



In practical work, with limited statistics, it is almost always necessary to perform averages 
over more than a single phase-space cell. Let Q m be such a cell (e.g. a single rapidity 
interval of size 5y) and divide the phase-space volume into M non-overlapping cells fl m 
of size 5Q, independent of m. Let n m be the number of particles in cell fl m . Different 
cell-averaged moments may be considered, depending on the type of averaging. 



Normalized factorial moments |pl] , |22|l , which have become known as vertical moments, 
are defined asi 



Ai 



— y 



(n m (n m 1) 



(n r 



-9+1)) 



(2.68) 
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[2.69) 



The full rapidity interval AY is divided into M equal bins: AY = M8y; each yi is within 
the (fy-range and (n m ) = p m 5y = J 5y pi(q)dy. 

One may also define normalized horizontal moments by 



F?(5y) = 



M 



[n r , 



9+1)) 



m=l 



(n m y 



^2.70) 



Here and in the following we consider rapidity space for definiteness 



2. 1 . DEFINITIONS AND NOTATION 



19 



with n m = £ m n m /M; (n m ) = (n)/M; n = J2 m n m- 

Horizontal and vertical moments are equal if M — 1. Vertical moments are normalized 
locally and thus sensitive only to fluctuations within each cell but not to the overall shape 
of the single-particle density. Horizontal moments are sensitive to the shape of the single- 
particle density in y and further depend on the correlations between cells. To eliminate 
the effect of a non-flat rapidity distribution, it was suggested to either introduce correction 
factors [23] or use "cumulative" variables which transform an arbitrary distribution into 



a uniform one [24, 25] 



Likewise, cell-averaged normalized factorial cumulant moments may be defined as 



K q {8y) 



M(5yY 



M 

E 



L Sy 



Cgjyi, ■•■■>y q ) 



(2.71) 



They are related [26] to the factorial moments b 



F 2 
F 3 
F 4 



l + K 2 , 

1 + 3K 2 + K 3 , 

1 + QK 2 + 3Kl + AK 3 + K A . 



^2.72) 



In F 4 and higher-order moments, "bar averages" appear. They are defined as AB = 
J2A m B m /M. 

m 

Besides factorial and cumulant moments, other measures of multiplicity fluctuations 



have been proposed. In particular, G-moments f27]] — known in statistics as frequency 
moments H] — were extensively used to investigate whether multiparticle processes possess 
(multi)fractal properties p^,^9|. G-moments are defined as 

M M 

G q = ^2'p q m , Pm = n m /n, n = ^2n m . (2.73) 

m=l m=l 

Also here, n m is the number of particles in bin m, the absolute frequency; n is the total 
multiplicity in an initial interval and M is the number of bins at "resolution" M. Bins 
with zero content ( "empty bins" ) are excluded in the sum, so that q can cover the whole 
spectrum of real numbers. For q negative, G q is sensitive to "holes" in the rapidity 



distribution of a single event. Note that p m in (|2.73|) is not a probability but a relative 
frequency or "empirical measure" in modern terminology. For small n, G-moments are 
very sensitive to statistical fluctuations ("noise"), especially for large M. This seriously 
limits their potential. In attempts to reduce this noise-sensitivity, modified definitions 
have been proposed in BT 



The higher-order relations can be found in 
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2.1.6 Multivariate distributions 



The univariate factorial moments F q characterize multiplicity fluctuations in a single 
phase-space cell and thus reflect only local properties. More information is contained in 
the correlations between fluctuations (within the same event) in two or more cells. This 
has led to consider multivariate factorial moments. For non-overlapping cells, the 2-fold 
factorial moments, also called correlators, are defined as: 

F m = (ngnB), (2.74) 

where n m (n m >) is the number of particles in cell m (cell m'). A normalized version of the 
two-fold correlator is discussed in |21] and defined as: 



F pq = ^ffl . (2.75) 

P Q 

For reasons of statistics, these quantities are usuallv averaged over many pairs of cells, 
keeping the "distance" (D) between the cells constant^. This averaging procedure requires 
the same precautions regarding stationarity of single particle densities as for their single- 
cell equivalents. 

Multi-fold factorial moments are a familiar tool in radio- and radar physics and in 
quantum optics There, they relate to simultaneous measurement of photo-electron 



counts detected in, say M, time-intervals, or in M space points, leading to a joint prob- 
ability distribution PM{ni,n 2 , . . . ,um)- The importance of multi-fold moments derives 
from the fact that, e.g. in the simplest case of two cells, Fu = (n m n m /) is directly re- 
lated to the auto-correlation function of the radiation field and obeys, for small cells, 
the Siegert-relation [|JT|, whatever the statistical properties of the field. The higher-order 



moments are sensitive to higher-order correlations and to the phase of the field. 

Factorial moments and factorial correlators are intimately related quantities. In terms 
of inclusive densities one has: 

^ pq = In dyi "' dyp In dVp+1 ' ' ' dVp+q Pp+q ( yi > ' ' ' ' yp] Vp+l} " ' ' yp+q ^ ( 2,76 ' ) 

where p p+q is the inclusive density of order p + q. The integrations are performed over two 
arbitrary (possibly overlapping) phase-space cells fli and fi 2 , separated by a "distance" 
D. 

It should be noted that the definition ( |2.76|) is more general than (|2.74j) . For Qi = Q 2 
or D = 0, ( 2.76 ) reduces to the correct definition of F 2 whereas ( |2.74j ) is, in this case, 



equal to (n ) and misses the so-called "shot-noise" term —{n). 



5 In one-dimensional rapidity space, D is denned as the distance between the centers of two rapidity 
intervals; in multidimensional phase space a proper metric must first be defined 
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Factorial moments and factorial correlators of the same order are thus seen to differ 
only in the choice of the integration domains. Note that for p ^ q, definition ( |2.76 ) is not 
symmetric in p and q and a symmetrized version is often used in experimental work: 

F$ = {F pq + F qp )/2. (2.77) 

From (|2.76|) follows that F\\ is directly derivable from measured two-particle correla- 
tion functions or from appropriate analytical parametrizations. Higher-order correlators 
involve higher-order density-functions which, in general, are unknown. 

We now turn to a discussion of multivariate factorial cumulants. For M non-overlapping 
cells, we introduce the M-variate multiplicity distribution Pjw(ni, . . . , um) and the corre- 
sponding moment- and factorial-moment generating functions: 



M(zi, ...,z M ) 
G(zi, . . . , zm) 



E £■ 

ni=0 ri2=0 
oo oo 

E E 

ni=0 ri2=0 



n M =0 



(2.78) 



Y, (l + zx^-'-il + ZM)^ P M (m,...,n M ) ,(2.79) 

n M =0 



from which the M-variate moments are easily obtained by differentiation: 

'_d_y ( a v m 
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M(z u ...,z M ) 



zi=—zm=0 



9M 



G(*i, 



(2.80) 
(2.81) 



Zl=—ZM=0 



The multivariate (ordinary) cumulants n qi ... qM and multivariate factorial cumulants 
f qi ...q M are likewise obtained by replacing M(.) and G(.) in ( |2.80[ ) and ( |2.81| ) by their 
respective natural logarithms |j32"| . The same expressions serve to extend the relations 
between univariate moments and cumulants to their multivariate counterparts. 

For M = 2 and non-overlapping cells, one has the identity [cfr. fl2.54 )]: 



EE 

1=0 m=0 



{z l ) l {z 2 ) m 



lm 



ex P (E E Um} flm) 



where F< 



oo 



Z=0 m=0 

1 and /oo is defined equal to zero. It follows that^ 



(2.82) 



Fn = /11 + /01/10, (2.83) 

F12 = /i 2 + /oi/20 + 2/io/n + /oi/i 2 , (2-84) 

F13 — /13 + /01/30 + 3/11/20 + 3/01/10/20 + 3/io/i2 + 3/ 10 /n + /01/10) (2.85) 
-£22 = /22 + 2/10/21 + /02/20 + /01/20 + 2/01/12 + 2/fj + 4/oi/lo/ll 

+ / 02 /i 2 o + /01/10 • (2-86) 

6 See also || 
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Similarly, expanding the logarithm in 

» » {- Zl )\- Z2 r . . ^ (-z 1 y(-z 2 r g, \ r9S7 s 

E E jf^j ^ = lo S ( E E ^ ^Jml, (2.87) 

i=0m=0 " b - \l=0m=0 L ' '"" I 

in powers of s and t and identifying coefficients, the reverse relations follow: 

fu = F±i — F iF w , (2.88) 

fi2 = F12 — F iF 2 o — 2Fio-Fn + 2F iii , (2.89) 

fi3 — F 13 — F01-F30 — 3FnF 2 o + 6F01F10-F20 — 3F 10 F 12 + 6F^ F U — 6F Q iFf , (2.90) 
/ 2 2 = F22 — 2F10-F21 — -^02-^20 + 2-F 01 F 2 o — 2F01-F12 — ^Fi X 

+ 8F 01 F 10 F 11 +2F 02 F? -6F* 1 F? . (2.91) 

The quantities F 0i , F i0 , f 0i and /jo are equal to the single-cell factorial moments 
and factorial cumulants, respectively. Expressions for Fji {fji) are obtained from the 
corresponding expression for Fij (fy) by permutation of the subscripts. By definition, 
foi = Fq\ and / i is equal to the average multiplicity in cell 2. 

It may be noted that the bivariate relations reduce to the univariate ones ( 2.56Q by 
simply amalgamating the indices. For example, from 

F12 = fu + /01/20 + 2/io/n + /01/10, (2.92) 

one recovers, by summing the indices 

Fs = fs + 3/i/2 + A 3 - (2.93) 

It is shown in Q (Sect. 13.12) that the above relations, while seemingly complex, have 
in fact a surprisingly elegant structure, rooted in simple algebraic properties of completely 
symmetric functions. Further discussion on this point and other useful properties may be 
found in [R2 . 



Extensions to more than two cells is straightforward, in principle, but involves tedious 
algebra. 



2.2 Poisson- noise suppression 

To detect dynamical fluctuations in the density of particles produced in a high-energy 
collision, a way has to be devised to eliminate, or to reduce as much as possible, the sta- 
tistical fluctuations — noise — due to the finiteness of the number of particles in the counting 
cell(s). This requirement can to a large extent be satisfied by studying factorial moments 
and their multivariate counterparts. It forms the basis of the factorial moment technique, 



known in optics, but rediscovered in multi-hadron physics in [21,22]. The method rests on 
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the conjecture that the multi-cell multiplicity distribution Pm(wi, • • • , %) can be written 

as 

P M (ni, ■ - • , n M ) = / dpi . . . dp M P P (pi, Pm) \{ , — exp (-p m 5). (2.94) 

J m=l n m- 

The Poisson factors represent uncorrelated fluctuations of n m around the average p m 5 = 
(n m ) in m-th interval; 5 is here the size of the interval. This can also be written as: 

P M {n l ---n M ) = (\{ ^L_exp(-(n TO ))) p , (2.95) 

m=l lim - 

where the outer brackets mean that an average is taken over the probability distribution 
of the densities p m , which are subject only to dynamical fluctuations. If these are absent, 
P p (pi, . . . , Pm) is simply a product of (^-functions. 

The formulae ( |2.94| - [2~95| ) are formally identical to the expression for the multi-interval 
photo-electron counting probability distribution in quantum optics and based on the fa- 
mous Mandel formula [ |3"4"| , |3"5] ] . The latter relates the probability distribution of the number 
of detected photo- electrons to the statistical distribution of the e.m. field. 

In optics, p m has the meaning of a space- or time-integrated field intensity. The 
ensemble average is calculated from the field density matrix which describes its statistical 
properties. 

Equations ( p.94|) -( pl95|) express -Pm^I; ■ ■ ■ > u m) as a linear transformation of 
P p (pi, ■ ■ ■ , Pm) with a "Poisson kernel". This transformation is known as the 



'Poisson Transform" of P p 



The Poisson-transform of a single- variable function f(x) is the function f(n) (n inte- 
ger) defined by the linear transformation 

_ POO J,™ 

f(n)= / dxf(x)—e-*. (2.96) 
Jo n\ 

A trivial example is the function 5(x — p) whose transform is the Poisson probability 
distribution. The Bose-Einstein distribution 

/(") = (1 /^n+i (- = 0> !'•••), (2-97) 

is obtained as the Poisson-transform of the exponential function (l/p)exp(—x/p). 

For suitably-behaved functions, the inverse Poisson-transform exists. It is closely 
related to the Laplace-transform of f(x). Several practical methods have been developed 
to determine the function f(x) from its Poisson-transform. Besides methods based on 
series expansions, the inversion problem may be reduced to an inverse moment problem. 
This follows from the equality between the factorial moments of f(n) and the ordinary 
moments of f(x), as further discussed below. 
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A table of useful transforms for probability distributions and further mathematical 



properties can be found in pi ] 



From the basic Poisson transform equation (|2.95|) it is easily seen that the multi-fold 
factorial moment generating function has the simple form 

M 

G(z u ...,z M ) = (J] exp {z jPj 5)) p , (2.98) 
i=i 

where the statistical average is again taken over the ensemble of densities pi, ■ •• ,Pm, as 
indicated by the subscript. 

On the other hand, the (ordinary) moment generating function of the densities is given 

by: 

Q( Zl ,...,z M ) = J P p (p 1 ,...,p M y^ + -- +pM ^dp 1 ...dp M (2.99) 
W / p 



Comparing ( gg8J ) and ( |2~TUq ), it follows that: 



G(zi, ...,z M ) = Q(8piZ!, p M z M S) . (2.101) 

This equation implies that the normalized multi-variate factorial moments of the multi- 
plicity distribution 

K nu = i \ Fqi "T \ ( 2 - 102 ) 

qi " qM (m)ft . . . {n M )iM ^ 1 

are equal to the normalized multivariate (ordinary) moments of the relative density fluc- 



tuation Pm/(pm}- This is the "noise-suppression" theorem |^T],^j. It assumes that the 



noise is Poissonian (cfr. fl2.94p ) and that the number of counts in all intervals (the total 



multiplicity) is unrestricted!. 

The property of Poisson-noise suppression has made measurement of factorial mo- 
ments a standard technique, e.g. in quantum optics, to study the statistical properties 
of arbitrary electromagnetic fields from photon-counting distributions. Their utility was 
first explicitly recognized, for the single time-interval case, in |J7| and || and later gen- 
eralized to the multivariate case in [32]. The authors of || further stress the advantages 



of factorial cumulants compared to factorial moments, since the former measure gen- 
uine correlation patterns, whereas the latter contain additional large combinatorial terms 
which may mask the underlying dynamical correlations (however, see the discussion in 
Sect. 2.1.2). 



7 If the sum over all intervals of the number of counts is fixed, a slightly more complicated relation 
can be obtained if the noise has a Bernoulli (multinomial) distribution |Bj| 
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Multivariate factorial cumulants are derived from the (natural) logarithm of the facto- 
rial moment generating function. Taking logarithms of both sides of (|2.101|) , one finds that 
the multivariate normalized factorial cumulants of the counting distribution are equal to 
the multivariate normalized ordinary cumulants of the densities [p5]. This relation, there- 
fore, extends the noise-suppression theorem to cumulants. This property is exploited in 
many fields from quantum optics |32| to radar-physics and astrophysics (see e.g. 



2.3 Sum-rules 

In an interesting a-model analysis of factorial correlators |pS|| , scaling relations are de- 
rived between single- variate and 2-variate factorial moments which are independent of the 
dimension of the phase space. The result is stated as follows: If a correlator Fn(D,5) is 
effectively independent of 5 in a range 5 < D < 5$, then 

F n (D) = 2F 2 (2D) - F 2 {D). (2.103) 

Here, 5 is the interval size and D the distance between the intervals. 

Similar types of relations — or sum-rules — are well-known in optics since the early 
1970's. They are exploited in so-called Multi-Cathode and Multiple-Aperture Single- 
Cathode (MASC) photo-electron counting experiments (see e.g. 0,f|l! and refs. therein). 

Consider again the multivariate multiplicity distribution Pm(«i, • • • , Um) giving the 
joint probability for the occurrence of n\ particles in a cell fli, . . . , n M particles in cell 
Qm, with Qi fl Qj = 0, Wi,j and % ^ j. Let n be the number of particles counted in the 
union of the M cells, 

M 

n=Y. n m ■ (2.104) 

m=l 

The probability distribution of n is given by 

n n 

p ( n ) = " ' X] PM(n u . . . , n M )Sn,n 1 +-+n M ■ (2.105) 
ni=0 njvf=0 

Define the single- variate factorial moment generating function 

oo 

^) = ^(l + 2 ) n ?( n ). (2.106) 

n=0 

The function g(z) can be expressed in terms of the multivariate generating function ( |2.79| ) 

as: 

g(z) = G m (z 1 ,...,z m )\ zi=Z2= ... =zm=z . (2.107) 
Equation (|2.107|) allows to express factorial moments of n in terms of the multivariate 



factorial moments of {ni, . . . , Um}- Application of the Leibnitz rule 

k m = E — j — y — f f^)" m ■ ■ ■ (f] ak f M (z) 

ai la 2 ! . . . a k \ \dzj \^ Z J 
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to the function 

f(z) = f 1 (z)---f M (z) 

leads immediately to the relation 



p, = Y F^ M \ — -. (2.108) 

The summation is over all sets {cij} of non-negative integers such that 



M 

E a i = 



Formula ( 2.108|) may be looked upon as a generalization of the usual multinomial theorem 
for factorial moments'! 

Likewise, taking the natural logarithm of both sides of Q2.107 ), one obtains an identical 
relation as Q2.108|) among single-variate and multivariate factorial cumulants. 

As an example, for two rapidity bins (M = 2) of size 5 separated by a distance D, one 
finds: 



+ F$, (2-109) 

The factorial moments F 0i are determined from the single cell (marginal) counting 
distribution, whereas the univariate factorial moments F q are obtained from the sum of 
the counts in the two cells. 

The relations derived in |39|] follow immediately from (|2.109|) by considering two adja- 
cent cells and normalizing properly. Since the derivation of ( |2.108| ) is completely general, 
it obviously holds irrespective of the dimension of phase space. 



F 2 


p(2) 
— -^02 


+ 2F 1 ? ) + ^ 2) 


F 3 


p(2) 


+ 3(F 1 ( 2 2) + F 2 ( 1 2) 


F, 


- ^( 2 ) 

— r 04 


+ 4(F$ + F 3 (2) 



The relations (|2.109| ) are trivially extended to more than two cells. They allow to 
measure high-order correlators by varying the distances between the cells. In optics 
and radar-physics, they are typically used in determining spatial coherence properties of 
arbitrary e.m. fields. 



2.4 Scaling laws 

A major part of this paper is devoted to recent experimental and theoretical research on 
possible manifestations of scale-invariance in high-energy multiparticle production pro- 
cesses. This work centers around two basic inter-related notions: intermittency and frac- 
tality. A review of the experimental data accumulated over the last years will be given in 
Chapter 4. Theoretical work is discussed in Chapter 5. 



See also @ 
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In particle physics, intermittency is defined, in a strict sense, as the scale-invariance 
of factorial moments ( 2.68|) -( 2.70|) with respect to changes in the size of phase-space cells 



(or bins) say 8y, for small enough 5y: 

F q (5y) oc {6y)-+* (5y - 0). (2.110) 

The power <p q > is a constant at any given (positive integer) q and called "intermittency 
index" or "intermittency slope". The form of ( |2.110| ) strictly implies that the inclusive 
densities p q and the connected correlation functions C q become singular in the limit of 
infinitesimal separation (8y — *> 0) in momentum space. 

Inspired by the theory of multifractals, scaling behaviour of the G-moments (|2.73|) has 
also been looked for in the form 

G q (Sy) oc (5 V yM (fy-0). (2.111) 

To describe the inter-relation of the two proposals, we briefly discuss the formalism of 
fractals. 

Power-law dependence is typical for fractals [f|2], i.e. for self-similar objects with a 



non- integer dimension. These range from purely mathematical ones (the Cantor set, the 
Koch curve, the Serpinsky gasket etc.) to real objects of nature (coast-lines, clouds, lungs, 
polymers etc). For reviews see [p5|-|4"5]] . 



The fractal dimension D F is defined as the exponent which provides a finite limit 

< limiV(e)e DF < oo (2.112) 



for the product of t° ¥ and the minimal number of hypercubes N(e) of linear size I = e 
(Kolmogorov definition) or I < e (Hausdorff definition) covering the object when e — > 0. 

To a physicist, the definition becomes more transparent if one considers the relation 
between the size / of an object and its mass M as a scaling law: 

Moc/ Df . (2.113) 

For usual objects Dp coincides with the topological dimension (for a line Dp = 1, for a 
square D F = 2 and so on). The condition e — > means in practice that such a law should 
hold in some interval of "rather small" e-values. 

The probability Pi(l) to be in a hypercube iVj(Z) is proportional to l Dp at small /. 
Therefore, for a fractal the mean value of the g-th order (ordinary) moment is given by 

(p?(Z)> oc l qDF {D F = const) . (2.114) 

Multifractals generalize the notion of fractals, since for these holds 

£p? (0 = (pf _1 (0) oc , (2.ii5) 
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where 

r(q) = (q-l)D q . (2.116) 

The D q are called the Renyi dimensions f46[ , |27| and depend on q (generally, for multifrac- 
tals they are decreasing functions of q). 

Sometimes it is more convenient to characterize multifractals by spectral properties, 
rather than by their dimensions. 

Let us group all the boxes with a singularity a (pi(l) ~ l a , I — > 0) into a subset S(a), 
where a is called the local mass dimension. The number of boxes dN a (l) needed to cover 
S(a) is 

dN a (l) = dp(a)/" /(a) , (2.117) 

where f(a) is the fractal dimension of the set S(a) related to the Renyi dimension. For 
the sum of moments one obtains: 

Ni(l) 



Ep?(0« U P {a)l aq - m • (2.H8) 

i=l J 



Dg = T mm a (aq - /(a)) = (aq - /(a)) (2.119) 



From ( |2.118| ), one gets by the saddle-point method: 

— -min a (aq - /(a)) = — 1 
q-1 q 

with a defined as 

rl f 

q{6c). (2.120) 



d/ 
da 



(a=a) 

The notion of Renyi dimensions D q generalizes the notion of fractal dimension Dq = 
Dp, information dimension D\ and correlation dimension D2 = v. A Renyi dimension, 
therefore, is often called a generalized dimension. 

The difference between the usual topological dimension D (i.e. the support dimension) 
and the Renyi dimension is called the anomalous dimension (or codimension) 

d q = D-D q . (2.121) 

The multifractal method is a widely used tool in many branches of physics and science in 
general (cfr. Q||,|7|). 

A direct relation may be established between the exponents of factorial and generalized 
moments at comparatively low values of q, much smaller than effective multiplicities 
contributing to the sum: 

<J> q + r(q) = (q-l)D. (2.122) 
Then the exponents are related to Renyi dimension and to codimension as 

r(q) = (q-l)D q (2.123) 
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and 

<P q = (q - l)d q . (2.124) 
According to the general theory |48| , ^9| , there exists "a class of multifractals exhibiting 



universal properties" . They are called universal multifractals and are classified by a Levy 
index < < 2 which allows the codimension to be expressed as 

d Q = —^—-^^- (Ci = const). (2.125) 
\i — 1 q — 1 

The Levy index is also known as the degree of multifractality (/U = for mono-fractals). 
Values fi < 1 correspond to so-called "calm" singularities, values fi > 1 correspond to 
"wild" singularities. 

One can proceed further and try to analyse experimental data at two different levels 
of bin-splitting. For that purpose, it was recently suggested to study Double Trace 

Moments (DTM). The procedure is, first, to sum up z/-th-order moments of multiplicity 
distributions at some bin-splitting level O within bins belonging to a single bin of one of 
the previous steps (having bins of size A) and then to calculate their g-th moments at 
that level 

Tr, oc £(£<)« oc A~ K ^-\ (2.126) 
a e 

It is claimed 0] that "the DTM-technique provides a robust estimate of \x and C\ for 



universal multifractals. According to the theory of universal multifractals [5(]], one 
should observe the following factorizable behaviour of "double" exponents K (q, v): 

K{q,v)=v*K{q,l) , (2.127) 

where fi is the same Levy index as in ( |2.125| ). 

Experimental results on multifractals and generalized multifractals, as well as some 
theoretical implications are discussed in Subsect. 4.7.7 and in Chapter 5. 



2.5 Bunching-parameter approach 

A simple mathematical tool alternative to the normalized factorial moments (2.68-2.70) 
is the bunching-parameter approach, suggested for high energy applications in ||52| . In 
order to reveal spiky structure of the events, it is only necessary to study the behaviour 
of the probability distribution near the multiplicity n = q by means of the "bunching 
parameters" 

Vq{oy) = 7 p2 77 v , q>l . 2.128 

As is the case for the normalized factorial moments, the bunching parameters r] q are 
independent of 5y if there are no dynamical fluctuations. For example, rj q = 1 for all q 
for the case of a Poissonian probability distribution. 
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As the F q (5y), the P q (5y) can be averaged over a number M of bins. Assuming 
approximate proportionality of n m and Sy at 5y — > and Po(Sy) ^ 1 for <5?/ ^ 0, one 
obtains 



or 



rj 2 (5y) oc (<fy) A 

77 ? (5y) ex (fy)-"" (2.130) 



with 



(3 q = d q (q-l)+d q _ 2 (q-3)-2d q „ l (q-2) , q>2. (2.131) 
Expressing d q in terms of the Levy-law approximation (2.125), 

ft = ^ + ( ,-y-2 (9 -iy (2132) 

In case of monofractal behaviour (/x = 0), f3 q = for q > 2. In the limit of the log-normal 
approximation (/i = 2), on the other hand, /3 9 = c? 2 and all bunching parameters follow 
the same power-law. 

The Levy-law approximation allows a simple description of multifractal properties of 
random cascade models using only one free parameter /j. In the bunching-parameter 
approach, one can make an approximation of the high-order bunching parameters to 
obtain a simple lin ear expression for the anomalous fractal dimensions d q , still maintaining 
the number of free parameters at one. 

Assuming that high-order bunching parameters can be expressed in terms of the 
second-order one as 

Vq (5y) = [r) 2 (5y)Y , q>2 , (2.133) 

the linear expression becomes 

d q = d 2 (l - r) + d 2 r| . (2.134) 

The use of bunching parameters is interesting, because it gives a general answer to 
the problem of finding a multiplicity distribution leading to intermittency: according 
to (|2.128| ), any multiplicity distribution can be expressed as 

pjm = p.M [PlW/ T oW18 fewr 1 -', *>i. (2.135) 

1- 1=2 



The possible forms of multiplicity distributions with multifractal behaviour of d q ( p. 121 



are discussed in [53,p4l 
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2.6 The wavelet transform 



An increase of factorial and cumulant moments with decreasing bin sizes reflects a widen- 
ing of a multiplicity distribution, i.e. an increase of multiplicity fluctuations in individual 
events. This phenomenon can be studied by other methods, as well. In particular, the 
so-called wavelet transform seems to be suited for that purpose. 

The wavelet transform is of particular importance in pattern recognition. This is a 
more general problem than the fluctuation study itself, since it involves the analysis of 
individual event shapes, not only the event ensemble, and may become of interest in the 
analysis of very-high multiplicity events. 



It is shown |55j that, for pattern recognition, the wavelet transform is about two orders 
of magnitude more efficient than ordinary Fourier analysis. 

An application of wavelets to multiparticle production processes has been proposed 
in ||56|1 . The main principle of the wavelet transform is to study the dependence of fluc- 
tuations on the phase-space bin size by the so-called difference method. One considers 
the difference between the histogram of an individual event at a definite resolution to the 
corresponding histogram at a (e.g., twice) finer resolution. Proceeding step by step, one 
is able to restore the whole pattern of fluctuations. 

Let us explain how this procedure can be applied to an individual event. We consider 
the one-dimensional projection of the event onto the rapidity interval Ay. Any n-particle 
event can be represented by the histograms of particle densities p = dn/dy at various 
resolutions. The simplest information is obtained from the value of the average density 
(p) = n/AY. To consider the forward-backward correlations, one splits the rapidity 
interval Ay into two equal parts and gets the forward and backward average densities 
(Pf,b) = 2 n f,b/ Ay, where rif^ are the forward (backward) multiplicities with rif + rib = n. 
Proceeding further to the J-th step, we approximate the event in terms of the histogram 
with 2 J bins. 

Let us construct now the difference of the two histograms described above. Namely, 
we subtract the average density from the forward-backward histogram and get another 
histogram with positive ordinate at one side and negative at the other, demonstrating the 
forward-backward fluctuations in the event. 

Splitting the forward and backward regions further into equal halves, one gets the 
histogram at J = 2. Its difference from the forward-backward histogram at J = 1 
reveals the fluctuations at finer resolution. Iterating to higher values of J, one studies 
how fluctuations evolve at ever finer resolution. The set of difference histograms is called 
the wavelet transform of the event. The above procedure corresponds to the so-called 
Haar-wavelet transform. Those interested in mathematical details are referred to 



The wavelet transform provides direct information on the evolution of fluctuations at 
different scales, i.e. on the dynamics of individual high-multiplicity events revealing their 
clustering (and subclustering) structure. A generalization to factorial (and cumulant) 
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wavelets is possible ||56|| . The simplest cascade models show such remarkable properties 



of wavelet transforms |56 as (quasi) diagonalization of their correlation density matrices, 
scaling exponents etc. It is interesting to note that the equations for the generating 
functions of wavelet transforms [56| look very similar to the "gain-loss" equations (in 



particular, to QCD equations) discussed at the end of Chapter 5. All those features are 
yet to be studied. 

The very first application to experimental data is presented in 
spectra of JACEE events are studied. 
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Chapter 3 

Experimental survey on correlations 



In this chapter, we review experimental results on "classical" correlations, a subject with 
a long history in particle physics. It was instrumental in establishing fundamental con- 
cepts of hadrodynamics, such as short-range order, which are an essential ingredient of all 
popular Monte-Carlo models of hadronization. With the exception of Bose-Einstein inter- 
ferometry, the field lay dormant for several years, but was revived with the introduction 
of generalized concepts. The data cover a variety of multiparticle-production processes 
ranging from e + e~ annihilation to nucleus-nucleus collisions. 

In Chapter 4, we shall review material on factorial moments and related quantities, 
obtained since 1986. At that time, a pioneering suggestion was made to investigate the 
patterns of particle density fluctuations in multihadronic events: the intermittency idea. 
Measurement of factorial moments opened a way to establish possible scale-invariance 
and fractal behaviour in hadrodynamics. 

Interest in correlation functions received a vigorous boost when their intimate con- 
nection with factorial moments was realized (see Chapter 2). Both are now explored in 
parallel with novel techniques. These offer promising perspectives towards a long overdue 
unified approach to correlation phenomena, including Bose-Einstein interferometry. 

Another obviously related subject, the phenomenology of multiplicity distributions 



59j, is not explicitly covered here. Multiplicity distributions inspired many early ideas 
on scale-invariance and phase-transition analogies in multiparticle production, such as 
Koba-Nielsen-Olesen scaling |60] and the Feynman- Wilson liquid picture ||61|| . However, 



the major part of the data relate either to full phase space or to sizable portions of it. 
It remains an interesting task for the future to explain the "large-scale" properties of 
multiplicity distributions in terms of correlation function behaviour at "small distances" , 
the main subject of this paper. Of course, the factorial moments discussed in Chapter 4 are 
just another representation of multiplicity distributions and their increase with decreasing 
bin size reveals the evolution of the multiplicity distribution. 
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3.1 Rapidity correlations 



The study of correlation effects in particle production processes provides information on 
hadronic production dynamics beyond that obtained from single-particle inclusive spec- 
tra. Correlations in rapidity y, as defined in Sect. 2.1, have been studied in various ex- 
periments on e + e~, lepton-nucleon, hadron-hadron, hadron-nucleus and nucleus-nucleus 
collisions. Strong ^-correlations have been observed in all experiments in one form or 
another, depending on the specific form of the correlation function, type of interaction, 
kind of particles, the kinematic region under consideration, etc. The main conclusions 
were (for early reviews see Q ) : 



1. Two-particle correlations are strong at small interparticle rapidity-distances [j/x — 2/2 1 
(see Fig. 3.1). 

2. They strongly depend on the two-particle charge combination. 

Rapidity correlations are now being studied with renewed attention. One reason is 
that their structure at very small rapidity distances is directly related to self-similar 
particle-density fluctuations (intermittency) , a topic to be covered in Chapter 4. 



3.1.1 Correlations in hadron-hadron collisions 



In Fig. 3.2 the pseudo-rapidity correlation function 6*2(771, 772) as defined in ( p. 25 ) is given 



for r/i = 0, as a function of r/ 2 = Vi f° r t ne energy range between 63 and 900 GeV [|65| . 
Whereas C 2 (0,r/) depends on energy, the short-range correlation C*s defined in (|2.36| ) 
does not strongly depend on energy and has a full width of about 2 units in pseudo- 
rapidity. The function Cl is not a two-particle correlation, but derives from the difference 
in the single particle distribution function for different multiplicities. As can be seen in 
Fig. 3.2b, Cl is considerably wider than C$ and increases with energy (the 63 GeV data 
are from lp6|| ). 

In Fig. 3.3, the semi-inclusive correlation C^^r/i, 772) for pp collisions at 900 GeV |J7) 
is compared to the UA5 Cluster Monte Carlo (MC) GENCL |j8|, as well as to the 
FRITIOF 2 £j| and PYTHIA [7T| Monte Carlos, for charge multiplicity 34 < n < 38. 



The Cluster MC is designed to fit just these short-range correlations, but also FRITIOF 2 
is doing surprisingly well (see however Subsect. 4.4.4). 



At lower energy, the NA23 Collaboration [[TT| has studied the short-range correlation 
of charged particles in pp collisions of ^Js = 26 GeV in terms of -^2(2/1,2/2) defined in 
Q2.32| ). Only events with charge multiplicity n > 6 are used. The positive short-range 
correlations are in agreement with those found earlier at ^/s = 53 GeV |72|. 

The NA23 data are compared to single-string LUND |73| and to a two-chain Dual- 
Parton Model (DPM) JF4| in Fig. 3.4. The one-string model (without gluon radiation) 
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does not at all describe the short-range rapidity correlation in the data. The two-chain 
model does better, but remains unsatisfactory. Somewhat better but still insufficient 
agreement is obtained by renormalizing the MC events to the experimental multiplicity 
distribution (not shown). The effect of Bose-Einstein correlations in the (++) and ( — ) 
data is found to be insignificant, as may be expected for data integrated over transverse 
momentum px and azimuthal angle (p. Obviously, more chains, possibly with higher p^, 
are needed to explain short-range order with fragmentation models, even below ^/s pa 30 
GeV. 



NA22 results for C 2 (0,t/ 2 ) and C 2 (0,t/ 2 ) (Eqs. ^25), gjjp for vr+p and K+p collisions 



at y/s=22 GeV J75[ are compared with FRITIOF 2, a 2-string DPM and QGSM [[76 



predictions in Fig. 3.5a,b. FRITIOF and 2-string DPM largely underestimate the cor- 
relation. QGSM reproduces C 2 (0, y 2 ) very well and even overestimates C 2 + (0,y 2 ) an d 
C 2 (0,7/2)- K has been verified that the differences between QGSM and FRITIOF or 
DPM are not due to the different treatment of tensor mesons (only included in the latter 
two). 

In Fig. 3.5c, FRITIOF and QGSM are compared to the NA22 data in terms of the 
short-range contribution C$ (0,7/2)- The (H — ) short-range correlation is reproduced rea- 
sonably well by these models. For equal charges, however, the strong anti-correlation 
predicted by FRITIOF is not seen in the data. QGSM contains a small equal-charge cor- 
relation due to a cluster component, but still underestimates its size. Similar discrepancies 
are also observed in semi-inclusive (fixed multiplicity) data for each charge combination 
(not shown here). They are even larger than in the inclusive data, also in the QGSM 
model. 

From this brief survey, we conclude that in hadron-hadron collisions two-particle cor- 
relations are badly reproduced and generally underestimated in currently used models. 



3.1.2 Correlations in e + e and /i + p-collisions 



Fig. 3.6 shows K 2 ~ (7/1,7/2) and K 2 (7/1, 7/2) for muon-nucleon interactions at 280 GeV/c 
TTfl . A steep peak is seen at 7/1 = y 2 = for K 2 ~, with two shoulders along the diagonal 
2/i = l/2- On the other hand, K 2 is below for most of the distribution, but we shall 
see that the most impressive correlation is in fact coming from y\ pa y 2 , just for this case. 
As in hadron-hadron collisions, correlations are strong and depend on the two-particle 
charge combination. 

Fig. 3.7 shows K 2 (yx,y 2 ) in /i + p interactions at 280 GeV/c with y\ e [—0.5,0.5], the 
hadronic invariant mass W in the interval 13 < W < 20 GeV and for n > 3 WM, together 



with the NA22 non-single-diffractive M + p sample, n > 2 [75]. Correlations in /i + p seem 
smaller than in NA22, but one has to consider a possible energy dependence. Indeed, 
extrapolating from the energy dependence of K 2 (0, 0) published in |78| , one finds quite 
similar values for /i + p at 22 GeV and M + p in NA22. 
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In Fig. 3.8a,b we compare the function ^2(0, y) for the NA22 non-single- diffract ive 
M + p sample (charge multiplicity n > 2) [|75| with that for e + e~-annihilation at the same 

energy (i/i=22 GeV) J75|] . The values of if 2 (0, are larger for (++) pairs than for ( ) 

in meson-proton (M + p) reactions; for ( ) and (-| — ) pairs they agree with K 2 for e + e~ 

annihilation in the central region. 

A comparison of the correlation functions for e + e~-annihilation and non-single-diffract- 
ive M + p collisions throughout the full kinematic region with y\ e [—1,0] is shown in 
Fig. 3.8c for charged pairs. The e + e~ data are given at ^Js = 14 and 44 GeV [ ff"9[ . 
At y 2 = yi, the 22 GeV M + p correlation lies between the e + e~ results. The shape is, 
surprisingly, more symmetric than in e + e~. 

For /i + p [fr7],f78l and e + e~ collisions ||79"HsTll , the LUND-type Monte Carlo is reported 



to reproduce the majority of the experimental distributions. In jnj it is shown that 
this is mainly due to the inclusion of hard and soft gluon effects. However, important 
underestimates of -fT 2 (2/i, 2/2) are still observable, in particular in the central and current 
fragmentation regions. For e + e~ ||80|| , this is shown in Fig. 3.9, where -^2(2/1, 2/2) is com- 
pared to the LUND model (JETSET 7.2 PS) as a function of y\ — y 2 (dotted line), for the 
full sample (upper plots) and for a two-jet sample (lower plots). In all cases, the LUND 
model underestimates the correlation at y\ — y 2 — 0. In general, the disagreement be- 
comes smaller when Bose-Einstein correlations are included (full lines). The main feature 
to note is that correlations are much weaker in the two-jet sample than in the full sample. 
Furthermore, correlations are larger for y < (left plot), i.e. in the hemisphere opposite 
the most energetic jet, than for y > (right plot). These two observations, again, point 
to hard gluon radiation as the main source of two-particle correlation in e + e~ collisions. 

A systematic test of analytic QCD calculations and of QCD Monte-Carlo models 
for two-particle correlations has been performed by OPAL ||81|| . The authors study the 
function 

fl(£i,6) =#2(6, 6O + 1 (3.1) 

with £ = ln(l/xp), xp = 2p/E cm being the Feynman variable, i.e. the particle momentum 
p in the cms normalized to half the cms energy E cm . In Fig. 3.10, R is plotted as a function 
of (£1— £2) fo r (£i+£2) centered at the values 6, 7 and 8, respectively. Fig. 3.10a proves that 
a next-to- leading order calculation (full lines) is better than leading order (dashed), 
but still overestimates the overall level of the correlation for any reasonable value of A. 
Since the next-to-leading correction is large, still higher-order terms are needed. It is 
therefore likely that a satisfactory analytical treatment of correlations, even at the parton 
level, will not be obtained in the very near future. 

Higher-order effects are, in an average sense, included in the existing Monte-Carlo 
models. In Fig. 3.10b, the same data are compared to the coherent parton shower models 



JETSET PS |73f, HERWIG §3j and ARIADNE §4J. The latter gives an excellent fit to 
the data, JETSET lies slightly below (within uncertainty of parameters), but HERWIG 
considerably above. The agreement of JETSET could only slightly be improved by includ- 
ing Bose-Einstein correlations. As far as incoherent parton shower models are concerned, 
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none of the various versions of COJETS [|85| gives a particularly good representation of 
the correlation data. 

All the models were tuned on the OPAL data in terms of event shapes and generally 
describe single-particle distributions. It is clear that correlations allow better and more 
discriminative tests than more integrated quantities. 

We have mentioned the difficulties string-hadronization models experience in predict- 
ing like-sign correlations in hadron-hadron collisions. It is important to verify if the 
otherwise successful e + e~ models are also able to reproduce correlations between charge- 
separated systems such as (H — ) and (±±) particle pairs. 



3.1.3 Charge dependence 

How Cs and C$ depend on the charge of the pairs is shown in Fig. 3.11 for the combinations 



(++) and (H — ) in NA22 [73]. The short-range correlation is significantly larger 



for (-| — ) than for ( ) and (++) combinations. This is also seen in the EMC data \7%. 

Resonance production is a likely explanation of this difference. For like charges, a small 
enhancement is seen near y\ m y% m above a large negative background. This is possibly 
due to Bose-Einstein interference. 



3.1.4 Charge-multiplicity dependence 

The multiplicity dependence of 0^(0, y) for the (-1 — ) combination is shown in Fig. 3.12 
||75|| . Near the maximum at y = the correlation function is approximately Gaussian 
and narrows with increasing n. In Fig. 3.13a are presented the values of (7^(0, 0) as a 
function of n for three charge combinations. Within errors, C2 (0, 0) is independent of n, 

but consistently higher for (H — ) and ( ) than for (++)• The reason for the difference 

between ( ) and (++) probably lies in the positive charge of both beam and target. 

On the other hand, an increase of (\vi ~ Vzl) with l/(n — 1) is found [^(| when 



averaging over a region \r]\ < 2 (Fig. 3.13b). Since C2 becomes smaller when moving 
away from the center, and that may happen faster for higher than for lower n, this is not 
necessarily in contradiction with the data in Fig. 3.13a. 



3.1.5 Transverse momentum dependence 

The search for density fluctuations, described in later sections, has revealed the impor- 
tance of correlations in multi-dimensional phase space. It is, therefore, of interest to gain 
insight into the transverse momentum (px) dependence of rapidity correlations. Early re- 
sults on this topic can be found in JS7J. Recent data on -^(0, 2/2) for all particles and 
for particles with smaller or larger than 0.3 GeV/c, plotted in Fig. 3.14, indeed reveal 
a strong sensitivity to transverse momentum. The correlation function is largest, and 
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stronger peaked, near y 2 



effect was noted already in [87 



for pt < 0.3 GeV/c, in particular for ( )-pairs. A similar 

The data of Fig. 3.14 were fitted with the functions 



fi = cexp[-(y-y ) 2 /2a 2 ] (full line) 
f 2 = aexp(— b\y\) (dashed) , 



(3.2) 
(3.3) 



with c, yo, a, a and b as free parameters. Even though for low the data point at y 2 = 
lies systematically above the curve, K 2 (0,y 2 ) is well fitted by the Gaussian fi but not by 
the exponential f 2 , in this one-dimensional projection on rapidity. 

Changing to the variables Xi = (y 2 + yi)/2 and x 2 = (y 2 — yi)/2, a steepening is 
observed at small x 2 (not shown). For like-charge pairs, this becomes particularly sharp 
when the bin size is reduced to 5x 2 = 0.1. For the latter, C 2 (xi,x 2 = 0) increases and 
both a Gaussian and an exponential can fit the correlation function. 



3.1.6 Strange particles 

In string-fragmentation models, first-rank hadrons are formed from neighbouring quark- 
antiquark pairs tunnelling out of the vacuum. The hadronic final states, therefore, show 
short-range order due to local flavour conservation. Using stable mesons only, this charac- 
teristic property is difficult to study experimentally because of the large qq combinatorial 
background. What is needed is a flag identifying the qq pairs created together. A suitable 
choice is strangeness since the number of ss-pairs per event is small and the combinatorial 
background strongly reduced. 

Good strangeness identification is available for e + e~ annihilation in the TPC detector 
at yfs = 29 GeV ||89|| . This collaboration observes significant short-range K + K~ correla- 
tions in y, well reproduced by the LUND model and by the Webber QCD model. 

In hadron-hadron collisions, strange particle pairs have been studied by the NA23 
Collaboration |K|. The distribution in the rapidity difference Ay for two K°'s is given in 



Fig. 3.15a, for a K° and a A in Fig. 3.15c. The results are compared to the single-string 
LUND model. As is the case for non-strange particles, the model slightly underestimates 
the rapidity correlation. 



3.2 Azimuthal correlations 



In interactions of unpolarised particles, no distinguished direction exists in the plane 
transverse to the beam and the distribution in the azimuthal angle ip is uniform. Still, 
a two-particle correlation exists also in <p and is visible in the distribution W(Aip) of 
Atp = \(fi — (p 2 \, the azimuthal angle between two particles, Aip e (0, 7r). The azimuthal 
correlation may depend on the charge of the particles in the pair, on the rapidity distance 
Ay = \yi — 2/2 1 between these particles and on their transverse momentum. 
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The first experiments to extensively study two-particle correlations as a function of 
both rapidity and azimuthal angular separation [p6| , |9~l| already showed that the correlation 
at small rapidity distance is strongest when the two particles are produced in the same 
or opposite directions in transverse momentum (see Fig. 3.16). The correlation- length 
in rapidity is larger towards Aip = tt than towards Aip = 0. Furthermore, significant 
differences in the shape of the joint rapidity and azimuthal correlation functions have 



been observed for pairs of like and unlike pions |91 



In Fig. 3.17, the distribution W(A<p, Ay), normalized to unity, is shown as a function 
of Aip, for all charge combinations, in the intervals Ay < 1, 1 < Ay < 2 and 2 < 
Ay < 3 H§B |. A horizontal line at the average value 1/tt corresponds to a flat distribution 



in Aip. The distribution is influenced by conservation of transverse momentum, by the 
decay of resonances (mainly for unlike-sign particles) and by Bose-Einstein correlations 
(for like-sign particles). In all cases, W is larger than 1/tt for Aip > tt/2 and has a 

maximum at Aip = tt. Except for ( ) pairs at Ay < 1, the W function is smaller than 

1/tt for Aip < tt/2. Such a global anti-correlation follows from transverse momentum 
conservation. 

Model predictions are shown in Fig. 3.17 for FRITIOF 2 (dot-dashed), two-string 
DPM (full) and multi-string QGSM (dashed). The comparison with the data shows that 
it is much easier to account for azimuthal correlations at large than at small Ay. At 
small Ay the models differ from each other and from the experimental data. The QGSM 
shows somewhat better agreement with experiment than the other models. This is a 
consequence of the multi-string structure of QGSM, where strong azimuthal correlations 
in a single string are destroyed, with the result that the A^-dependence is weaker than 
in two-string models. 

Differences between experiment and all models exist at small Aip and Ay < 1, in 

particular for ( ) pairs. Bose-Einstein correlations, not included in the models, may 

explain this disagreement. The influence of Bose-Einstein correlation can also be observed 
in the (++) combination, but is smaller because of the influence of the (positive) beam 
particle. 

Azimuthal distributions are shown in Fig. 3.18 for particles with Ay < 1, for p^ < 0.30 
GeV/c and for p^ > 0.30 GeV/c, together with model calculations. A comparison of these 
figures reveals that azimuthal correlations have a strong pp-dependence. Large positive 
azimuthal correlations exist at small Aip and Ay < 1 for like-sign particles with small pr- 
As the transverse momentum of particles increases, the peak at Aip = tt becomes more 
pronounced. This is reproduced by the models and reflects momentum conservation. 



For AA pairs, an azimuthal correlation has been observed in MARK II at 29 GeV [[52 
Similar K + K~ correlations are seen in the exclusive hh final state K~p — > 
pK+K-K-TT+vr- at 32 GeV/c [(JJ. 

Azimuthal correlations between (H — ) and (++, ) charge combinations have been 



studied in /xp collisions W7\ for \Ay\ < 1 and \Ay\ > 1. The distribution W(Aip) is 
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described fairly well by the LUND model including primordial and gluons, except that 
for | Ay | < 1 it slightly underestimates the anti-correlation for (H — ) and overestimates it 
for (++, ). 

In the azimuthal correlation of K° pairs (Fig. 3.15b) and of K°A° (Fig. 3.15d) studied 
by NA23 |Xj , the data tend to show pairs of small Aip not present in low-pr LUND (solid 
line) . 

By the same collaboration, the azimuthal correlation is studied [71 in terms of the 
asymmetry parameter 



B = [N(A<p > tt/2) - N(Atp < r/2)]/N M (3.4) 

for hadron pairs with 

a) opposite charge (h + h~) 

b) equal charge (h + h + + h~h~) 

c) possibly opposite strangeness (A°1i + ,xa < —0.2) 

d) no opposite strangeness (A°h~,XA < —0.2), 

for Ay < 2 and for Ay > 2. No azimuthal correlation is seen for Ay > 2 in all cases and 
for Ay < 2 in case of no common qq pairs (h + h + + h~h~,A°h~). For h + h~ and A°h + , 
the parameter B is compared to low-pr LUND and DPM predictions in Table 3.1. 



Table 3.1 Asymmetry parameter B 







Experiment 


LUND 


DPM 


A°h^ 


~(Ay < 2) 


0.18 ±0.03 


0.30 ±0.01 


0.19 ±0.01 


h+h" 


-(Ay<2) 


0.066±0.003 


0.126±0.002 


0.106±0.002 



The parameter B is strongly overestimated in single-string low-px LUND and still too 
large in the two-string DPM. Furthermore, B increases with the sum of the transverse 
momenta (Fig. 3.19) but less strongly than in the models. 

The azimuthal correlation has also been studied for cc pairs in DD production. An 
asymmetry is indeed observed in 7r~p collisions at 360 GeV/c ||94|| . Also there, the LUND 
model overestimates the effect. 



As shown on 7r~N interactions at \fs = 26 GeV also NLO perturbative QCD 
calculations overestimate the azimuthal asymmetry for DD pairs (Fig. 3.20a). Agreement 
can be obtained with a model where a (Gaussian shaped) transverse component is 



added to the incoming parton momentum before performing the NLO perturbative QCD 
calculation (Fig. 3.20b). 
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3.3 Correlations on the parton level 



The OPAL collaboration |97j has compared hadronic azimuthal correlations to coherent 
and incoherent shower models (Fig. 3.21). The coherent models JETSET PS with angular 



ordering |7§, HERWIG |83[ and ARIADNE |8J describe the azimuthal correlations in 
hadronic Z° decays, but the incoherent models JETSET PS without angular ordering [73 
and COJETS g| fail for (p > tt/2. 

The hadronization of a quark-antiquark pair at high virtuality is currently thought to 
proceed via parton showering Q. QCD implies that this parton showering be coherent. 



The coherence can be incorporated into Monte-Carlo programs as angular ordering [[99 
whereby for each successive branching the gluon is emitted at a smaller angle. 

Furthermore, the idea of local parton-hadron duality (LPHD) [|100|1 suggests that fea- 
tures at the parton level survive the fragmentation process. We can, therefore, expect 
that the coherence of the parton radiation will be reflected in angular ordering of the 
observed particles. 

As a method particularly sensitive to angular ordering, particle-particle correlations 



PPC and their asymmetry PPCA ||101| , |102[| are examined in a way analogous to the study 



of energy-energy correlations [103] 



1 n 1 

ppc(x) = Tr-PEnfMx-Xv)) (3-5) 

X i<j n 



PPCA( X ) = PPC(180°- X )-PPC(x) , (3.6) 

where Xij is the full spatial angle between tracks i and j, ( ) is the average over all events 
in the sample, n is the number of charged tracks in an event, and Ax is the bin width. 
The function Su n {x ~ Xij) is 1 if Xij an d X are m the same bin and otherwise. 

At y/s = M z , the fraction of two-jet events is very high. For two-jet events, particles 
in different jets will in general be separated by an angle x greater than 90°. The PPC for 
X > 90° can, therefore, serve as an indication of what the PPC within a jet (x < 90°) would 
be in the absence of angular ordering. By forming the asymmetry, these 'uninteresting' 
correlations are effectively subtracted. The effects of angular ordering should, therefore, 
be more directly observable in the PPCA than in the PPC. Note, however, that the sign 



convention following ||103|| leads to a negative sign for a positive correlation. 

Figures 3.22a and b show the PPCA distribution of L3 data (corrected for detector 
effects ||102|| ) compared to coherent and incoherent Monte-Carlo models, respectively. 

In Fig. 3.22b we see that for x ~ 60° JETSET 7.3 PS without angular ordering 
(incoherent) disagrees strongly with the data, while being in fair agreement at larger 
values of x- COJETS is seen not to reproduce the data over the entire angular range. On 
the other hand, in Fig. 3.22a, the coherent Monte Carlo models, JETSET with angular 
ordering, HERWIG, and ARIADNE all reproduce the data reasonably well over the full 
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angular range. Note that the disagreement of the incoherent models can not be due to the 
Bose-Einstein effect. Turning this effect off in the non- angular ordered JETSET model 
does not raise but lower its PPCA points. So, the data from the L3 experiment strongly 
disfavour the incoherent models. 



3.4 Three-particle rapidity correlations 



Whether dynamical correlations exist beyond the two-particle correlations discussed so 
far is of crucial importance for much of the present search for scaling phenomena in 
multiparticle processes, a subject treated in Chapter 4. With conventional techniques, 
this question is not easy to answer and beyond the sensitivity of many experiments. 

Nevertheless, three-particle correlations in rapidity have been looked for in a number 
of experiments |63| , [75| , |f04Hf 6fl . The third order normalized factorial cumulant is defined 



as [cfr. (F32D 1: 

ts t \ i . . I da da da 

K3(yi,y2,V3)=Ca(y 1 ,y 2 ,y3)/ 3 , , , , (3.7) 



°inei d 2/i d 2/2 dy : 



3 



, . 1 d 3 a n 1 da da da 

^3(2/1,2/2,3/3) = ,11 + 2 3 A a a ( 3i 

<r inel dy x dy 2 dy 3 crf nel dyi dy 2 dy 3 

1 d 2 a da 1 d 2 a da 1 d 2 a da 



a Li dyidj/2 dy 3 of nel dy 2 dy 3 dj/i af nel dytdy 3 dy 2 
with a inc i = J]a n . 

n>8 

The C$(yi, y 2 , y 3 ) correlation function is determined as a sum of topological correlation 
functions: 

^(2/1,2/2,2/3) = Y, p «C { a\vi,y2,yz) , (3-9) 

n>8 

C 3 n \yi,y 2 ,y 3 ) = p { 3 n) (2/1,2/2, 2/3) - ^3 n) (2/i,2/2, 2/3) , (3.10) 



^3 (2/1,2/2,2/3) = P2 (2/1, 2/2) Pi (2/3) +P2 (2/2,2/3)Pi (2/1) +P2 (2/i,2/s)Pi (2/2) 

/a) , 



-2pS n) (2/i)M n) (2/ 2 )pS n) (2/3 



1 1 d 



a 



P3 n) (2/i,2/2,2/3) = ,/ oq , ~ , , (3.11) 

n(l, 2, 3) a n dy!dy 2 dy 3 

where n(l, 2, 3) is the mean number of three-particle combinations in events with charge 
multiplicity n. 
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The corresponding normalized function is defined as: 

#8(2/1,2/2,2/3) = Cs{yi,y2,y,)/Y,Pnpt\yi)tt\y2)pt\yz) ■ (3.12) 

n 



Because of small statistics, three-particle correlations were not observed in pp interac- 
tions at 200 GeV/c at FNAL |63|. In K~p interactions at 32 GeV/c fL04] , three-particle 



correlations were considered using £3(2/1,2/2,2/3) and Kg(yi,y 2 ,y 3 ). No positive short- 
range correlation effect was observed. Correlations in the form of K have been observed 
in the central region by the ISR experiment for n > 8 ||106| . 

Fig. 3.23 from NA22 shows K 3 (0,0,y) and Ks(0,0,y) for the combined M + p sample 
at 250 GeV/c |f75| . Also shown are the values of -^3(0, 0, y) obtained in pp-interactions at 
i/i=31-62 GeV [ |106|| (lines). Inclusive three-particle correlations K 3 (0,0,y) are indeed 
seen in the NA22 data. They are strongest when a third particle partially compensates 
the charge of a pair of identical particles. There are, however, no correlation effects visible 
in the function K$(0, 0, y). In FRITIOF and QGSM, three-particle rapidity correlations 
are absent in both -^3(0, 0, y) and -Ks(0, 0, y). 

Recently, a factorization of the normalized three-particle correlation function has been 
proposed | |107| - |109( | under the form of a "linked-pair" structure: 



#3(2/1,2/2,2/3) = K 2 (yi, 2/2)^2(2/2, 2/3) + #2(2/1, yz)K 2 {y 3 , 2/2)- (3.13) 

The comparison of the prediction of Q3.13T ) to the data is given in Table 3.2, for 
n > 2, at a resolution of 0.5 rapidity units. At this resolution, the linked-pair ansatz is in 
agreement with the measured three-particle correlation within two standard deviations. 
Note, that 2/-correlations are much stronger for low-px particles and that the linked-pair 
ansatz continues to hold. 

With the accuracy presently attainable for three-particle correlations, it is obvious 
that studies of still higher-order correlation functions require better methods. The most 
successful ones will be discussed in Chapter 4. 



Table 3.2 The 3-particle correlation function compared to the prediction from the 
linked-pair ansatz, for non-single diffractive data (n>2). 
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3.5 Summary and conclusions 

1. The main contributions to the correlation functions C2 and 63 come from the mixing 
of events with different multiplicity and different single-particle density, but some 
effect remains in the so-called short-range correlation part. 

2. 6*2(0,1/2) increases much faster with increasing energy than its short-range contri- 
bution. 

3. The short-range correlation is significantly larger for (H — ) than for the equal-charge 
combinations, and is positive over a wider rapidity range in 6*2(1/1,1/2 —Vi)- 

4. The correlation functions 0^(0, 1/2), contrary to 6 2 n ^(0, 1/2), are similar for different 
multiplicity n, except that C^" 1 " - becomes narrower with increasing n. 

5. In hadron-hadron collisions, the correlation functions depend strongly on transverse 
momentum and are largest for small-pp particles. Consequently, correlations are 
stronger in multi-dimensional phase space than in a lower dimensional projection, 
such as rapidity space. Further implications of this observation will be discussed in 
Sect. 4.3. 

6. In the central cm. region, and at comparable energy, the correlation strength ob- 
served in M + p collisions at y/s=22 GeV is of similar magnitude as in e + e~ collisions 
and as in /xp collisions, if the trend of the latter is extrapolated to W = 22 GeV. 
Model predictions for e + e~ and /ip interactions slightly underestimate the correla- 
tion strength but give, nonetheless, clear evidence that (hard) gluon effects are the 
main source of correlations in rapidity space. 

7. Combinatorial background can be suppressed by studying the correlation of strange 
particles. Data are scarce, but support the conclusions drawn from data on non- 
strange particles. 

8. The UA5 cluster Monte Carlo and FRITIOF describe 6^(7/1, r] 2 ) at CERN-Collider 
energies, at least in the charge multiplicity range 34 < n < 38. At lower energies 
(20 < y/s < 30 GeV), the single-chain LUND model shows a strong anti-correlation 
among like-charge particles. The two-string FRITIOF model and DPM predict 
negative values for 6 2 (yi,y 2 ) or K 2 (yi,y2) in the central region for like charges. 
They are positive but far below the data for unlike-charge pairs. QGSM reproduces 
6*2(1/1,1/2) and 62(3/1,2/2) for all charge combinations, but cannot account for the 
short-range part 63(2/1,1/2). 

9. Positive correlations are observed at large values of the azimuthal angle Aip, as 
expected from transverse momentum conservation. The correlations among like- 
charge particle pairs at small values of A<p and Ay, where Bose-Einstein effects 
should contribute, are significantly larger than predicted by FRITIOF 2, DPM and 
QGSM. The deviations are stronger for particles with small transverse momentum. 
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10. In general, short-range correlations in e + e~ annihilation are reasonably well de- 
scribed by the LUND- and Webber-type models. To the contrary, in models for 
hh collisions which contain only one or two strings without additional pr effects or 
gluons, correlations in rapidity are underestimated, those in azimuthal angle over- 
estimated. 

Models such as LUND and DPM are known to underestimate the height of the "sea- 
gull" wings (the particle average transverse momentum as a function of Feynman- 
x) ||110||, a signal of semi-hard interactions. The models neglect such processes. This 



may partially explain why the models fail in both instances. 

11. The distribution in the interparticle opening angle of e + e~ collisions at LEP favours 
models with coherent parton showering. 

12. Three-particle correlations are now observed in all charge combinations. They are 
particularly large for low-pp particles. Within two standard deviations, they satisfy 
the linked-pair ansatz. No short-range contribution K$ is observed in three-particle 
correlations. Other methods are needed to study higher-order correlations. 
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Chapter 4 

Multiplicity fluctuations and 
intermittency 



4.1 Prelude 



The study of fluctuations in particle physics already has a long history going back to 
early cosmic-ray observations. To our knowledge, Ludlam and Slansky ff f f[ | were the first 
to advocate analysis of event-to-event fluctuations in hadron-hadron collisions. Compar- 
ing rapidity distributions of single events with the sample averaged distribution, they 
put in evidence strong clustering effects in longitudinal phase-space, indicating "a re- 
markably structured phase-space density" |TL2| . Fluctuat ions in individual events were 
also considered in the context of Reggeon theory in the important paper establishing the 



AGK-cutting rules ||113| . 



Early evidence for large concentrations of the particle number in small rapidity regions 
for single events were reported in cosmic ray experiments ||114| - |116|| and in pN collisions 
at 200 GeV beam momentum ||1 1 7\\ . A further number of high density "spikes" in rapidity 
space have been reported during the last decade. Fig. 4.1a shows the notorious JACEE 
event ||118|| at a pseudo-rapidity resolution (binning) of Srj = 0.1. It has local fluctuations 
up to dn/Srj ~ 300 with a signal-to-background ratio of about 1:1. The NA22 event ||119 
of Fig. 4.1b contains a "spike" at a rapidity resolution 5y = 0.1 of dn/dy = 100, cor- 
responding to 60 times the average density in this experiment. UA5 [|120|| has reported 
"spikes" in dn/drj up to 30 (10 times average) as early as JACEE, but found these to 
be in agreement with a short-range cluster Monte Carlo. Also EMU-01 ||121|| sees events 
with dn/drj = 140 satisfactorily explained by FRITIOF. 

From an experimental point of view, there is little doubt that events with large local 
density fluctuations exist. The real question is whether these are of dynamical or merely 
statistical origin, whether the underlying probability density is continuous or intermittent. 

Early attempts were made to answer the question of non-statistical fluctuations em- 
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ploying transform techniques |[L22|| , but these were not followed up. The problem resur- 
faced in the work of Bialas and Peschanski [pi, 22|, who suggested that spikes could be a 
manifestation in hadron physics of "intermittency" , a phenomenon well-known in fluid- 
dynamics. The authors argued that if intermittency occurs in particle production, large 
density fluctuations are not only expected, but should also exhibit self-similarity with 
respect to the size of the phase-space volume. 

Ideas on self-similarity and fractals in jet physics had earlier been formulated in ||123|, 
124 ], rephrased in the language of QCD branching processes in [ 125|| and in a simplified 



form in ||126|| . For soft hadronic processes, fractals and self-similarity were first considered 
127|| and their quantitative measures in ||128| , |129| . 



in 



In multiparticle experiments, the number of hadrons produced in a single collision is 
small and subject to considerable "noise" . To exploit the techniques employed in complex 
system theory, a method must be devised to separate fluctuations of purely statistical 
origin, due to finite particle numbers, from the possibly self-similar fluctuations of the 
underlying particle densities. The latter are the quantities of physical interest. A solution, 
already used in optics and suggested for multiparticle production in PH, W^i consists in 
measuring suitably normalized factorial moments of the multiplicity distribution in a given 
phase-space volume. 



4.2 Normalized factorial moments 



4.2.1 The method 



The method proposed in |2l],|22| consists in measuring the dependence of the normalized 
factorial moments F q (6y) defined in ( p . 6 8| - p770|) as a function of the resolution 5y. For 
definiteness, 5y is supposed to be an interval in rapidity, but the method generalizes to 
arbitrary phase-space dimensions. 

In Sect. 2.2 we have pointed out that the scaled factorial moments enjoy the property 
of "noise-suppression". It is easily verified that this crucial property does not apply 
to ordinary moments (n q )/(n) q (cfr. Sect 4.7 below). High-order moments further act 
as a filter and resolve the large n m tail of the multiplicity distribution. They are thus 
particularly sensitive to large density fluctuations at the various scales 5y used in the 
analysis. 



As proven in |2T],|22|], a "smooth" (rapidity) distribution, which does not show any 
fluctuations except for the statistical ones, has the property that F q (5y) is independent 
of the resolution 5y in the limit 5y — > 0. This follows directly from ( [2.100| ), if P p is a 
product of (^-functions in p m {m = 1, . . . , M) centered around (p m ). On the other hand, 
if dynamical fluctuations exist and P p is "intermittent" (i.e., regions of fluctuations exist 



at all scales of y), the F q obey the power law (|2.110| ). Equation (|2.110|) is a scaling law 
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since the ratio of the factorial moments at resolutions L and 



R = F g (£)/F q (L) = {LI If* 



(4.1) 



only depends on Ljl. 

As mentioned in Sect. 2.4 and Subsect. 5.2.2, the "intermittency indices" <p q (slopes in 
a double-log plot) are related [|130| , |131| , |28| to the anomalous dimensions d q = <fi q /(q — 1), 



a measure for the deviation from an integer dimension. 

We noted in Sect. 3.4 that the experimental study of correlations is difficult already 
for three particles. The close connection between correlations and factorial moments 
(Subsect. 2.1.4) offers a possibility to measure higher-order correlations with the factorial 
moment method at smaller distances than previously feasible. The method further relates 
possible scaling behaviour of such correlations to the physics of fractal objects. Despite the 
advantages, it should be remembered that reliable data can only be extracted if factorial 
moments are averaged over a large domain of phase space. This holds the danger of 
obscuring important dynamical effects. 

The definition of "intermittency" given in ( [2.1 10| ), has its origin in other disciplines^. 
It rests on a loose parallel between the high non- uniformity of the distribution of energy 
dissipation, for example, in turbulent intermittency and the occurrence of large "spikes" 
in hadronic multiparticle final states (Sect. 4.1). In the following we use the term "in- 
termittency" in a weaker sense, referring to the rise of factorial moments with increasing 
resolution but not necessarily according to a strict power law. 

The suggestion that normalized factorial moments of particle distributions might show 
power-law behaviour has spurred a vigorous experimental search for (more or less) linear 
dependence of In F q on — In Sy. Within a surprisingly short time (one-dimensional) anal- 



yses were performed for e+e" P3-p9fl, jxp |140| , vA jl41j , hh p2| - |148| , hA |149| - p4| 



and AA ||149| , |150| , |155|41~61"f collisions. With respect to the original objective, the early 
one-dimensional work has remained inconclusive, but valuable information and experi- 
ence was accumulated. Much more promising insight has come from studies in two- 
and three-dimensional phase space. This is discussed in Sect. 4.3. Further extensions 
of this approach, concentrating on improved integration methods and differential studies 
in Lorentz-invariant variables have lead to further clarification of the issues involved in 
intermittency. These very recent developments are presented in Sects. 4.8-4.10. 



4.2.2 Results on log-log plots (in one dimension) 

In this and the next few sections we review experimental results and model predictions 
obtained from one- dimensional studies. Due to the vast amount of data available, we limit 
ourselves to an illustration of the major characteristics of factorial moment behaviour in 
various processes and at various energies. 



For a masterly expose of this subject see 



132 1 
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In Fig. 4.2a, logF 5 is plotted |2T|,^2[ as a function of -log Srj (rj is the pseudorapidity) 
for the JACEE event. It is compared with an independent emission Monte-Carlo model 
tuned to reproduce the average 77 distribution of Fig. 4.1a and the global multiplicity 
distribution, but has no short-range correlations included. While the Monte-Carlo model 
indeed predicts constant F 5 , the JACEE event shows a first indication for a linear increase, 
i.e. a possible sign of intermittency. 

Further examples are given in Fig. 4.2b for KLM ||149|| , again showing an roughly linear 
increase for drj < 1 (— In £77 > 0) instead of the flat behaviour expected for independent 
emission, and in Figs. 4.2c and d for UAl ||143|| in terms of Srj and 5(f), respectively. 

Anomalous dimensions d q fitted over the range 0.1 < 8y(8rj) < 1.0 are compiled in 
Fig. 4.3 |[L62|| . They typically range from 0.01 to 0.1, which means that the fractal (Renyi) 
dimensions D q = 1 — d q are close to one. The d q are larger and grow faster with increasing 
order q in /ip and e + e~ (Fig. 4.3a) than in hh collisions (Fig. 4.3b) and are small and almost 
independent of q in heavy-ion collisions (Fig. 4.3c). For hh collisions, the g-dependence is 
considerably stronger for NA22 (y/s = 22 GeV) than for UAl (y/s = 630 GeV). 



4.2.3 Model predictions 
Hadron-hadron collisions 



A comparison to NA22 data on slopes <p q (Fig. 4.4a) shows [142] that intermittency 



is absent at y/s = 22 GeV in a two-chain DPM and underestimated by FRITIOF. In 



Fig. 4.4b, PYTHIA is seen to stay below the UAl data ||143|| , even after inclusion of Bose- 
Einstein interference for identical particles. The UA5 cluster Monte Carlo GENCL, able to 
reproduce conventional short-range correlations (at least in a certain range of multiplicities 
cfr. Fig. 3.3), follows the data down to a resolution of Srj ~ 0.3, but completely fails for 
smaller Srj. 

Also, a multi-chain version of DPM including mini-jet production has been compared 
to NA22 and UAl data. The slopes are found to be too small by at least a factor of 2 [ |163|| . 

With respect to intermittency analysis, the situation may improve with the introduc- 
tion of ECCO [|164j| , an eikonal cascade model based on geometrical branching, which now 
can account for strong fluctuations, in particular in higher dimensions (Sect. 4.3 below). 
However, the present version of ECCO is still less refined than the more conventional 
models with respect to other observables. 

The above examples show that present models for multiparticle production in hh 
collisions are unable to reproduce the magnitude and the growth of factorial moments 
with increasing resolution. From the discussion in Chapter 3, it is evident that model 
predictions for correlations in general are quite unreliable. The two-particle correlation 
function, measured by F 2 , also determines to a large extent the higher-order factorial 
moments (cfr. Eq. |2.72| ) because of the weakness of genuine high-order correlations. It 
is, therefore, mandatory to improve the models before evidence for "new physics" at very 
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small (rapidity) separation can be claimed. We return to this important question in later 
sections. 



hA and AA collisions 

The intermittency indices are much smaller in hA and AA collisions than in hh collisions, 
and the event samples are much smaller. Model comparisons are, therefore, less conclusive 



than in hh collisions. FRITIOF is found too low in NA22 jTBlj for 7T+/K+ on Al and Au 



at 250 GeV/c, in E802 J[B(J for central 16 OAl and 16 OCu at 14.6 A GeV/c, in WA80 [TUT 



for SS and Au at 200 A GeV/c, and in NA35 |L50[ for pAu, OAu, SAu and SS at 200 A 
GeV/c. In WA80 it is shown that rough agreement can be obtained by renormalization 
to the leftmost point of FRITIOF on the log-log plot (essentially the shape of the overall 
multiplicity distribution) to the data. NA35 shows that agreement can be obtained by 
adding Bose-Einstein interference for like-charged particles (for a detailed analysis of the 
influence of BE correlations see further below). 



Lepton-hadron collisions 



140 



are compared to what is expected from an extrapolation 

108| . At small 5y, the data are con- 



In Fig. 4.5a EMC data 

of conventional short- and long-range correlations 
sistently above these expectations. As Fig. 4.5b shows, the slopes 4> q in the same data 
are considerably larger than predicted by the Webber and LUND models. Similarly, 
Fig. 4.5c shows too low lnF 3 from LUND, not only for uNe but also for the "simpler" 
|T4I . 



z/D? interactions 



We tentatively conclude that also presently used lepton-hadron models as such are 
unable to reproduce the intermittency observed in this process. 



e + e annihilation 

The annihilation of e + e~ into hadrons is by far the best understood of all multihadron 
reactions. Creation of hadrons is traditionally pictured as a multi-step process comprising 
a "hard" parton evolution phase, described by perturbative QCD — the parton shower — 
and a non-perturbative colour-confining soft hadronization phase (Fig. 4.6). The former 
is a cascade process of nearly self-similar type, and is expected to show characteristics 
typical of a fractal object ||123| , |124| , |126|| . In fact, already in 1979, in a discussion of QCD 



jets, it was stated ||124|| that "the resulting picture of a jet is formally similar to that of 
certain mathematical objects, known as fractals, which look more and more irregular and 
complex as we look at them with a better and better resolution". The expectation is, 
therefore, that parton showers should exhibit intermittency at the parton level. However, 
this is not sufficient to guarantee "intermittency" at the hadron level. It is indeed difficult 
to imagine how the "re-shuffling" of the parton momenta during the hadronization phase 



52 



CHAPTER 4. MULTIPLICITY FLUCTUATIONS AND INTERMITTENCY 



with, e.g., the formation of hadronic resonances and their subsequent decay would preserve 
the (supposedly singular) nature of the correlations. A local parton-hadron duality type 
of explanation is not satisfactory either, since "it is merely a name for a mechanism that 
is not at all understood" [ |165|| . 

To describe the hadronization phase, all present Monte-Carlo codes rely in last instance 
on a large amount of e + e~ data at different energies and are carefully tuned to these. It 
came, therefore, as a surprise that a first (indirect) analysis [|133|1 of HRS results, shortly 



followed by TASSO data [134], revealed deviations from model predictions quite similar to 



those observed in lh and hh collisions (Figs. 4.7a,b). More recently, CELLO ]|135|| and, in 
particular, the LEP experiments ||136| - |138(| , claim "reasonable" agreement with the parton 
shower version of the LUND Monte Carlo (Figs. 4.7c,d). Nevertheless, new DELPHI data 
now show, with ten times larger statistics, significant deviations even with a "re-tuned" 
version of the Monte Carlo (Fig. 4.7d). 

The origin of intermittency in the models is not quite as clear as is often stated. Indeed, 



comparison of the factorial moments on parton and hadron level in Figs. 4.8a,b ||166 
shows that in (standard) JETSET the increase of lnF q at small 5y is not due to the 
parton shower, but to hadronization! Only if the parton shower is allowed to continue 
down to very low Ql values (Fig. 4.8c,d for Ql=0A GeV 2 ), implying local hadron-parton 
duality, is intermittency becoming visible also at the parton level. It has been verified 
that the influence of Qq is, of course, much less important at 1 TeV. 

On the other hand, intermittency seems to be fully developed on the parton level 
already at 91 GeV in the Webber model, and is in fact smeared out by hadronization [167]. 

The sensitivity to the cut-off in the perturbative QCD cascade and the role of hard 
and soft phases has also been discussed in terms of the dipole radiation model [ |168| J. 
Intermittency can be increased in the soft phase by an increase of the n/p ratio, also 
required from direct measurements by NA22 ||169| , EMC [|170|[ and in hA collisions \\L71\\ . 
The direct pions resolve the underlying parton structure better than the more massive 
resonances. From a tunnelling production mechanism, these pions are expected to have 
smaller than other particles, a property presently neglected in the MC programs. 
A Goldstone-like mechanism causing additional soft direct pion production at break-up 
points has recently been suggested by the LUND group ||172|| . 

For further progress, additional studies are needed. 

• One should identify the true causes of intermittency in present Monte-Carlo models, 
preferably on more sensitive distributions, such as those to be discussed below. This 
should reveal the influence of hard and soft gluon emission at high energies where parton 
showering is fully developed and dominates over the soft phase. 

• Intermittency is also particularly sensitive to the exact treatment of the soft phase. 
This phase can be studied with high statistics at lower energies where parton showering 
is less important. 
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4.2.4 A warning 

Before going into the necessary further detail, we should mention the influence of possible 
experimental biases. On purpose and by its very definition, the higher factorial moments 
are sensitive to a small number of events in the tail of the multiplicity distribution in 
small phase-space bins. 

Moments can be reduced by limited two-track resolution, by track losses from lim- 
ited acceptance or bad reconstruction, or simply due to truncation of the multiplicity 
distribution in a finite event sample. 

Moments can be increased due to double counting of tracks (track match failures), 
Dalitz decays and nearby 7-conversions or K°/A decays. A dangerous increase comes from 
the commonly used "horizontal" averaging, where a constant average (pseudo-) rapidity 
distribution is assumed over the range AY. Contrary to first belief, this problem is not 
completely solved by the correction method proposed in |E3] ! 



Further influence is to be expected from the choice of the sample, e.g., inelastic or 
non-diffractive, cuts on multiplicity, cuts on p^, all events or only those with n > n in 
AY, etc., the size and position of AY, the 8y region chosen for the fit and the correlation 
of errors. 

Many of these effects have been studied in a number of experiments and we refer to 
these and to ||173|, |174|| for more details. 



4.3 Higher dimensions 
4.3.1 The projection effect 

So far, we have discussed factorial moments derived from one- dimensional distributions in 
rapidity or pseudorapidity. The analysis can evidently be extended to other ID variables, 
such as the azimuthal angle ip in the plane perpendicular to the beam or event axis, 
or the particle transverse momentum (pr)- Given sufficient statistics, distributions can 
be analysed in two- and three-dimensional phase-space domains. Common choices are 
(Ay, Aip), (Ay, A lnp T ), (A^Alnpr) and (Ay, Aip, A lnp T ). 

Fig. 4.9a gives an example of ID-results from UAl ||143|| showing that intermittency is 
also present in ip. The intermittency effect is larger when two-dimensional cells (Ay, Ap) 
are studied than in ID (Fig. 4.9b,c). This is particularly pronounced in e + e~ annihilations 
(Figs. 4.10a,b), the measured slopes (j) q being about six times larger in 2D than in ID. 



These observations are now understood to imply that intermittency "lives in 3D" p4| , |175|| . 
Projection onto lower- dimensional subspaces dilutes the effect and leads to flattening 
of the factorial moments. This is most pleasantly demonstrated by the fact that one 
can enjoy a continuous (two-dimensional) shadow of a tree, in spite of the self-similar 
branching of this tree in three dimensions. 
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The projection-effect is convincingly illustrated in Fig. 4.10a,b. The lines in sub- fig. a) 
are fits by a 2D a-model; the curves in sub-fig. b) are the projections onto rapidity-space 
and show considerably less increase and even a flattening for by —>■ (note the difference 
in scale). Nevertheless, Fig. 4.10a still shows saturation of Fi at large M even in 2D, an 
indication that an analysis in three dimensions may be required. 



4.3.2 Transformed momentum space 

To study intermittency in three-dimensional phase space, one faces the additional difficulty 
that the particle density is all but uniform in the usual single-particle variables y, <p and 
Pt- The distribution in is in fact falling exponentially. Uniformity of the density is, 
however, an explicit assumption in the derivation of the power law Q2.110] ). Violation of 
this condition renders an intermittency analysis useless. 

To circumvent this problem, the authors of have proposed to use domains in 

a transformed momentum space with (almost) constant density. This is accomplished by 
a transformation of the original variables y, ip and lnpx to "cumulative" variables. Thus, 
for a single variable, say y, one defines the new variable X(y) as 

x(v) = JL-pMW u 2) 



For higher dimensions, it is assumed in [24| that the single particle density factorizes as 

Pi(v,<P,Pt) = P a (y)f>b(v)Pc(PT) • (4.3) 

Under this rather strong hypothesis, one can transform each of the three variables inde- 
pendently. The method proposed in [25[ does not assume factorization but is technically 



quite involved. In practice, the two techniques give satisfactorily similar results [|176| 



Data on F2 in various dimensions are shown in Fig. 4.11 for e + e ||136|| and hh col 



lisions fL42l, |14"3|1 . In all cases, the data behave more power-like in 2D than in ID. From 



Fig. 4.11a, it is also evident that JETSET PS remains in good agreement with e + e data 
in higher dimensions. 

At variance with power-law behaviour expected from intermittency, NA22 finds that 
the 3D factorial moments show an upward bending (Fig. 4.11c). This effect persists after 
exclusion of Dalitz decays and 7-conversions. A rise faster than power law is also observed 
in 3D for collisions of various projectiles with Au by NA35 (Fig. 4.12a) [|150|| . Following a 
suggestion in [|177|| , NA35 finds that the normalized factorial cumulant K% = F2 — 1 shows 



much better linearity in a log- log plot than F2 itself (Fig. 4.12b). 

This observation, in fact, furthers considerably our understanding of the intermittency 
phenomenon. In [ |177| , |178| the author has compared 3D data on F2 at ^fs ~ 20 GeV for 
/ip ||140|| , 7r/Kp |[L42|| , pAu, OAu and SAu ||150|| collisions using the parametrization 



F 2 = 1 + c(M 3 f 2 + d , (4.4) 
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where M 3 is the number of 3D phase-space cells. The second term in ( |4.4| ) is equal to 
The constant d accounts for long-range correlations, known to exist in hh collisions. 

The comparison of (|4.4|) to the data is shown in Fig. 4.13; the parameters are given 
in the figure caption. The parameter d is negligible for /xp and heavy ion collisions, but 
non-zero for meson-proton and pA collisions, in agreement with expectations. The most 
noteworthy result, however, concerns 02, which is seen to have a value in the range 0.4-0.5 
for all processes. This is remarkable in various respects. 

Firstly, if confirmed in further studies, and in particular for e + e~ annihilation, it 
strongly suggests that the resolution-dependence of F 2 exhibits a high degree of "uni- 
versality", is independent of specific details of the production process and thus reflects 
general features of hadronization dynamics. 

Secondly, such universality is at variance with the hitherto accepted idea that the 
factorial moments and the anomalous dimensions become the smaller the more complex 
the collision process, due to an increasing inter-mixing of production sources [|131| . 

Thirdly, if "universality" continues to hold in high energy e + e~ annihilation, one must 
revise the commonly expressed opinion that the perturbative parton evolution, and in 
particular hard-jet emission, is the primary cause of the rise of factorial moments at high 
resolution. Needless to say, it would be most interesting to verify systematically the 
universality conjecture in other reactions and for three-particle correlations. 

The experimental success of expression (fOJ) becomes quite intriguing when one real- 
izes that the volume 5 ~ M~ 3 of a phase-space cell (for sufficiently large M) is in fact 
related to the invariant mass Mj nv of the two-particle system or to Q 2 , the square of 
their four-momentum difference. The form ( |4.4|) implies that the two-particle correlation 
function behaves as a power law in M- mv or Q 2 . The data, therefore, seem to tell that an 
intermittency analysis should be performed in (Lorentz-invariant) multiparticle variables, 
rather than single-particle variables. This will be further discussed in Sect. 4.8-4.10. 



As mentioned in Subsect. 4.2.3 above, ECCO ||164|| has some success in describing 
the NA22 data on fluctuations in varying scales of resolution, in particular when the 
analysis is done in three dimensions (Fig. 4.14). The basis of this model is geometrical 
branching for soft production at low px- The geometrical aspect of hadrons, i.e. the 
fact that they are extended objects, puts the impact parameter R in a pre-eminent role. 
The fluctuation in R from event to event leads to fluctuations in pi and explains the 
non- vanishing intermittency in lnpx reported by NA22. The (stronger) intermittency in 
rapidity can be generated only with a singular splitting function for branching in rapidity 
space. Since there is no branching in (p in the model, intermittency is nearly non-existent 
in this variable. Still, the long-range correlation due to px conservation leads to a decrease 
of F q at large bin size, a feature also observed by NA22. 
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4.3.3 A generalized power law 

It has been pointed out [|179|| that the one- dimensional moments follow the generalized 
power law 

F q ex {g{5y)f« , (4.5) 



in multiplicative cascade models. In (|4.5| ), g(5y) is a general function of 5y. Expressing g 
in terms of F 2 , one finds the linear relation 

\nF q = c q + {<P q /<p 2 )lnF 2 , (4.6) 

from which the ratio of anomalous dimensions is directly obtained. This intriguing relation 
has successfully been confirmed by experiment, not only in one dimension, but up to 
3D |[24|| . Moreover, the ratios 4> q /4>2 are found to be largely independent of the dimension 
of phase space (Fig. 4.15a) and of the type of collision (Fig. 4.15b). 

The ratio of the anomalous dimensions d q (= <p q / (q — 1)) and d 2 are shown in Fig. 4.16b 
as a function of q. The g-dependence is claimed to be indicative of the mechanism causing 
intermittent behaviour. For a (multiplicative) cascade mechanism, in the log-normal 
approximation (long cascades), the moments satisfy the relation [pi|,|2"2"| 



dq _ 4> Q 1 _ q 
d 2 02 q - 1 2 



(4.7) 



However, the use of the Central Limit Theorem for a multiplicative process, such as in the 
a-model, is a very crude approximation ||180|| particularly in the tails. As argued in [p.81 



a better description might be obtained if the density probability distribution is assumed 
to be a log-Levy-stable distribution, characterized by a Levy-index \i. In that case 
generalizes to 

rl 1 «M _ n 

(4- 



d q _ 1 <f - q 



d 2 2"-2?-l 
For n = 2, the Gaussian case, ( |4.8p reduces to ( |4.7|) . 

The multifractal behaviour characterized by (|4.7|44.8| ) reduces to a mono-fractal be- 



haviour ranssi 



d f = 1 (4-9) 
d 2 

for jj, = 0, implying an order-independent anomalous dimension. This would happen if 
intermittency were be due to a second-order phase transition. Consequently, monofractal 
behaviour might be a signal for a quark-gluon-plasma phase transition. 

The data are best fitted with the Levy- law solution with \x = 1.6. This value is 
inconsistent with the Gaussian approximation, and also definitely higher than expected 
for a second-order phase transition. 

The validity of the dimension-independent generalized power behaviour has been ques- 
tioned in a recent NA22 analysis ||142|| shown in Fig. 4.16a. While a fit to the combined 
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data on all variables and dimensions (full circles), as well as a weighted average over all 



individual fits give \i values in rough agreement with those of p4 |, the 3D-data have /z > 2, 
not allowed in the sense of Levy-laws. 

Even larger values of /i, ranging from 3.2 to 3.5, have been found for /zp deep-inelastic 
scattering in ip.81 ]. According to |5(],|51[], this is evidence that the procedure to obtain 



the Levy-index is used outside its domain of validity. An allegedly more general method, 
based on Double Trace Moments (to be discussed in Subsect. 4.7.7) indeed yields /z-values 
within the mathematically allowed boundaries. However, we shall see that the latter 
method may be criticized on other grounds. A possible way out is self-affinity to be 
discussed in Subsect. 4.3.5, below. 



The linear d q /d2 behaviour in Fig. 4.16a) and b) gives some justification for ( 2.134 ). 
Fig. 4.16c) and d) show [|52|] the slope r of ( [2.134j ) for a number of experiments. All 
experiments, except perhaps SAg/Br, show multifractal behaviour (r > 0). 

Despite the confusion, it remains a noteworthy experimental fact that the facto- 
rial moments of different orders obey simple hierarchical relations of the type ft4.6|) . 
This means that correlation functions of different orders are not completely indepen- 
dent but are somehow interconnected. Such situations are commonly encountered in 
various branches of many-body physics (see e.g. [|107| , |109| , |184|| ), but a satisfactory link 
with particle phenomenology, let alone QCD, remains to be established. Nevertheless, 
on a simple example it was recently shown ||185|| that a linear relation between lni 7 ^ 
and lni 7 ^ can be obtained if the connected correlation functions are assumed to be of 
a factorized Mueller- Regge power-law form in two-particle invariant-masses squared Sy, 
i.e. C*3(l,2,3) oc (si2) 1_Ql (s23) 1 ~" 2 + cycl. perm.. Note that this Regge-form has the 
"linking" structure of (3.13). 



4.3.4 Thermal versus non-thermal phase transition 
Second order phase transition? 

A simple model that can provide some hint on the nature of a second-order phase tran- 
sition is the Ising model in 2D ||186||. Its intermittency behaviour has been studied both 



analytically and numerically ||182| , |187|1 . The anomalous dimension is found to be d q = l/£ 



independent of q. Based on that finding, it has been conjectured that intermittency may 
be monofractal if due to a QCD second-order phase transition ||183|| . However, as men- 



tioned in Subsect. 4.3.3 above, all types of interactions, including heavy-ion collisions, 
show multifractal behaviour. 

Of course, the Ising model is very simple and the above conjecture has little basis. 
In ||188|| , intermittency is, therefore, studied in the framework of the Ginzburg-Landau 
theory also used to describe the confinement of magnetic fields into fluxoids in a type II 
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superconductor. In the model, the anomalous dimension is not constant, but follows 

d 



d-2 



(9-1) 



v-l 



v = 1.304 , 



(4.10) 



with v being a universal quantity valid for all systems describable by the GL theory, 
independent of the underlying dimension or the parameters of the model. This is of 
particular importance for a QCD phase transition, since neither the transition temperature 
nor the other important parameters are known there. 

In quantum optics, 7 production at the threshold of lasing is discribable as a second- 
order phase transition. Indeed, a photo-count experiment | |189| | has verified ( |4. 1U[ ) to high 
precision. On the other hand, the current NA22 data on particle production in hadronic 
collisions give v = 1.45 ± 0.04 [ 190|| , heavy-ion experiments v = 1.55 ± 0.12 ||188|| and 
v = 1.459 ±0.021 |5§ . 

For a first-order phase transition, all d q are zero and no intermittency would be ob- 
However, it has been shown in [ |191|| that in a generalized GL model, a first- 



served 183 



order phase transition combined with the quantum optics analogy of lasing at threshold 
can lead to intermittency behaviour in some regions of the parameters, with approximately 
the same intermittency indices as a second-order phase transition. 



Non-thermal phase transition? 

Of course, the phase transition does not need to be thermal, i.e., the new phase need 
not be characterized by a thermodynamical behaviour. Such a transition could, e.g. take 



place during a parton-shower cascade and has been formulated in ||192|| for a number of 
"ultra-soft" phenomena, including intermittency. It leads to the co-existence of different 
phases, in analogy to different phases of the spin-glass systems. The examples of the 
JACEE event (Fig. 4.1a), which contains many "spikes" and "holes", and that of the 
NA22 event (Fig. 4.1b), which consists of just one spike, indicate that such a possibility 
may be more than just a speculation. 

The condition for the existence of such different phases of a self-similar cascade is that 
the function 

\ = (<P q + l)/q (4.11) 



has a minimum at some value q = q c (not necessarily an integer) [ |45] , |193| - |195| . The 



regions q < q c and q > q c are dominated by numerous small fluctuations and rare large 
fluctuations, respectively. In the terminology of |195||, the system resembles a mixture 



of a "liquid" of many small fluctuations and a "dust" of high density. We see either the 
liquid or the dust phase, depending on whether we probe the system by a moment of 
order q < q c or q > q c , respectively. 

In Fig. 4.17a, X q is compiled | |195| from KLM, EMC and NA22 as a function of the order 
q. The low pr NA22 data ||142|| (px < 0.15 GeV/c) indeed show a marked minimum with 



q c between 3 and 4, while the uncut data have not saturated at q < 5. Following ||195|| , the 
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X q behaviour has been studied by a number of heavy-ion experiments [ |153| , |154| , |156| , |159| . 
While a saturation, but no clear minimum is seen by experiments stopping their analysis 
at q = 5 or 6, a minimum is now observed for 4 < q c < 5 in central C-Cu collisions at 4.5 
A GeV/c, where the analysis is carried to q = 8 [|159|| (Fig. 4.17b). 



The observation of a minimum in the A„-distribution suggests a phase transition |E5 



193| . |194j| , but according to the interpretation ||195|| it is merely the "apparatus" changing 



from a sensitivity for the dominating small fluctuations at q < q c to an insensitivity for 
those at q > q c . The two phases could coexist without a transition being necessary. 

So, phase transition or not, two phases seem to coexist and it will be a challenge to 
find their physical interpretation in terms of the theory of strong interactions. 



4.3.5 Self-affinity 

Comparing log-log plots for one phase-space dimension, one notices that the liaF q saturate, 
but at different F q values for different variables y, ip or lnpx- The saturation in one dimen- 
sion can be explained as projection effect of a three-dimensional phenomenon. However, 
also in three dimensional analysis the power law (2.108) is not exact. In Fig. 4.11c, the 
3D data are seen to bend upward. 



It has been shown in | |196| | (see also [fL97[| ) that this can be understood by taking the 



anisotropy of occupied phase-space (longitudinal phase space ||198j| ) into account. In view 
of this phase-space anisotropy, also its partition should be anisotropic. In other words, 
the density fluctuation in phase space should be self-affine rather than self-similar [[L99|| . 



If the phase-space structure is indeed self-affine, it can be characterized by a parameter 



called roughness or Hurst exponent [199], defined as 



Hij = In Ai/ In Xj (0 < H, tj < 1) (4.12) 

with Aj (i = 1, 2, 3; Xi < A2 < A3) being the shrinkage ratios in the self-affine transfor- 
mations 

Sxi -> 5xi/Xi , (4.13) 

of the phase-space variables Xj. 

The Hurst exponents can be obtained from the experimentally observed saturation 
curves of the one-dimensional lnF q (5xi) distributions. Using the NA22 curves for y and 
hipr (Fig. 4.11c), a Hurst exponent of H ytPT = 0.516 ± 0.015 is obtained for these two 
variables, in agreement with self-affinity (H < 1) rather than self-similarity (H = 1). 

The upward bending for F q in the three-dimensional self-similar analysis is then easy 
to understand: Performing a self-similar analysis, phase space is not shrunk according 
to the self-affine dynamical fluctuation. So, the real dynamic fluctuation cannot be fully 
observed and the corresponding F q comes out smaller at intermediate scales. At very 
small bins, however, this difference between self-affine and self-similar space shrinkage 
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disappears and the F q values obtained approach each other. As a consequence, the slope 
on the log-log plot has to increase at small bin sizes and the self-similar analysis grants 
an upward bending if the underlying structure is self-affine (i.e. corresponds to a power 
law). 

On a self-affine Monte-Carlo branching model exactly reproducing the NA22 d q jdi 
values of Fig. 4.16a, this upward-bending effect is shown to cause the apparent violation 
of the Levy stability \i < 2 described in Subsect. 4.3.3 [ [200)1 ■ 



4.4 Dependences of the effect 



4.4.1 Charge dependence 

A mechanism known to cause correlations at small distances in phase space is Bose- 
Einstein interference between identical particles 
status of this field we refer to ||203| , |204| 



mm 



For reviews of the present 
From the outset it must be realized, however, 
that the conventional Gaussian- or exponential-type parametrizations of the Bose-Einstein 
effect lead to a saturation at Sy — > and not to the power law ( p.llOj )! 

In p02| | it is argued that the slopes should be roughly a factor 2 larger for identical 
particles than for all charges combined. The experimental situation is less than clear, 
in particular for ID analyses. Contrary to the prediction, TASSO and DELPHI see less 
intermittency for identical particles. EMC finds an enhanced effect for positive but not 
much for negative particles in a one-dimensional analysis, and very similar slopes in a 3D 
analysis. NA22 observes an enhancement for negatives, but not for positives. UAl sees 
no difference, whereas NA35 sees an increase. 

CELLO finds Bose-Einstein interference necessary to explain the residual difference 
between data and JETSET 7.2, but needs an un-physically large strength-parameter A to 
obtain agreement. In the DELPHI analysis, Bose-Einstein interference is insufficient to 
explain the difference between data and models, even with an un-physically large value 
of the coherence parameter A. 



Following a suggestion in [205] , higher-order Bose-Einstein correlations have been stud- 
NA22 [ED7J and DELPHI p08| . In this study, "correlation functions" 



ied by UAl 
of order q, 
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R q {Q%) = N q {Ql)/Nf G {Ql), 



(4.14) 

,5) 



are defined as ratios of the distribution of like-charged g-tuplets (q = 2, 3, . . 
and a distribution of reference (background) g-tuplets N q BG {Q 2 q7T ) obtained from random 
event-mixing. The variable Q 2 is defined as a sum over all permutations 



Ql = Ql2 + Qls + ---Q 



(9-1)3 



(4.15) 



of the squared four- momentum difference Qf, = — {pi — Pj) 2 of particles i and j. Note 
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that the functions ( f4.14| ) are normalized inclusive densities and not correlation functions 
in the proper sense (cfr. Sect. 4. 8). 

The UAl data are shown in Fig. 4.18. A good fit is obtained if in the expansion 
of R q {Q\ L) suggested in ||205|| , Gaussians (dashed curve) are replaced by exponentials in 



Q q7T (solid curve). Since low Qf- pairs are lost due to limited two-track resolution in the 
detector, the data at the smallest Q?- have to be regarded as a lower limit. A power law 
as expected from intermittency cannot be excluded. 

UAl has further studied the distributions R 2 for all-charged-(cc), (±±)- and (H — )- 
pairs as a function of Q 2 {= Q 27T ) (Fig. 4.19) These results have important implica- 



tions. The charge-dependence of "intermittency", controversial in single-particle variable 
analyses (see before), is now quite clear in invariant-mass variables (Q 2 = M 2 nv — 4m 2 ). 
The data for R^ ± (dashed) has a much stronger Q 2 -dependence than R 2 ~ and effectively 
determines the small-Q 2 behaviour of R 2 C . This is unambiguous evidence that intermit- 
tency at small Q 2 is predominantly due to like-sign particle correlations. It does not 
necessarily imply, however, that Bose-Einstein interference is the sole cause. 

In |77j it is shown on EMC data that, especially in 3D, F 2 deviates much more from 
LUND Model predictions than F 2 + ~ . The LUND Model version used does not include 
Bose-Einstein correlations. The deviation from the data is indicative for the importance 
of this effect. 

Bose-Einstein interference must thus play a significant role at least for small Q 2 . This 
seems in contradiction with claimed successes in e + e~ annihilation of parton shower Monte 
Carlos which neglect Bose-Einstein interference. 

Finally, we reiterate our remark that a "conventional" Bose-Einstein effect with expo- 
nential or Gaussian Q-dependence is incompatible with intermittent power-law behaviour. 
We return to this point in Subsect. 4.8.4. 



4.4.2 Transverse-momentum dependence 



An interesting question is whether semi-hard effects ||179|| , observed to play a role in the 
transverse- momentum behaviour even at NA22 energies ||110| |, or low-px effects | |192| , |209| 
are at the origin of intermittency. A first indication for the latter comes from the most 
prominent NA22 "spike" event | |119| |, where 5 out of 10 tracks in the spike have px < 0.15 
GeV/c. 

In Fig. 4.20a, NA22 data on lnF q versus — In 5y are given for particles with transverse 
momentum px below and above 0.15 GeV/c, and with px below and above 0.3 GeV/c. 
For particles with px below the cut (left), the F q exhibit a stronger (^-dependence than 
for particles with px above the cut (right). 

NA22 does not claim straight lines in Fig. 4.20a, but uses fits as an indicative measure 
of the increase of lnF q over the region 1 > 5y > 0.1. In the upper half of Fig. 4.20b, the 
fitted anomalous dimensions d q are compared to those obtained in the full px-range. The 
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restriction to particles with px < 0.15 or 0.30 GeV/c indeed leads to an increase of d q ; 
a decrease is observed for px > 0.15 or 0.30 GeV/c. This observation is confirmed by 
IHSC PBj 



FRITIOF predictions are given in the lower part of Fig. 4.20b, again for all tracks 
and for tracks with restricted px- It is known [|142j| that FRITIOF gives too small slopes 



for factorial moments integrated over px- Here, one notices that it also fails to reproduce 
their p T -dependence. 

UA1 has a bias against tracks with p T < 0.15 GeV/c, but gives the dependence of 2 
on the average transverse momentum px of the event (Fig. 4.20c) [ 210[| . The data show 



a remarkable decrease of 02 with increasing px and, after passing through a minimum at 
px ~ 0.5 GeV/c, a slight increase at higher px values. Lower px events correspond to soft 
processes, while higher px ones correspond to events with hard jet subprocesses. Both 
types of events have higher slopes 02 than their mixture at intermediate px values. (See 



further in pTTJ for a possible connection to the multiplicity dependence to be described 



in Subsect. 4.4.4 below.) 

Fig. 4.20c also contains the results obtained from Monte-Carlo events generated with 



PYTHIA 5.6 ||70|| . At low px values, the PYTHIA 2 values are strongly suppressed as 
compared to those of the data. 

We conclude that the intermittency observed in NA22 and UAl data is enhanced 
at low transverse momentum and is not dominated by semi-hard effects. Hard effects 
dominate in high energy e + e~ and lh collisions. Data on the px-dependence of factorial 
moments in these processes should help in clarifying the origin of intermittency. The 



effect of px-cuts on e + e data has been studied by DELPHI |[L36|| . One-dimensional data 
are shown in Fig. 4.21 and provide several important pieces of information: 

i) The log-log plot for low-px particles shows less saturation (i.e. stronger intermit- 
tency) than for larger px particles. So, intermittency is strongest in the px-region where 
hard gluon effects are weakest! 

ii) A discrepancy between data and models (only indicative in Fig. 4.7d above) is 
observed in the interval 0.255 < px < 0.532 GeV/c. This looks surprising at first, but we 
shall show in Subsect. 4.4.4 that the intermittency effect can be stronger for individual 
mechanisms than for a mixture. 

iii) The factorial moments are larger for px > 0.532 GeV/c than for px < 0.255 GeV/c, 
opposite to the trend of the NA22 data (Fig. 4.20a). Also this seems contradictory, but it 
should be realized that for NA22 transverse momentum refers to the beam axis, which is 
usually close to beam and target jet-axes. In the e + e~ analysis, px is calculated relative 
to the global event axis which differs from the direction of individual jets. 



4.4. DEPENDENCES OF THE EFFECT 



63 



4.4.3 Dependence on jet topology 

In their recent analysis, DELPHI [ |136|| selects 2-jet and 3-jet events using the JADE/EO 



invariant-mass algorithm |[212j1 , with resolution parameter values y cut = 0.04 and 0.01, 
and with additional cuts to clean the 2-jet and 3-jet sample. At large bin sizes, factorial 
moments rise faster with decreasing bin size (and are, therefore, larger) in 3-jet than in 
2-jet events. This is compatible with the (large bin-size) behaviour expected from hard 
gluons. At small bin sizes the increase is similar for 2-jet and 3-jet events. 

In 3-jet events, factorial moments were calculated for tracks belonging to jet 1, jet 2 
and jet 3 ordered in energy. The rapidity was defined with respect to the individual jet 
axis. As seen in Fig. 4.22, intermittency is weakest in jet 3 and strongest in jet 2. The 
deviation from JETSET is also strongest for jet 2. 



4.4.4 Multiplicity (density) dependence 



In general, a decrease of the intermittency indices <p q is found with increasing energy, 
in particular for hh, hA and AA collisions. As seen in Fig. 4.23a, a strong multiplicity 
dependence of the intermittency strength is observed for hh collisions by UA1 [ |143|| . 
The trend is opposite to the predictions of the models used by this collaboration. This 
decrease of the intermittency strength with increasing multiplicity is usually explained as 
a consequence of mixing of independent sources of particles ||131]| . The cross-over of data 
and FRITIOF at intermediate multiplicity explains the apparent success of FRITIOF in 
Fig. 3.3, for multiplicities close to 30 as being accidental. 

with 



Mixing of emission sources leads to a roughly linear decrease of the slopes Y , t 

. This is indeed 



increasing particle density (p) in rapidity ||108| , j213| , [214|| : <p q oc (p) 



observed by UAl p3| . Multiple emission sources are present in multi-chain Dual Parton 
models. The calculated slopes indeed depend linearly on multiplicity but are too small 
by a factor of two | 215|| . Similarly, the model studied in |[216|| with independent emission 
at fixed impact-parameter finds decreasing (p q with increasing multiplicity. 

Also here, a study of the multiplicity dependence in e + e~ data and JETSET allows 



interesting comparisons. In fact, the LEP results ||136|1 suggest little or no n-dependence, 
except for the lowest multiplicities, where the slope is largest and also the difference with 
JETSET PS is the largest. 

Fig. 4.23a helps in explaining why intermittency is so weak in heavy-ion collisions 
(cfr. Fig. 4.3): the density (and mixing of sources) is particularly high there. In Fig. 4.23b 
EMU01 ||155|| , therefore, compares 02 for NA22 (hp at 250 GeV) and heavy-ion collisions 
at similar beam momentum per nucleon, as a function of the particle density. Whereas 
slopes averaged over multiplicity are smaller for AA collisions than for NA22 in Fig. 4.3, 
at fixed (p) they are actually higher than expected from an extrapolation of hh collisions 
to high density and even grow with increasing size of the nuclei. This may be evidence 
for re-scattering (see ||141|| ) or another (collective) effect, but, as shown by HELIOS ||157 
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and recently confirmed by EMU-01 ||155| , one has to be very sure about the exclusion of 
7-conversions before drawing definite conclusions. 

We conclude this section with an additional warning. In Subsect. 4.3.2 we mentioned 
the Fialkowski "universality-conjecture" and noted that it is incompatible with the "mix- 
ing" hypothesis usually invoked to explain the multiplicity dependence of factorial mo- 
ments and slopes. A different explanation of the multiplicity dependence may therefore 
be needed, especially since intermittency and Bose-Einstein effects are now known to be 
closely related. 



4.5 Factorial cumulants 



Normalized factorial cumulant moments, first introduced in M and recently studied in [26 



are defined in ( |2.71| ) as integrals over the background subtracted correlation functions. 



They share with factorial moments the property of "noise suppression" . The normalized 
factorial moments F q can be expanded in terms of normalized cumulant moments K q as 
given in (|2.72|) . This expansion has been found to converge rapidly |26j . The terms in the 



expansion correspond to contributions from genuine q, (q — 1) . . . , 2-particle correlations. 

An analysis of factorial cumulant moments is presented in |26| . Roughly, it is estimated 
that K 2 ~ 0.6,^3 ~ 0.7, #4 < 1.0, K 5 < 1.5 for UAl data at y/s_ = 630 GeV. (The 
inequalities for K4 and K5 are due to the approximation AB by A ■ B in ( f2.72j ) since 
no direct measurements of these averages exist.) Clearly, the two-particle contribution 
to factorial moments is large, but higher orders are not negligible. At the energy of the 
NA22 experiment K 2 is small (~0.2), but K 3 is significantly larger (~0.45). 

From ( 2.72|) it is seen that the contribution F^ to F q from two-particle correlations 



alone can be expressed as 

F 3 {2) = 1 + 3K 2 (4.16) 

F 4 {2) = l + 6fG + 3Kf ; (4.17) 
the contribution Fjp from two- and three-particle correlations to F± as 

F 4 (3) = 1 + 6K 2 + 3KI + AK 3 . (4. 18) 

The difference F q — F^> is a measure for the importance of higher-order correlations. 

Fig. 4.24a shows a cumulant-decomposition of F 3 and _F 4 in UAl-data |[217||. The 



differences between the curves indeed indicate large contributions from genuine higher- 
order correlations. Similar results are observed for NA22 ||142|| in Fig. 4.24b, for p =2 



and 3 and q =3 and 4, in one-, two- and three-dimensional phase space (transformed 
y,y — <p and y — ip — lnpx)- I n general, the difference increases with increasing InM 
(decreasing bin size). This means that the contribution of higher-order correlations to 
the factorial moments increases at higher resolution. An exception are factorial moments 
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in the variable (p, for which only two-particle correlations are found to be non-zero (not 
shown) .11 

The situation is completely different in heavy-ion collisions where, with present accu- 
racy, K q for q > 2 (Fig. 4.24c). The factorial moments are completely dominated by 
two-particle correlations |[217| , |219| , [220|| , implying that higher-order F q contain little or no 
further dynamical information for this type of collisions (see Eq. |2.56|) . 



Using the linked-pair ansatz (see further Sect. 5.1.1.), higher-order cumulant func- 



tions can be expressed as products of K 2 (see also ||221|| for an interpretation in terms of 
independent superposition of sources) 

K q = A q Kt\ (4.19) 

with free constants A q . 

For a negative-binomial (NB) multiplicity distribution, K 2 = 1/fc and the linking 



parameters are fixed numbers given by A^ = (q — 1)! | 109 |. A necessary condition is 



stationarity, i.e. constancy of 1/k. This works well for UAl. For NA22 |142j| , A q is 
observed to increase with decreasing bin size. Approximately constant A q « (q — 1)! are 
found if the data are averaged only over a narrow rapidity region (—0.75 < y < 0.75) and 
the most prominent spike event is excluded. The linked-pair ansatz may thus be a valid 
approximation for high-order correlations in small phase-space domains but not for the 



average over phase space. This would be consistent with the well-documented fact |fL3 
that the negative binomial is often a good parametrization of multiplicity distributions in 
restricted 5y-intervals. 

We shall come back to cumulants and genuine higher-order correlations in Sect. 4.10, 
where they are studied by means of a largely improved methodology. 



4.6 Factorial correlators 



4.6.1 The method 



The moments defined in ( |2.68| - |2T70| ) measure local density fluctuations in phase space. 
Additional information is contained in the correlation between these fluctuations within 
an event. This correlation can be studied by means of the factorial correlators defined in 
( |2.75| ). Correlators are typically calculated at a given 5y for each combination mm' of bins 



with size Sy, and then averaged over all combinations separated by a given bin-distance 



2 Absence of genuine higher correlations has been reported in [218], but at far too low statistics. 



66 



CHAPTER 4. MULTIPLICITY FLUCTUATIONS AND INTERMITTENCY 



D. This is illustrated below. 



H — I — I — I — I — I — I — I — I — I — I — I — I — I — h 



m 



m 



h- 



AY 



In the simple intermittency model (a-model) described in pi, Z2"|, F pq depends on D 
but not on Sy and follows the power law 



F pq oc (AY/ Dy 



The powers <p pq (slopes in a log-log plot) obey the relations |f2T|. [22 



J pq — Yp+q ~ Yp 



(4.20) 



(4.21) 



where the first equality sign is due to the a-model proper, the second to the log-normal 
approximation. According to ( 4.20|) (fin = (f> 2 , so that ([4.21 ) can also be written in the 
form 

(4.22) 



11 



4.6.2 Results 



Preliminary results for pseudorapidity resolution Srj > 0.2 have been reported by the 
HELIOS Collaboration [222 1 . There, however, multiplicities n m had to be estimated from 



the transverse energy -&r,m i n bin m and the average transverse energy (E^) per particle: 
n m = £?T, m / (Et) rounded to the nearest integer. The first direct measurement is from 
NA22 | |223] ]. The In F pq are shown as a function of — In D in Fig. 4.25a-d, for four values of 
Sy > 0.1 (corresponding to M — 10, 20, 30 and 40). Statistical errors (estimated from the 
dispersion of the F pq distribution) are in general smaller than the size of the symbols. F pq 
can be measured up to third order in p and q for 8y=0A binning (Fig. 4.25a). For Sy —0.1, 
the analysis is possible to first and second order only (Fig. 4.25d). The smallest possible 
value for D being equal to the bin size Sy, Fig. 4.25a extends to D = 0.4 and Fig. 4.25d 
to D — 0.1. In all cases, an increase of lnF pq is observed with increasing — InD. Very 
similar results have recently been reported by EMC [140|, EMU-01 [ |155| ] and in | |152| . 

In Fig. 4.26a, the lnF pq are compared at fixed D = 0.4 for four different values of Sy. 
The dashed lines correspond to a horizontal line fit through the data. In agreement with 
the a-model, the F pq indeed do not depend on Sy. Also this result has been confirmed 
on EMC data M and in |T51 



The (^-independence of correlators holds exactly in 
the a-model. Nevertheless, Fig. 4.26b shows that the Sy independence is also valid in 
FRITIOF. For the particular value of D = 0.4, this even happens at very similar values of 

In fact, this property is far from unique to the a-model, but holds 



lnFp g as in the data. 



approximately in any model with short-range order | 224 |. 
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For Fn, the 5y independence is easily derived from a parametrization of the two- 
particle density, integrated over two regions of size 5y separated by D. Using exponential 
short-range order | ]109| ], this gives 

F n -1 oc le- D / L /(e a -l)(l-e- a ) , (4.23) 
a z 



where L is a correlation length and a = Sy/L. According to ( 4.23|) , Fn becomes indepen- 



dent of Sy for a < 1. Since e D I L — > 1 as D — > 0, this form also leads to the deviations 
from ( |4.20| ) observed clS db bending in Fig. 4.25. 

Because of the bending, fitted slopes (p pq have no meaning, except as an indication for 
the increase of F pq in a restricted range. The slopes for two values of 5y are compared 
to FRITIOF predictions in Fig. 4.27a and 4.27b, respectively. As observed earlier for 
the case of univariate moments [142], the FRITIOF slopes are too small also for the 



correlators. This is not surprising since the model does not succeed in reproducing even 
the lowest-order (i.e. two-particle) rapidity correlation function (Chapter 3). 



4.6.3 Interpretation 

Factorial correlators have been analysed in using a suitable parametrization of 
K 2 (2/1,2/2) an d the linked-pair ansatz ||107|| for higher-order correlations. The relations 
Q2.83| - f2~!86| ) of Sect. 2.1.6 then allow to express all correlators in terms of K 2 for arbi- 



trary (p,q). Note that the expressions for F pq contain many lower-order "combinatorial" 
terms which effectively dominate and mask the contribution from genuine (p + g)-order 
correlations. 



A basically similar analysis is presented in p24j], inspired by techniques used in 



quantum-optics. The two analyses have no difficulty to describe basic features of the 
NA22 data, including the sum-rules discussed in Sect. 2.3 which were claimed to be a 
unique test of random cascade models. Fig. 4.28a compares NA22 data on F 2 (5y) and 



Fn(D) with the calculations from [ 224 j . F 2 (5y) is used to fix the parameters of K 2 (5y) 
(assuming stationarity); Fn(D) follows after integration over the appropriate rapidity- 
domains. With the linking-ansatz of ||109|| all other correlators are calculated without fur- 



ther assumptions. An illustrative example is shown in Fig. 4.28b which compares F\ 2 (D) 
from NA22 to the prediction. The agreement is excellent in all cases. This observation is 
confirmed in ||152| . 



According to Q4.21]) , the ratio 4> pq / (j) 2 is expected to grow with increasing orders p and 
q like their product pq. In Figs. 4.27c and 4.27d this is tested for 8y—0A and Sy=0.2, 
respectively. In both cases, the experimental results lie far above the dashed line corre- 
sponding to the expected (f> pq /4> 2 = pq. Since the dependence of ha. Fpq on — hiD is not 
strictly linear, this comparison depends on the range of 5y and D used to determine <p 2 
and <f) pq . In Fig. 4.27d one, therefore, compares a number of fits. Slopes are smaller when 
the upper limit in D is reduced, but do not reach the a-model prediction (dashed line). 
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is 



It can be verified that, at least for the higher orders, the discrepancy with ( [4.21 
mainly due to the second equal sign, derived from a log-normal approximation to the 
density distribution. In a recent paper [fL80||, this approximation has been shown to be 



valid if the density fluctuations are weak or if the density probability distribution is log- 
normal. The NA22 data demonstrate that none of these conditions is fulfilled. 

We conclude that the correlators F pq increase with decreasing correlation length D, but 
do not really follow a power law for D<1. For fixed D, the values of F pq do not depend on 
the resolution 8y, a feature expected from the a-model, but also reproduced by FRITIOF 
and approximately true in any model with short-range order. When the increase of the 
correlators is roughly approximated by a straight line in a restricted interval, the powers 
4> pq increase linearly with the product pq of the orders, but are considerably larger than 
expected from FRITIOF and from the simple a-model. 

The extension of single- variate factorial moments to the multivariate case offers better 
insight into the complicated nature of the correlations. However, the original expectation 
that correlators would help in clarifying the issue of intermittency is not borne out by 
present data. Simple but reasonable models for higher-order correlation functions which 
use the experimental 2-particle correlations as input, have no difficulty in reproducing the 
behaviour of factorial correlators measured e.g. by NA22. 



4.7 Multifractal analysis 

Power-law dependence of normalized factorial moments on the resolution 5 (bin size) is a 
signature of self-similarity in the fluctuation pattern of particle multiplicity. It suggests 
that the probability distribution P(p, 5) of the particle density p has fractal properties. 
For simple Widom- Wilson |[225|| type scaling, P(p, 5) is of the form 

P(p, 5) ~ 8-ep*{plF) , (4.24) 

where (3 and v are critical exponents. All q-th order moments of p (g = 1, 2, . . .) obey power 
laws in S with inter- related exponents depending on q and on (/3, v). This characterizes a 
simple or mono-fractal. Another possibility is a multifractal behaviour, in which P(p, S) 
obeys a relation of the type (cfr. [j226|| ) 

\nP(p,5)/\n5 = /(a), a = lnp/ln<5. (4.25) 



Multifractals, first introduced in ||227|| represent infinite sets of exponents — the multi- 
fractal spectrum — which describe the power- law scaling of all moments of P(p,5). In 
principle, knowledge of the multifractal spectrum is completely equivalent to knowledge 
of the probability distribution. 

Unlike geometrical or statistical systems, multiparticle production processes pose spe- 
cial problems if a multifractal analysis is to be considered. The most obvious one is the 
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finiteness of particle multiplicity in an event at finite energy. Self-similarity, if existent, 
therefore cannot persist indefinitely to finer and finer scales of resolution. 

In multiparticle production P(p, 5) is not directly accessible. At best one can con- 
struct, for a single event of multiplicity n and for given 5, a frequency distribution which 
approaches P(p, 5) only for n — ► oo. For any finite (and usually small) n, the frequency 
distribution and its moments will be subject to statistical fluctuations. 

Since the data sample contains a large number of events, it is obviously recommended 
to consider the event average. This averaging, however, supposes ergodicity. The appli- 
cability of the multifractality concept can, therefore, only be justified a posteriori. 



4.7.1 The method 

A multifractal analysis is based on the the properties of G-moments whose definition is 
given in (|2.73|) . The moments G q (or more often \nG q ) are obtained for each individual 
event at a specified resolution 5y ~ 1 /M and then averaged over the event sample!. 

In the theory of multifractals, the G-moments share with the scaled factorial moments 
the property that self-similar density fluctuations lead, in principle, to scaling behaviour 

G q oc (5y) T " for 5y -> . (4.26) 

In a fractal analysis (see also Sect. 2.4), one therefore determines the slope 

on a double-logarithmic plot, after averaging over all events in the sample. 

A multifractal spectral function is introduced via a Legendre transform defined as 

f(a q ) = qa q -r q , (4.28) 

with 

a.-*j (4.29) 

being the Lipschitz-Holder exponents. The spectral function f(a q ) is a smooth function, 
concave downward, with its maximum at q = 0. It gives a quantitative description of the 
density fluctuations in the dense and in the sparse regions, corresponding to its left and 
right wing, respectively. A wide spectrum reveals a non-smooth density distribution. 

The generalized (Renyi)-dimensions are given by 

D q = -1^ r q = -1^ [qa q - f(a q )} . (4.30) 



3 Note that analyses based on (G q ) or on (In G q ) in general differ and probe different aspects of the 



system under study [ 228 1 
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4.7.2 Experimental results 



G-moments have been studied in a number of experiments ||140| , |152| - |154j , p~56| , |158| , [L59 
g29|-p3l [p^-p38|. As an illustrative example, we show UAl results pi in Fig. 4.29a 



where the event average (\nG q ) is plotted as a function of the resolution (M = 2 M ). 
Starting at a value of for /i = according to definition ( |2.73| ), the moments grow for 
q < 1 and fall for q > 1 as \i increases. The slopes decrease and (InGq) tends to saturate 
for large /i. The saturation is due to an increasing number of bins with content n m = 
or 1, as M becomes large, 



G q (M) 



n 



n 1 q for M — ► oo 



(4.31) 



Fig. 4.29b shows (r q ) and (a q ) for \i = 1 and 2 (small M). The corresponding spectral 
function (f(a q )) is given in Fig. 4.29c as a function of (a q ). The fact that (f(ac q )) does 
not degenerate into a single point implies multifractality in hadron production, at least 
for large bin size 5y = ^jf (small /x). However, for smaller bin sizes, the function turns 
over (i.e., bends upward) and falls into the non-physical region above the dashed line (not 
shown) . 



4.7.3 Universality 

From Fig. 4.29c, it is clear that (/(a)) depends on /i. It also depends on the cms energy 
y^s (or the multiplicity n = 2 V ). In |30| it is conjectured that G- moments (at fixed q) 



show universality in £ = fi — v, however. The latter quantity is directly related to the 
average particle multiplicity per bin, n/M = 2 U ~^ = 2~^. Using a branching model, the 
authors of ]30[ derive the universality relation 

r ff (0 = lnG,(M, v) - lnG 9 (z/, v). (4.32) 

It expresses the scaling behaviour of a function of two variables (/i, v) in terms of a function 
of one variable only. The function r g (£) determines the G-moments as functions of \x for 
all values of v. The validity of ( |4.32|) is claimed to be a strong evidence for self-similarity. 

Fig. 4.30a shows T q as a function of £ for q = ±5. All data points are close to universal 
lines, thus indeed indicating universal behaviour. 

A further prediction is that also (f{ot q )) is universal for fixed £. Fig. 4.30b demon- 



strates that this is not confirmed by the UAl data [233]. In the EMC data [140] the left 
branch shows universality, but not the right one. Besides being more sensitive to univer- 
sality breaking than r 9 (£), the function (f(a q )) also reveals more clearly shortcomings in 
the models. For PYTHIA and GENCL this is illustrated by Fig. 4.30b. 
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4.7.4 Modified G-moments 

As stated earlier, G-moments have the advantage that not only spikes are included in the 
analysis, but also non-empty valleys (for q < 0). Disadvantages are that the moments 
saturate at by — > when the content of non-empty bins approaches unity and that 
statistical fluctuations are not filtered out (see also |239| , |240| ). 



In ||241|| a modified definition of the G-moments is proposed in an attempt to circum- 



vent the problem of statistical noise. "Truncated" G-moments are defined as 



M 



G q =Y, P q m^ m - q) , (4.33) 



m=l 



where is the usual step-function equal to 1 for n m > q and zero otherwise. For very 
large multiplicity n (as in a macroscopic statistical system), n/M 3> q and ( |4.33| ) is in 



practice identical to fl2.73|) . In particle physics, n is a relatively small number and the 



G function exerts a crucial influence on the G-moments. It imposes non-analytical cut- 
offs at positive integer values of q. With the help of a Monte Carlo (ECCO) based on 
the Geometrical Branching Model, the authors show that ln(G g ) now exhibits a linear 
dependence on InM for q > 1, without saturation. 

For q > 1, the linearity of ln(G g ) with InM has been verified on /xp, pp and e + e~ 
data [|230| , |233| . The slopes r q turn out to be very similar in all three reactions and 



roughly equal to r q = — 0.9(q — 1). 



4.7.5 Bernoulli trials and G-moments 

Before we conclude the discussion of experimental characteristics of G-moments, it is of 
interest to inquire in more detail about the dynamical content revealed in multifractal 
analyses. This is best done in a comparison to a model without dynamics. 

Let P(n) be the probability distribution for observing n particles in an initial wide 
interval AY. Let this interval be subdivided into M smaller intervals of size 8y = AY/M, 
each of which contains n m particles (n m = 0, 1, . . . , n; m — 1, 2, . . . , M) with J2m n m — n - 
Assume that for every subdivision of AY the n particles are independently distributed 
over the intervals with probability 1/M. For fixed n, the joint occupation probability in 
M cells is given by (see also 



n' / 1 \ n 

P n ( ni ,..., nM ) = P(n) -^-L— (_j 5^-n). (4.34) 
With ( |4.34| ) the G-moments ( |2.73| ) at fixed n are given by 

n 

G q {n,M) = Mn- q B(n, 1/M; i), (4.35) 



8=1 
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where B(n, 1/M; n) is the binomial distribution. For integer q > 1 one obtains 

G q (n,M) = Mn~ q £S^n yl . (4.36) 



3=1 



<S^ is a Stirling number of the second kind (cfr. ( |2.67| )) and = n(n — 1) . . . (n — j + 1). 
In the Poisson limit of the binomial (n large and 1/M small with n/M fixed), ( 4.36|) 
simplifies further to 

G q {n,M) = Mn~ q {n/M) 3 . (4.37) 



In inclusive analyses the average over P(n) has to be taken in ( [4.36| - fj~3~7l ). This introduces 
the (inverse) moments {n°~ q ) and (n^/n q ) of the multiplicity distribution in AY. 



Numerical studies indicate that ( |4.35| - [4~37l ) reproduce and explain many of the multi 



fractal and universality properties seen in the data. Here we can only give a few examples. 



With respect to the structure of (|4.35| ) it should be noted that the binomial distribution 



is in fact a multifractalo in the sense of ( f4.25|) [ |242| . [226|| . Consequently, "proper" , but quite 



trivial binomial multifractal behaviour will be seen if the data are noise-dominated. This 



is the case in practically all analyses referred to before. This point was recognized in [ [239 
but its full consequences were not further studied. 



From ( |4.37| ) follows immediately that the function T g (/i, v) defined in ( |4.32| ) depends 



only on the ratio n/M = 2~^. Thus, the parameter-free function T q (fi, v) is indeed, but 
trivially, universal in the Poisson limit. It describes accurately the data in Fig 4.30a. 
Being a purely mathematical property of noise, it is not surprising that the usual Monte- 
Carlo models also show this type of universality. T q (fi, v) ceases to be universal in the 
(more general) binomial case, although the deviations remain small in cases of practical 
interest. This probably explains the universality breaking observed in e + e~ Monte-Carlo 
simulations at 1 — 10 TeV in [243]. 



The above considerations can be extended to the modified G-moments defined in 
( f4.33| ). In particular T q (fi, v) remains universal in the Poisson limit. Further numerical 



properties of modified G-moments are illustrated on Fig. 4.31 (a-c). The results shown are 
based on ( |4.37|) further averaged over n with a negative binomial distribution truncated 
at 0. They depend only on (n) in Ay and on the NBD parameter k. 

Fig. 4.31a demonstrates that the modified G-moments can be well approximated by 
power laws. The pseudo-linearity extends over a much larger interval in M than for 
usual G-moments obtained from ( |1.37| ). The improved linearity is due to negative terms 
in the expressions of truncated moments of the Poisson (or binomial) distribution. The 
calculations displayed in this figure coincide (up to an overall normalization factor) nearly 



exactly with the EMC data for M > 8 shown in ||230|| . This proves that the "clear 



4 This is easily verified using Stirling's approximation to n! which is notoriously accurate even for quite 
small n. 
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asymptotic power-law behaviour of (G q ) characteristic for a self-similar system" |[230f | is 
in fact due to Bernoulli noise. 

Fig. 4.31b shows the "Ochs-Wosiek" plot for modified G-moments. Here again, the 
quasi-perfect linear relation between lnG q and lnG2 is seen to be a characteristic of 
Bernoulli trials. This linearity property holds in fact for any combination hiG q and 
lnG q i. The exponents r q derived from power-law fits to the "data points" in Fig.4.31a 
are shown in Fig. 4.31c. The line is a fit with the form r q = —C(q — 1). The slope 
C is a slowly- changing function of the NBD-parameter k with a value around 0.9, as 
experimentally observed in Subsect. 4.7.4 above! 



4.7.6 Evaluation of noise and connection between F q and G q 

The self-similar property of multiparticle production at high energy can, in principle, be 
investigated by F-moments and by G-moments. The power-law behaviour of the scaled 
F-moments provides evidence for a self-similar cascading process of dynamical origin. The 
G-moments, as an ingredient of fractal theory, are designed to describe the multifractality 
aspect of high multiplicities. In the real environment of high energy collisions, however, the 
multiplicities are rather low and the G-moments are dominated by statistical fluctuations. 

The F- moments are defined for integer powers q > 1, the G- moments for all real powers 
q. In order to establish a connection to the F-moments, the powers q are restricted to 
integer values of q > 1 here also for the G-moments. 

The number n m of particles per subdivision 5y = Ay/M has to be equal to, or larger 
than q (n m — q + k, k — 0, 1, . . .) for the F- and G-moments. Functions 

B (M) = { Q q+k (M,n) 
n q 

are defined from the number of bins Q Hm {M, n) containing n m = q + k particles in an event 
of multiplicity n in the total phase-space region Ay, normalized by n q and averaged over 
all events. They express the basic fractal structure of the data, if they show a power-law 
behaviour of the form 

B qjk (M) oc M X ^ M) . (4.39) 

In order to suppress statistical fluctuations, the G-moments can be defined as the event 
average over (4.33), or, equivalently, as 

oo 

(G q (M)) = Y,B q , k (M)(q + k) q . (4.40) 

k=0 

They are proportional to M~ Tq for large M. The F-moments are defined as 
(FJM)) = M- 1 Y / (^(n m -l) (n m -g + l) 

\ <a n jla {n/M)q i y ) 
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or, equivalent ly, as 

(F q (M))^M^J2B q , k (M)^ft . (4.42) 

fc=0 

They are proportional to M^ q for large M (note, however, that (4.41) is different from 
the form (2.68) of the F-moments generally used). 

When (4.33) and (4.40) to (4.42) are applied to the data they should show a power-law 
behaviour for large M, if there are fractal structures present in the data. 

The dynamical contribution to the G-moments can be expressed by 

{Gq)dyn = (i^ M(1 ^ ' (443) 

where (G q ) st can be determined by distributing the n particles of an event randomly in 
Ay. The randomization procedure destroys short-range particle correlations, but does 
not alter the Bernoulli nature of the particle repartition in smaller bins discussed in 
Subsect. 4.7.5. As a result, this method does not eliminate the binomial, noise-induced 
multifractal behaviour, but just gives its behaviour. 

When (G q ) st is equal to (G q ), a trivial "dynamical" effect remains: a flat dn/dy leads 
to a probability 1/M for a particle to be in a given bin and {G q ) dyn = M 1 ^ q . 

The dynamical contributions to the slope r q can be expressed by 

T^ a = r q -rf + q~l , (4.44) 

where r^* is the statistical part of the slope. Subtracting the statistical contribution from 
r q gives 

T Q ~ '■ q • q 

which can be directly compared to the slopes <p q obtained from the F-moments 

Fig. 4.32a gives a comparison [[241|| of <p q (crosses) and (q — 1 — T^ yn ) (full circles) from 
ECCO simulation results and shows that the deviation of r^ yn from q — 1 is indeed close 



r q - 7? = r„ dyn - q + 1 « -(f) g (4.45) 



to the deviation of (j) q from zero. This observation gains support from the UAl [233 



hA [ |152| ] (Fig. 4.32b) and AA |156| ] analysis. The remaining difference can be attributed 



to the difference in the definition of (F q ) and (G q ). 

Fig. 4.32a, however, also shows that r^ yn cannot be simply replaced by r q (open circles) 
in a quantitative analysis. This is in agreement with the observation of Subsect. 4.7.5, 
but has not been taken into proper consideration in recent experimental application on 
lh p30| , hA pi and AA [|3§ collisions. 



To summarize the present experimental findings, the data indicate that the multifrac- 
tal spectral function f(a) has, at least for large bin sizes, the properties expected from the 
theory of multifractals. The function f(a) is very sensitive to violations of universality 
and to details of present Monte-Carlo models. However, with the methods used so far, the 
multifractality analysis breaks down at finer resolution. The finite multiplicity effect— 
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statistical noise — overwhelms and it is difficult to disentangle it from dynamical features. 
An advantage is that the (G q ) can probe holes in the distribution, not just spikes. A 
recent extension of the definition of the G-moments filtering out high multiplicities can 
claim some success in extracting the dynamical component. The advantage is that (G q ) dyn 
lends itself more readily to (multi)fractal interpretation and direct extraction of the Renyi 
dimensions according to (2.114), while (F q ) is more closely related to the correlation func- 
tion. Fig. 4.32 can serve as a rough link between the two. Whereas a higher-dimensional 
analysis could be useful also here, the low average multiplicity even in the highest energy 
experiments presently precludes further progress in this direction. 



4.7.7 Universal multifractals 



At first sight interesting approaches, recently applied |5T| to obtain the degree of multi- 
fractality (or Levy-index) fi in multiparticle production, are the methods of Probability 
Distribution Multiple Scaling (PDMS) and Double Trace Moments (DTM) 0. In the 
first method, the fundamental scaling law is written in terms of the probability for the 
number of particles n m in bin m at resolution M to be larger than a certain threshold 
n th = M\ 

P(n m (M) > M 7 ) oc M" c(7) . (4.46) 

The statistical function 0(7) is the codimension function describing the sparseness of large 
intensities n m . Like the factorial moments ( |2.68| - pTTQ] ) or the extended G-moments ( |4 . 3 3| ) , 
the DPMS method is a straightforward filter for spikes of large n m . 

The PDMS method is closely related to Large Deviation Theory, a topic in probability 
theory and of much theoretical interest in statistical mechanics |[244j| . Equation (|4.46| ) 



expresses a Level-2 Large Deviation property, describing deviations of the "empirical 
measure" n m from the infinite sample probability density; 0(7) is related to a generalized 
entropy. 



Fig. 4.33a shows data |51[] at different charge multiplicity {n = 6 and 14 are given as 
examples) presented in a double- logarithmic plot, for various threshold values n t h = M 1 . 
In spite of limitations on statistics and on multiplicity n m , a region of linearity can be 
seen for all multiplicities and thresholds. This is claimed to be evidence for PDMS. 

When 0(7) is smaller than the topological dimension D of the embedding space, it 
is possible to define a function D{p/) = D — 0(7) corresponding to the classical fractal 
dimension. If a sample of N s events is used in the analysis (instead of one event), 0(7) 
can become larger than the topological dimension since different events can contribute to 
the same bin m. In Fig. 4.33b the function 0(7) is shown for the same multiplies. The 
dotted line corresponds to c(%) = D + D s , where D s is the sample dimension defined as 
N s = M Db . For n = 20, e.g., this limit is crossed for a threshold of «-th=3 with N s = 15. 
Singularities of that type are called "wild" singularities (not arising from Poisson-like 
fluctuations). 

A parametrization of 0(7) in terms of two parameters is provided by the theory of 
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universal mult ifr act als ^ 



c(t) = Ci( r -^r + -j /WM^I 

V//'Ci //y 

= dexp^-l) /or/i = l, (4.47) 

with - + ^ = 1 and < u < 2; u is the Levy-index giving the degree of multifractality 
and C\ is the codimension of the average field. The two parameters can be obtained from 
fits to the results given in Fig. 4.33a, but turn out to be highly correlated. 

The Levy-index fj, can, however, be determined independently of C\ with the help of 
the Double Trace Moments ||50|| , a generalization of the G- moments. They are defined 



in ( |2.126| ). A DTM analysis of a/s = 16.7 GeV data yields [i- values ranging from 0.4 to 



about 0.9, increasing with multiplicity n. Such values are far from monofractality (// = 0), 
but considerably below /n ~ 1.6 obtained in Subsect. 4.3.3 from factorial moments. 

For multifractal theory, it is important to know whether the limit fi = 1 is crossed 
(signalling "hard unbounded" singularities) or asymptotically approached from below (in- 
dicating "soft bounded" singularities). This question cannot be answered at low energies 
and needs high-energy, high-multiplicity data. 

Extension to higher-dimensional space, though difficult in practice, is necessary since 
singularities can easily be washed out if a particular variable is not sensitive to them. 

In spite of the interesting potential of the "Universal Multifractal" idea, one should 
keep in mind that the method suffers from the same limitations as the multifractal method 
based on G-moments. This is easily illustrated by considering again the Bernoulli-trials 
model discussed in Subsect. 4.7.5, above. For the pure binomial noise (4.34), one has 

P B (n m (M) > nth) = h/ M (n t h, n-n th + l), (4.48) 



where I x (a : b) is the incomplete beta function fT6 



The logarithm of ( 4.48|) is approximately linear in InM for reasonably large M. Equation 



( |4.48[ ) not only has all the features of the data plotted in Fig. 4.33a, but even agrees 
numerically quite well. The co-dimension function 0(7) is, for this simple model, approxi- 
mately equal to n t ^, implying constant differences between the slopes for successive values 
of n t h. The data in Fig. 4.33b show exactly this property. 

Double Trace Moments are easily calculated in the Bernoulli model. We find that the 
"Levy-index" /i is a smoothly increasing function of the event multiplicity n crossing the 
"hard unbounded" value /i = 1 near 77 = 30. We conclude that the data in Fig. 4.33 are 
merely reflecting statistical noise. 

Dynamically useful information could possibly be extracted if Double Trace factorial 
moments were used instead of the usual moments. This is easily verified on the simple 
Bernoulli model and, of course, applies to G-moments as well. 
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4.8 Density and correlation strip-integrals 
4.8.1 The method 

A fruitful recent development in the study of density fluctuations is the density and 
correlation strip- integral method [ |130| , [245| , |246|| . By means of integrals of the inclusive 
density over a strip domain, rather than a sum of box domains, one not only avoids 
unwanted side-effects, such as splitting up of density spikes, but also drastically increases 
the integration volume (and therefore the accuracy) at a given resolution. 

Consider first the (vertical) factorial moments F q defined, for an analysis in one di- 



1 m (r) ig]\ -I m I n i dy i pq(y 1 ...y g ) 
Wv) = -b E iH = T7 E > . ; ; — • (4-49) 



M m=1 (n m )« M m=1 f Ut dy tPl (yi) ■ ■ ■ PM 



The integration domain VL-q = X)m=i thus consists of M (/-dimensional boxes VL m of 
edge length 8y. For the case q = 2, Qb is the domain in Fig. 4.34a. A point in the m-th 
box corresponds to a pair (yi, y 2 ) of distance \yi — 2/2 1 < 5y and both particles in the same 
bin m. Points with 1 2/1 — 2/2 1 < 5y which happen not to lie in the same but in adjacent 
bins (e.g. the asterix in Fig. 4.34a) are left out. The statistics can be approximately 
doubled by a change of the integration volume to the strip-domain of Fig. 4.34b. For 
q > 2, the increase of integration volume (and reduction of squared statistical error) is in 
fact roughly proportional to the order of the correlation. The gain is even larger when 
working in two or three phase-space variables. 

In terms of the strips (or hyper-tubes for q > 2), we define as (vertical) density integrals 

/ ILidyip q (yi,...,y q ) 

F*{6y) = -f* (4.50) 

/ TlidyiPiiyJ ■ ■ ■ PiiVg) 

and, similarly, the correlation integrals K^(5y) by replacing the density p q (y%, . . . ,y q ) by 
the correlation function C q (yi, . . . ,y q ). (Note that in the literature the term "correlation 
integral" is often also used for the F q (5y).) 

These integrals can be evaluated directly from the data, after selection of a proper 
distance measure (|t/j — yj\, [(yi — yj) 2 + (0j — 0,,) 2 ] 1 / 2 , or better the four-momentum 
difference Q 2 j = —{pi — Vj) 2 ) an d after definition of a proper multiparticle topology, 
the snake integral ||107|| , the GHP integral [ 130|| , or the star integral ||247|| as shown in 



Figs. 4.34c-e, respectively. 
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4.8.2 Results 



As an example, F±(8y) is compared to Ff(8y) (and Ff , F 3 S )) for the NA22 spike event ||119 
in Fig. 4.35a (no error-bars are shown, because it is one event). Depending on whether the 
prominent spike lies entirely in one bin or is split across two, F 4 shows large fluctuations. 
These are practically absent in Ff (Fig. 4.35b). Large improvement in one-dimensional 



(rj) and two-dimensional (rj — ip) analysis is also observed in | J248|| . 

How much the statistical errors are reduced can be seen on Fig. 4.36a where the NA22 
data |p.42|| are plotted as a function of — InQ 2 , with all two-particle combinations in an 



ra-tuple having Q 2 , < Q 2 ||130|| . The following observations can be made: 



i) the errors and fluctuations are indeed largely reduced, as compared e.g. to Fig. 4.20a. 

ii) with the (one-dimensional) distance measure Q 2 , the moments show a similarly 
steep rise as in the three-dimensional analysis (e.g. Fig. 4.11c). 

iii) Contrary to the results in rapidity, positives and negatives behave very similarly 
here (only negatives are shown in Fig. 4.36a), but are now much steeper than all-charged. 

iv) F 2 S is flatter for (H — ) than for all-charged or like-charged combinations. 

The first two observations demonstrate the strength of the new method and the advan- 
tage of using the proper variable. The second two observations directly demonstrate the 
large influence of identical particle correlations on the factorial moments. These results 
agree very well with results from the UAl collaboration [ I4H ] shown in Fig. 4.36b and 
with lh results j77l|249fl. 



Monte-Carlo simulations with FRITIOF 2 show the following (see Fig. 4.37 for the case 
ofF|). The default "plain" version is unable to describe the all-charged NA22 data, but a 
"biased" version (including misidentified Dalitz decay + 0.25% undetected 7-conversions) 
comes closer to the data. However, not unexpectedly, both versions fail completely in 
describing the like-sign data, where the model stays way too low. On the other hand, F% 
for the (H — ) combination is largely overestimated when 7-conversions are included, but 
saturates without. 



4.8.3 Transverse-momentum and multiplicity dependence 



As in Subsect. 4.4.2 (Fig. 4.20c), UAl |210|| has studied the 02 dependence on the average 
transverse momentum pi of the event, but now in terms of density integrals in Q 2 . In 
contrast to the strong decrease (and subsequent slight increase of ^2) with increasing pp 
observed for the one-dimensional analysis in Fig. 4.20c, a strikingly flat behaviour (and 
slight increase above 0.6 GeV/c) is observed for the data (full circles) in Fig. 4.38a. The 
discrepancy of PYTHIA (open circles) is even stronger here than in Fig. 4.20c. The slope 
4>2 starts at even negative values for small pi, but increases fast with increasing pi to 
reach values overestimating (fi 2 at pi ~ 0.5 GeV/c. 

A similar disagreement is observed for the multiplicity dependence in Fig. 4.38b. While 
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the UAl data (full circles) decrease with increasing n, PYTHIA predicts a strong increase. 
In Figs. 4.38 c) and d), it is shown that this violent discrepancy between PYTHIA and data 
is mainly due to like-sign pairs, so to the way Bose-Einstein correlations are incorporated 
into the model. 



4.8.4 Bose-Einstein correlations versus QCD effects 

Of particular interest is a comparison of hadron-hadron to e + e~ results in terms of same 
and opposite charges. This is shown in Fig. 4.40 for q = 2 UAl and DELPHI data 
in ||25U|] (note that in this figure the derivative of (|4.50|) is presented in small Q 2 bins). 



An important difference between UAl and DELPHI can be observed on both sub-figures: 
For "large" Q 2 {> 0.03 GeV 2 ), where Bose-Einstein effects do not play a role, the e + e~ 
data increase much faster with increasing 2 log(l/Q 2 ) than the hadron-hadron results. 
For e + e~, the increase in this Q 2 region is very similar for same and for opposite sign 
charges. At small Q 2 , however, the e + e~ results approach the hadron-hadron results. 
The authors conclude that for e + e~ at least two processes are responsible for the power- 
law behaviour: Bose-Einstein correlations at small Q 2 following the evolution of jets at 
large Q 2 . In hadron-hadron collisions at present collider energies only Bose-Einstein effects 
seem relevant. 

Since string fragmentation causes an anti-correlation between same-charged particles, 
it is of interest to compare e + e~ results to JETSET in terms of strip integrals for the 



different charge combinations, separately. This has been done in ||250|| and, indeed, the 
Monte-Carlo results level off at small Q 2 and fall below the data for the same-charge 
results, while they describe the opposite-charge data perfectly well (not shown here). 

The exact functional form of _F 2 S is derived from the data of UAl | |143| | and NA22 [ |142| ] . 
again in its differential form!, in Fig. 4. 40. Clearly, the data favour a power law in Q over 
an exponential, double-exponential or Gaussian law. 



If the observed effect is real, it supports a view recently developed in ||165f| . There, 
intermittency is explained from Bose-Einstein correlations between (like-sign) pions. As 
such, Bose-Einstein correlations from a static source are not power behaved. A power 
law is obtained i) if the size of the interaction region is allowed to fluctuate, and/or ii) if 
the interaction region itself is assumed to be a self-similar object extending over a large 
volume. Condition ii) would be realized if parton avalanches were to arrange themselves 
into self-organized critical states ||251|| . Though quite speculative at this moment, it is 



an interesting new idea with possibly far-reaching implications. We should mention also 
that in such a scheme intermittency is viewed as a final-state interaction effect and is, 
therefore, not troubled by hadronization effects. 

The effect on the factorial moments of adding Bose-Einstein correlations in FRITIOF 
is convincingly demonstrated for heavy- ion collisions in ||252|| . Because of the large number 



5 In fact in this differential form F^{Q 2 ) is identical to R(Q 2 ) usually used in Bose-Einstein analysis. 
The only difference is that it is plotted on a double-logarithmic plot, here. 
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of collision processes, other correlation effects are expected to play a much reduced role for 
this type of interaction and Bose-Einstein correlations, as a collective effect, can become 
the dominant source of non-statistical fluctuations. Also from these results it is clear that 
more than one fixed interaction-volume radius is needed to reproduce the experimental 
results. 

In perturbative QCD, on the other hand, the intermittency indices 4> q , are directly 



related to the anomalous multiplicity dimension 7 = (Qol s / , k) 1 I 2 [|253| - |257|| and, therefore, 
to the running coupling constant a s . In the same theoretical context, it has been argued 
254| - |257|| that the opening angle x between particles is a suitable and sensitive variable 



to analyse and well suited for these first analytical QCD calculations of higher-order 
correlations. It is, of course, closely related to Q 2 . 

A first analytical QCD calculation [ 254| , |255| is based on the so-called double-log- 
approximation with angular ordening 
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and on local parton-hadron-duality [|100j1 . A 
preliminary comparison with DELPHI data |[258|| gives encouraging results, even including 
an estimate for the running of the strong coupling constant a s . 



4.9 Correlations in invariant mass 

The previous section has illustrated the advantages of the correlation integral method with 
a "distance" measure directly related to the invariant mass of the particle system. The 
results give additional support to the Fialkowski conjecture, mentioned in Subsect. 4.3.2, 
from which could be anticipated that dynamical effects are most clearly revealed if the 
correlation functions and factorial moments are directly analysed in terms of Lorentz- 
invariant variables. 

Evidently, there are many arguments in favour of invariant mass as a dynamical vari- 
able rather than the single-particle variables often used in early intermittency studies. 
Resonances, the cause of most of the correlations among hadrons, and threshold effects 
appear at fixed values of mass; Bose-Einstein interference correlations depend on four- 
momentum differences; multiperipheral-type ladder diagrams are functions of two-particle 
invariant masses, and so on. 

The idea to study correlations as a function of invariant mass was, to our knowledge, 
first proposed in | p59| , f260[| . The authors introduce a method which is technically a dif- 
ferential version of the correlation integral method. It focusses directly on the correlation 
functions (cumulants) rather than on the inclusive density as in ( |4.50| ). Starting from the 
definition fl2.21|) , one defines the correlation function 

^ 2 (M inv ) = p 2 (M inv ) - Pl ® pi(M inv ), (4.51) 

obtained after integration (in a suitable region of phase space) of C2(pi,p2) over all vari- 
ables except M mv . Here, p2(M mv ) is the familiar normalized 2-particle invariant- mass 
spectrum. The "background term" piCg>pi(M mv ) is the integral of Pi(pi)pi(p2) with M mv 
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fixed. For the data shown below, it is obtained from "uncorrelated" ("mixed") events, 
built by random selection from a track pool. The same method is used in evaluating 
the denominator in ( [4.501 ). Higher-order correlations are obtained in a completely analo- 
gous manner. We further utilize the function K 2 (M- IWJ ) = ^M-Mjnv) / Pi ® Pi(^inv)' the 
normalized factorial cumulant of order two. 



The analysis in ||259|| , based on low statistics pp data at 205 GeV/c, demonstrates that 
K 2 (M[nv) an d ^2 O^inv) follow an approximate power law, written by the authors as 

K 2 (M inv ) = (M? ny r^-\ (4.52) 

The notation reminds of the interpretation of ( |4.52| ) in terms of the Mueller- Regge formal- 
ism (for details see |[259|| ). The power «x(0) is the appropriate Regge- intercept, X = R 
for non-exotic pairs and X = E for exotic ones. The ratio K 2 / K 2 was further seen to 
fall as M^ v , consistent with «r(0) — «e(0) = 1. Not relying on Mueller- Regge theory, the 
authors argued that most of the correlations at small M mv are due to resonance decays 
into three or more pions and to interference of amplitudes [|260|| . 



The results already obtained in |[259| 1 clarify several issues which have troubled the 
interpretation of intermittency data. Among others, they demonstrate that different 
charge-states should be treated separately since the M mv dependence is very different. 
This fact, obvious in M mv but much less so in rapidity, was not fully appreciated in 
early intermittency analysis and the crucial importance of like-sign particle correlations 
remained hidden in "all-charged" analyses. 

The method of |259| has now been applied by NA22 and DELPHI M . Fig- 
ure 4.41 shows data on K 2 (M[ nv ) for a combined sample of non-diffractive n + /K + p 
collisions at 250 GeV/c in the central cm. rapidity region — 2 < y < 2. K 2 ~ (M mv ) 
has a prominent p° peak, but is quite flat near threshold. The peak in the first bin 
of sub-figure (a) is attributed to contamination from Dalitz-decays and 7-conversions. 
K 2 -(M inv ) falls much faster. A fit of K 2 ~ (M? nv )~P yields /3" = 1.29 ± 0.04, 
/?++ = 1.46 ± 0.03, /?+- = 0.17 ± 0.02, in agreement with |259| and consistent with 
the relation «r(0) — o:e(0) = 1. 

NA22 also finds that cuts on transverse momentum or relative azimuthal angle Sip 
strongly affect the shape of i^^~(M mv ), but have little effect on K 2 (M mv ) for M mv < 
0.5 GeV/c 2 . This means that K 2 (M mv ) at small M mv is essentially a function of M mv 
(or Q 2 ) only, illustrating once more the advantage of M mv compared to other variables. 

The data in Fig. 4.41 confirm the conclusion of Sect. 4.8 that the correlations in like- 
charge systems are at the origin of the strong increase of factorial moments for small 
invariant masses. Whether Bose-Einstein effects are solely responsible for the differences 
between (tt ± tt ± ) and {jt + it~) pairs is not so evident. It suffices to consider ||260|| the 
contributions from decays of various resonances to realize that the M mv - dependence near 
threshold for "exotic" particle systems must be stronger than for "non-exotic" ones. In a 
dual Regge picture, such differences translate into very different values of the respective 
Regge intercepts as in ( |4.52|) . It remains, therefore, to be verified if the M mv -dependence 
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of the data can be explained as a superposition of a "standard" Regge-type power law 
and a conventional Bose-Einstein enhancement. 



As pointed out in |[263|| , there is a feasible way to test this and even to give access 
to the relative strength of BE interference and exotic like-charge irir interaction. The 
idea is that particle combinations exist which are either a) exotic, but not identical (e.g., 



K + 7r + or K 7i pairs) or b) identical, but not exotic (1 = 7r°7r° pairs). NA22 ||264|| and 



ALEPH ||265|1 data indicate that very-short-range correlations are indeed absent in the 



exotic Ktt channel. This supports Bose-Einstein correlations rather than exotic Regge 
behaviour, but the point deserves further investigation. 

The possibility that the correlation functions depend mainly on invariant mass has 
interesting further consequences. These were analysed in [|185|1 . Taking in ( 14.511 ) 

C 2 (M [nv ) oc (M?^- 1 BE(M inv ), (4.53) 

with BE a conventional Bose-Einstein factor, exponential in Q, good agreement is obtained 

with the NA22 second-order correlation integral data of ( )-pairs. Integrating the 

correlation function over all variables except 5y gives F 2 (5y) which also fits the data. 
Although C 2 does not explicitly depend on the transverse momentum of the particles, it 
turns out that F 2 (Sy,pTi,pT 2 ) is larger and more steeply increasing than F 2 (Sy) for small 
Sy and small px's. The opposite happens for large px's. This is the "low-px intermittency 
effect" seen in the NA22 and UAl data (cfr. Fig 4.20a and Subsect. 4.4.2). The explanation 
is simple: under the stated hypothesis, small pi for the two particles in a pair means, 
on the average, smaller invariant mass than for unrestricted transverse momentum and, 
therefore, larger and shorter-ranged correlations in rapidity. Enhanced intermittency 
follows as a consequence of kinematical cuts! The influence of the Bose-Einstein factor 
is easily checked in this simple model. It is found to be necessary in order to reproduce 
the correlation integral data and F 2 for restricted pi but has, as expected a priori, very 
little influence on the p^- integrated F 2 (5y). This explains early controversy over the role 
of Bose-Einstein effects in one-dimensional factorial moment analyses (Subsect. 4.4.1). 

A study of the invariant-mass dependence of the two-particle correlation function has 
for the first time given clear indications as to why the hadron-hadron Monte-Carlo models 
fare so badly when confronted with factorial moment data. For example, Fig 4.42 shows 
NA22 data for K 2 ~(M [nv ) and K$~(M- mv ) compared to FRITIOF. The predicted shape 
of K 2 ~(M[ nv ) is very different from the data, especially in the p° region. It shows an 
enhancement at low mass which causes the correlation function to drop much faster than 
seen in the experiment. In the model, this structure is traced back to reflections from rj, 
i]' and u resonances. The model also fails to describe K 2 (^q nv ) since correlations are 
very weak or even negative, except for a threshold enhancement due to //-decays. 

These examples suffice to demonstrate that FRITIOF (or rather JETSET) has serious 
shortcomings and is unable to reproduce two-particle correlations in invariant mass. For 
correlations in rapidity and azimuthal angle this was seen earlier (Sect. 3.1), but the rea- 
sons remained obscure, mainly because of the insensitivity of these variables to dynamical 
correlations at small mass. 



4.9. CORRELATIONS IN INVARIANT MASS 
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A study of the correlation function in terms of invariant mass clarifies the situation 
considerably. For NA22, the model was known to overestimate significantly the production 
rates of p° and 77 mesons | |266| | and presumably also those of rj' and uj for which no direct 



measurements exist (see also [|171|| for hA collisions). This is now seen to distort heavily 



the M^ nv -dependence of K 2 . Also Bose-Einstein low- mass enhancements, most likely 
responsible for the fast drop of K 2 (Mj nv ) in the threshold region, are not included 
in the FRITIOF model commonly used. Finally, we note that the values of K 2 in the 
considered mass interval are much smaller than the data. This is related to the width of 
the charged particle multiplicity distribution which is known to be too small in FRITIOF. 
It affects the global magnitude of factorial moments and cumulants. 

In Fig. 4.43a, the discrepancy is shown to be quite similar for (H — ) correlations in 



e + e~ collision ||262|| and JETSET. As in hh collisions, the correlation is underestimated 
in the mass region below the p°. This discrepancy can be cured by decreasing the rj' and 
p° production and increasing uo production in JETSET. 

Fig. 4.43b gives the like-sign correlation for the data and JETSET without BE cor- 
relation. For M inv < 0.6 GeV/c 2 . The experimental data are considerably higher than 
JETSET. This can be attributed to Bose-Einstein interference. However, it is striking 
that JETSET also predicts a strong rise towards threshold even without Bose-Einstein 
correlations. This is the tail of the QCD effect also seen in Fig. 4.39 and mainly due 
to multijet events. The difference between JETSET and data can indeed be removed by 
including BE correlations in the model, but the Q 2 cut used (Q 2 > 0.04 2 ) is too high to be 
able to distinguish a power law from an exponential or Gaussian, as is done in Fig. 4.40. 

To summarize, the above proves that the failures of models such as FRITIOF and JET- 
SET with respect to factorial moment and correlator data (Subsects. 4.2.3 and 4.6.3), are 
not necessarily due to "novel" dynamics. They are in first instance a consequence of a 
variety of defects — such as incorrect resonance production rates and absence of identi- 
cal particle symmetrization — which belong to "standard" hadronization phenomenology. 
These defects should be eliminated before "new physics" can be claimed. 

For e + e~ annihilation at LEP energies, we have found that models such as JETSET- 
PS are much more successful than for all other processes. Besides the evident fact that 
this process is much better understood theoretically, QCD effects dominate and the model 
parameters are much better tuned to the data. Still, serious, recently observed discrep- 
ancies of e.g. JETSET-PS with LEP measurements on particle and resonance production 
rates are a clear sign that hadronization in e + e~ is in fact less well understood than 
commonly stated and needs improvement. It could be rewarding to investigate carefully 
and differentially the invariant mass dependence of the correlations using the sensitive 
methods now available. 

Originating mainly from the low invariant mass region (typically < 1.5 GeV/c 2 ), it 
is not impossible that the observed correlations are quite independent of the process 
initiating the primary colour separation in the collision, being dominated by strong final- 
state interactions. This would explain "universality" in the sense discussed earlier. 
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Many authors argue that "intermittency" is somehow connected to (nearly scale- 
invariant) perturbative QCD-cascading. Others strongly contest this view on the argu- 
ment that QCD cascades have a limited extent even at LEP energies and are dominated by 
a very small number of "hard" emissions. In the former case, one may expect significant 
differences in the correlation functions at low mass for e + e~, on the one hand, and for hh, 
hA and AA collisions, on the other. Preliminary e + e~ data, mentioned in Sect. 4.8, seem 
to support the last opinion. Whatever the final outcome, if differences are found, they 
should be used to clarify the respective roles of perturbative and hadronization phases in 
the different types of collision processes. 



4.10 Genuine higher-order correlations 

Multiparticle production in high-energy collisions is one of the rare fields of physics where 
higher-order correlations are directly accessible in their full multi-dimensional character- 
istics, under well controlled experimental conditions. 

Three-particle correlations have been observed in the form of short-range rapidity 
correlations and higher-order Bose-Einstein correlations, but evidence for genuine higher- 
order correlations (i.e., after subtraction of all lower-order contributions) is very limited. 
While it was found to be completely absent in heavy-ion collisions, first evidence was 
given in Sect. 4.5 for their existence in hh collisions. 

The correlation integral method turns out particularly useful for the unambiguous 
establishment of genuine higher order correlations in terms of the normalized cumulants 



K q (Q 1 ), when using the star integration [ |247|| 



/ UidySn . . . Qi q C q (yi, . . . , y 2 ) 

K* g (Q 2 ) = ^ , (4.54) 

/ TlidyiQn . . . Q lq pi (yi) . . . p\ (y q ) 

with Qij = Q(Q 2 — Qlj) restricting all q — 1 distances Q\j to lie within a distance Q 2 of 
the position of particle 1. The star- integral method combines the advantage of optimal 
use of available statistics and minimal use of computer time. Since higher acccuracy is 
abtainable, dynamical structures in the correlations can be studied in greater detail than 
with conventional methods. 



Non-zero values of K*(Q 2 ) increasing according to a power law with decreasing Q 



2 



are indeed observed for E665 for third order [ 249 ] and for NA22 up fifth order ||267|| (see 
Fig. 4.44 for the latter). 



4.11 Summary and conclusions 
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1. Intermittency, defined as an increase of normalized factorial moments with increas- 
ing resolution in phase space, is seen in all types of collision. Intermittency is 
a 3D phenomenon. The anomalous dimensions are small (d q —0. 01-0.1) in a one- 
dimensional analysis, but the factorial moments are considerably larger, and their 
resolution-dependence more power-like, in two- or three-dimensional phase space. 
Self-similarity in the dynamics of multiparticle production is an attractive but not 
fully proven explanation. 

2. The factorial- moment method is very sensitive to biases in the data. These have to 
be studied in detail before final conclusions can be drawn. Because of its sensitivity, 
the method has in fact proven to be very helpful in detecting and tracing such biases. 

3. The logarithms of factorial moments satisfy a possibly dimension-independent lin- 
ear relation which allows to determine directly ratios of anomalous dimensions. The 
observed order-dependence of anomalous dimensions excludes a second-order phase 
transition (as treated in | 183|| ) as the origin of intermittency. Valid in random cas- 



cade models, the Ochs-Wosiek relation shows that correlation functions of different 
order are inter-related in a specific hierarchical structure. An explanation in terms 
of dynamics has not yet been found. 

In hadron-hadron collisions, factorial moments and intermittency indices depend 
strongly on transverse momentum and are largest for low transverse momentum 
hadrons. This effect, also seen in rapidity- and azimuthal correlations, needs further 
study in e + e~ collisions and in parton shower Monte Carlo's. There are serious 
indications that "low-pp intermittency" is a reflection of the very strong dependence 
of correlation functions for identical particles on invariant mass. This has to be 
examined in more detail. 

The multiplicity dependence in hh collisions agrees with what is expected from 
mixing of independent sources. However, the Fialkowski observation on possible 
universality casts some doubt on this type of interpretation. For a given density, 
heavy-ion collisions show more intermittency than hadron-hadron collisions, possibly 
as a result of Bose-Einstein interference or other collective effects. 

Factorial cumulants are direct measures of genuine higher-order correlations. These 
are present in hadron-hadron collisions, in particular for small phase-space domains, 
but seem to be absent in heavy-ion collisions. 

Factorial correlators reveal bin-bin correlations. The correlators F pq increase with 
decreasing correlation length D, but only approximately follow a power law for 
D < 1. For fixed D, the values of F pq are independent of resolution 8y, a property 
predicted in the a-model, but also shared by models with short-range order such as 
FRITIOF. The powers <ft pq increase linearly with the product pq of the orders, but 
are considerably larger than expected from FRITIOF and from the simple a-model. 
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8. A recent extension of the definition of G-moments filtering out high multiplicities, 
claims success in extracting the dynamical component. However, under the condi- 
tions prevailing in present hadroproduction experiments, multifractal and general- 
ized multifractal methods seem unable to overcome the overwhelming dominance of 
statistical fluctuations. 

9. The correlation (or density) strip integral strongly reduces statistical errors, as well 
as fluctuations due to splitting of spikes. Using the squared four-momentum Q\- as a 
distance measure, an increase similar to that found in three-dimensional analyses is 
observed. This increase is caused by correlations among like-charged particles. Bose- 
Einstein interference must contribute significantly to the intermittency effect but is 
not power-behaved in the conventional approach. Power-law behaviour in Bose- 
Einstein interferometry would imply a random superposition of "emission centres" 
with possibly fractal properties. Parton avalanches in a self-organized critical state 
are an intriguing possibility. 

10. The analysis of cumulants in terms of invariant mass, or related variables, reintro- 
duced recently after early work, has helped in clarifying several issues in intermit- 
tency. The reasons behind the dramatic failures of models for hadron-hadron col- 
lisions are clearly revealed. They are in first instance incorrectly predicted particle 
and resonance production rates and the near absence of correlations (or even pres- 
ence of anti-correlations) in identical particle systems with small invariant masses. 
These defects are not easy to cure in a consistent manner by simple parameter tun- 
ing and "new" physics may be needed to restore internal consistency in e.g. string- 
fragmentation models. 

11. The hadronization mechanism in hadron-hadron collisions is based on identical phys- 
ical principles and Monte-Carlo algorithms as those applied with apparently great 
success to "hard" processes, in particular to e + e~ annihilation. There is now grow- 
ing experimental evidence that the mentioned discrepancies, first seen in hadron 
collisions, are also present in such processes. These should first be removed before 
commonly expressed claims that "No new physics is involved in intermittency" can 
be accepted. 

12. The relative importance of perturbative and non-perturbative QCD contributions 
to hadron correlations in e + e~ remains controversial, also theoretically. This can be 
studied further with the new techniques now available. 



Chapter 5 

Theoretical description 



In parallel to the extensive experimental effort in quest for power-law behaviour, intense 
activity has developed on the theoretical side, to find acceptable explanations of the 
rapidly accumulating collection of data on factorial moments. The meaning of "inter- 
mittency" in multiparticle processes is still the subject of much debate and no definite 
consensus has as yet emerged. Let us remember once again that we are dealing here with 
the problem of evolution of particle number distributions (or multiparticle correlations) 
in ever smaller bins. 

A priori, the most direct road of attack starts from quantum chromodynamics (QCD), 
now firmly established as the theory of strong interactions. Unfortunately, since the prob- 
lem of confinement is unsolved, QCD can only be used as a guideline to build phenomeno- 
logical models for soft hadronic phenomena. While successful for e + e~ -annihilation, such 
models remain at present unsatisfactory for most other processes, and in particular for 
hh-collisions. The model's deficiencies are often invoked in support for claims of "new 
physics" , but also this matter is far from being settled. 

From the outset it is clear that phenomena such as "intermittency" are manifestations 
of dynamics in a, most probably, strongly non-linear regime of QCD. It is, therefore, 
quite likely that the observed phenomena are not very sensitive to the precise form of 
the Lagrangian, even though the general properties of the interacting fields (e.g., the 
vector nature of gluons) surely play a crucial role. Hence, a satisfactory description might 
be possible on the basis of quite general properties of non-linear systems as revealed 
by complex systems in many other branches of physics. This idea lies at the origin 
of various attempts to establish connections with models for turbulence, multiplicative 
cascade processes, "effective" field theory, the statistical mechanics of disordered systems, 
fractals, phase transitions of various kinds and others. 

A perturbative QCD approach to intermittency would provide a viable explanation if 
the hadronization of quark and gluon systems possesses the property of "early confine- 
ment" such that local parton-hadron duality would hold. Various efforts in this direc- 
tion have indeed established that parton (quark-gluon) avalanches exhibit (multi)fractal 
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properties. These follow from the (fortunate) fact that the process possesses Markovian 
properties. Further developments along this line will evidently improve our understanding 
of perturbative cascades. A direct connection with experimental observations is, however, 
not yet established. 

Numerous other interpretations of "intermittency" have been advanced and will briefly 
be described further below. 



5.1 Simplest approximations 
5.1.1 The linked-pair approximation 

If the rapidity distribution does not change appreciably within a bin interval (this is 
justified for small intervals, at least) one can rewrite ( [2.71| ), approximating it by 



1 M r 

K ^ 5y) ^Mi5vV^ /nWG/r, -..,!/,). (5-1) 

1 y) m=1 5y i 

Since there are no statistically independent contributions to the cumulant functions C q 
and K q , their arguments should be somehow linked. Studying the correlation of galax- 
ies |[L84|| , it was noted that K q can be decomposed into sums of products of two-particle 
correlation functions, Ki with overlapping arguments, in such a way that all multiparticle 
correlations are expressible as successive two-particle ones, so that the whole chain of 
particles becomes correlated. The last condition is necessary since we have learned that 
g-particle correlations (q > 2) are indeed present in the data. 



In this scheme, higher-order scaled cumulants [ |107| | are written as 



#3(2/1, 2/2, 2/3) = #2(2/1, 2/2) #2 (2/2, 2/3) > ( 5 - 2 ) 



3 

,J perm 

#4(2/1, 2/2, 2/3, 2/4) = y| #2(2/1, 2/2)^2(2/2, 2/s)#2(2/3, 2/4) , (5.3) 

perm 

etc. Here, all the permutations of indices 1, . . . , q are summed over; the number of terms 
is equal to the denominator of the factor in front of the sum. The numerator is an a priori 
arbitrary parameter for each order of correlation. 

Even now, the numerical integration in ( |5.1|) is hard to perform. For that reason it was 
suggested [ |107j | to assume translation invariance of the cumulants K g , and to use the strip 
approximation (i.e. instead of integrating over a set of hyper-cubes with linear size 8y, 
one integrates over a strip along the main axis Y and over the differences Q = (2/i+i — Vi))- 



In this way (|5.1|) reduces to a simple but approximate formula [E6f: 



K q « A q (K 2 ) q ~ 1 . (5.4) 
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Substitution of K q in ( |2.72| ) allows any factorial moment to be expressed in terms of the 
second moment so that, for example, 

F 3 = 1 + 3K 2 + A Z K\ . (5.5) 



As mentioned in Sect. 4.5, one can describe |26| the UAl and UA5 data on factorial 
moments at energies from ^=200 GeV to v /i=900 GeV with constant values of A q for all 
intervals 5y. Still, the intermittency indices derived from the linked-pair approximation 
remain somewhat below the experimental ones. At the lower energy of the NA22 experi- 
ment, one gets a much larger value of A% if the above estimate of F3 is used boldly [|142 



However, the assumption of translation invariance is not really justified there and one can 
hardly use ( |5.5| ) on data which are averaged over a large region of phase space. 

Apart from this problem, there are also more basic questions which remain unan- 
swered. For instance, even within the framework of the linked-pair ansatz, one could add 
to ( |5.2| ) a term containing a product of three K 2 s (loops or ring-graphs), as well as further 
terms with multiple links among the pairs. 

Although a dynamical justification for the linked-pair approximation is lacking at 
present, it should be remembered that similar approximation methods have proven their 
utility, e.g. in the theory of liquids. Within the Born-Bogoliubov-Green-Kirkwood-Yvon 
(BBGKY) hierarchy scheme, they allow to "close" the otherwise infinite sequence of equa- 
tions relating correlation functions of all orders. That the linked-pair ansatz may have 
more than accidental relevance is further indicated by the non-trivial fact that ( |5.4|) with 
A q = [q — 1)! corresponds to a multiplicity distribution in Sy which is of Negative Binomial 
type with /c-parameter k = 1/K 2 (5y) [ |109|| . This two-parameter distribution satisfacto- 
rily describes a large variety of (non-averaged) multiplicity data and also occurs as an 



approximation to the soft parton multiplicity distribution in QCD-jets [ |268| - |270f . Further 
extensions of the linked-pair approach beyond those of | 107 , |109| are treated in [[271 , 272 ]. 

The structure of many-particle correlations has also been analysed in partially coherent 
radiative systems ||201]| . This approach is closely related to the linked-pair ansatz ||109 | . 

Conformal theories, treated in connection with intermittency in [ 273j ], provide an al- 



ternative ||274|| to the linked-pair ansatz. In such theories the g-th order irreducible Green 
function is written as a product of two-particle ones to the power l/(q — 1). Taking into 
account the q(q — l)/2 permutations of all particle indices, one finds 

w/2 



K 

IV q ~ 



B q (K 2 



(5.6) 



instead of (|5.4|), which also fits the experimental data reasonably well [274]. 



5.1.2 The singularities of the correlation functions 

According to relations ( |2.68| )-( p.72| ), the singular behaviour of the factorial moments at 
small rapidity binning implies that the correlation functions are singular for small separa- 
tion of their arguments. In particular, the leading singularities of the correlation functions 
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p 2 and C 2 should coincide with the singularity of the corresponding factorial moment F 2 , 
if the formal mathematical limit Sy — > is considered. For 

F 3 ~(6y)-+> (Sy^O) (5.7) 

one should get 

C 2 ( yi ,y 2 ) ^C^(y 1 ,y 2 )\y 1 -y 2 \-^ + C?\y 1 ,y 2 ) (5.8) 

with (3 = 02, Cf ) is a regular function of (y\ — 7/2)1 while C 2 can contain non-leading 
singularities (less singular than the first term). 

A two-component model of this kind has been used in [|275| , where C^ (and similarly 
for the higher-order correlations) is chosen to be a regular function which results in a 
constant additive term to F 2 . The origin of the singular term has been ascribed in |[275 



to a phase transition. For a singular term to be dominant, it must overwhelm F 2 even 
numerically. However, the experimental data discussed above indicate that this is not the 
case. Linear dependence in a double-log plot is observed over quite a large background in 
the available region of 5y. This implies that in this region the integral contribution F 2 L ^ 
of the singular terms in ([5.8D to F 2 is rather small so that 



F 2 (L) « F 2 . (5.9) 
Such a situation provides ||275|, |178|] new possibilities with f3 different from 2 , since in 



that case one gets 

lnF 2 « ln(l + F 2 (N) ) + F 2 (L) /(1 + F 2 (N) ) • (5.10) 

It could even suggest a logarithmic dependence of F^ on by: F^ ~ log by, i.e. a 
logarithmic singularity of the correlation function for coinciding rapidities. However, 
such a behaviour is indistinguishable from a power-like one for small exponents (3 and the 
rather restricted range of rapidity intervals (e.g., 0.1 < by < 1) in which one usually looks 
for this dependence. Actually, if (3 log by <^ 1 one gets ||178|| 

ML) 

lnF 2 « ln(l + F 2 (N) ) + -^—{5 y )-P ^a 2 - <p 2 \nby (5.11) 

Pi 



where 



r (L) ML) 
a 2 = ln(l + F 2 (N) ) + and 2 = 0-^- ■ (5-12) 



' X,; — %r- and d> 2 = 0—^ 
Pi Pi 

Herefrom, one would expect the intermittency exponent to be much smaller than the 
corresponding strength of the singularity of the correlation function, i.e. 2 (3. 

The difference between logarithmic and power-like singularities may become observable 
for higher moments at small by as an upward curvature appearing on log-log plots for 
power-law dependence. This would be more noticeable for smaller Sy, for higher q and 
for larger values of (3. Still, one should not enter into the region of extremely small Sy, 
where the empty-bin effect may dominate and turns down all the curves. 
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The existing experimental data on factorial moments are not in contradiction with 
the above, although no clear signal of an upward curvature of higher moments is seen 
because of large irregularities appearing at small by and large q. These irregularities 
are suppressed if the method of correlation integrals ||130| , |245| , |107| , |T4"3"lj is used, where 
the binning procedure is not fixed but naturally follows the event structure. This has 
been discussed in Sect. 3.8. In fact, the singularities are still better exposed if factorial 
cumulants (|2.71) are used instead of factorial moments [177]. 



5.1.3 Intermittency and Bose-Einstein correlations 

One of the possible sources of increase of factorial moments at small bin sizes is the 
well-known attraction of identical Bose-particles (pions) when their momenta are very 
close. Therefore, one is tempted ||202|| to the extreme supposition that " intermittency" is 



governed by Bose-Einstein correlations, i.e. by symmetry properties of fields but not by 
their dynamics. As was shown in Chapter 4, the experimental data do indeed give some 
indications on the relevant contribution of such an effect. In general, the introduction of 
BE correlations tends to reduce the disagreement between experimental data and Monte- 
Carlo models. However, it was also shown there that the dynamical part is non-negligible 
and, consequently, is of main concern to us. 

Referring the reader to the more specialized review of the subject in |[276|| , we would 



like just to point out that the same physics effects can become more (or less) pronounced 
depending on the corresponding choice of the variable in which it is displayed. In par- 
ticular, it is shown in |276|1 that the rapidity variable is suitable to reveal the dynamical 



intermittency and to suppress the Bose-Einstein contribution to factorial moments. At 
the same time, when analysed as functions of squared 4-momentum transferred between 
pions, the factorial moments get the power-like increasing share of BE-correlations which 
shadows the true (dynamical) intermittency. Therefore, the power-like behaviour in that 
variable neither proves nor disproves dynamical intermittency. One should keep this in 
mind when looking at the corresponding experimental data. Surely, for quantitative es- 
timates Monte- Carlo calculations are necessary with full account of the indefiniteness in 
the description of the BE-effect itself. 



5.2 Dynamical approaches 
5.2.1 Various theoretical models 

From a theoretical point of view, experimental results are best approached via quan- 
tum chromodynamics (QCD). Attempts in that direction are further described in Sub- 
sect. 5.2.6. Unfortunately, the application of QCD to soft processes involving small mo- 
mentum transfers is quite limited since strong non-perturbative effects are involved (unless 
we use additional assumptions, as the local parton-hadron duality hypothesis generalized 
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to correlations of any order). Hence, we are compelled either to construct general relations 
like as of Chapter 2 and the previous section, or to develop phenomenological models that 
fit experimental distributions by adjusting a number of free parameters. 

By now, many phenomenological models have been proposed. Ideas inspired by QCD 
have been used in parton shower models and in their phenomenological counterparts: the 
dual topological model and quark-gluon string models [ff3| , [74], ff6| , |277| - |279|| ; in coherent 
gluon jet emission (Cherenkov gluons, in particular) |280|1 ; in the cold quark-gluon plasma 
model 



192 



Models of a still more phenomenological nature have been tried, such as clus- 
ter models f67] , |28lH283|1 , clan models R284H and narrow hadron jet emission |[L79|| . Whereas 
in all these models definite dynamical mechanisms are proposed for the origin of the fluc- 
tuations in multiparticle production, they still suffer from one important deficiency: they 
do not reveal the nature of the scaling laws observed for factorial moments at small bin 



sizes. 



From that point of view, one would prefer the random cascade models pT|,p2|,p85 



and/or the general approach of phase transitions ||275| , |193| , |286| , |287|| . While the cascade 
models rely on analogies with turbulence theories and lead to phase transitions, general 
considerations of the transition from parton to hadron phases of the process are based on 
important properties of strong coupling field theories reminding of QCD lattice compu- 
tations and the conformal group symmetry. Both approaches lead in quite a natural way 
to scaling behaviour of factorial moments and are preferred as heuristic tools. However, 
to date they cannot compete with phenomenological Monte-Carlo models in providing 
computational results comparable with experimental data at the same level of precision. 
We shall discuss them separately at some length later, together with important ideas 
of intermittency and fractality, but first we shall describe a variety of phenomenological 
models applied to the fluctuation problem. 

First of all, we should mention the furthest developed Monte-Carlo versions of QCD 
inspired models [ |?^[7^ ! [76| , ^ . p5| , p77j - p79|| of parton showers or quark-gluon strings. Qual- 
itatively, these models describe the behaviour of the two-particle correlation function 
C*2 (2/1,2/2) observed experimentally, showing the signature of short-range correlations. 
However, they cannot pretend to fit them quantitatively in hh collisions for all topologies 
and cut-offs. All the models predict noticeably smaller values of intermittency indices 
than the experimental data. Such models also fail to describe the probability distribution 
of maximum particle number per event at a given resolution 5y p8j[289|. 

As discussed in Subsect. 4.2.3, for e + e~-annihilation the situation is still controversial. 
Initially, the DELPHI-collaboration, working at the Z° peak, claimed agreement with 
the LUND parton shower model. Later, increasing statistics tenfold, it finds that the 
agreement only holds at the level of 10-20% (see Fig. 4.7d) ||136fl . However, in general the 
situation is better here than in hadronic reactions. In nucleus-nucleus reactions, these 
models fail to describe the damping of spectator particles observed in experiment [|290|| . 
Further, more detailed studies are needed. 



Thus, we see that fluctuations appear to be a stumbling block for phenomenological 



5.2. DYNAMICAL APPROACHES 



93 



models. They meet with difficulties when confronted with measured factorial moments, 
in particular in hadron-hadron collisions. 

In earlier days, no data existed on factorial moments for small bins and information on 
strong fluctuations in the number of particles inside such bins existed only for individual 
events. A distinctive feature of the fluctuations was the azimuthal symmetry of particles 
belonging to the spike. The whole event showed a noticeable ring of particles in the plane 
perpendicular to the collision axis. Precisely for this reason, the very first attempt to 
explain such fluctuations was based on an analogy with Cherenkov photon radiation. 

The hypothesis of coherent emission of gluon jets [ ESQ ] (involving, in particular, the 
Vavilov- Cherenkov mechanism) predicted that these jets should be emitted in a narrow 
pseudo-rapidity bin at rather large angles in the center-of-mass system of the colliding 
hadrons. All subsequent models did not predict any particular polar angle dependence for 
the dense groups of produced particles. This specific feature was experimentally verified 
in pp-inter act ions at 205 and 360 GeV energies ||291|| . It turned out that the distribution 
of centers of dense particle groups on the pseudorapidity axis contained several peaks 
superimposed on a fairly strong background. Consequently, the proposed mechanism of 
coherent jet emission in hadron interactions probably does exist but is not dominant. 
It could provide the only distinction [2UI] between pp and pp data available from ISR, 
but no analysis of that kind has yet been performed. The ring-like events were observed 
in earlier cosmic-ray experiments ||1 14| - ||1 1 7| and in recent studies of nucleus-nucleus 
collisions 
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Nevertheless, other mechanisms must be involved since intermittency is 
also observed in e + e~-annihilation, where the conditions for coherence seem unlikely to 
be satisfied. 



Large fluctuations may arise in an extended blob of a cold quark-gluon plasma [|192| . 



The appealing feature of such a model is the relationship between this phenomenon and 
the production of soft and 1ow-j>t hadrons, lepton pairs and photons. However, this 
model faces problems in explaining the large values of intermittency indices in electron- 
positron annihilation compared to hadronic and nuclear processes, since, contrary to 
present experimental results, it leads one to expect that the effect is largest in nucleus- 
nucleus collisions. 



Models based on clusters [|282|.|283f or clans ||284|| are more flexible in fitting factorial 



moments, since they involve several free parameters. It has been demonstrated [192, 293 



that the existence of clans leads to a power-law increase of factorial moments for smaller 
bin sizes. No quantitative comparison with experiment has been attempted, however. As 
to cluster models, in some cases they succeeded in describing the multiplicity distribution 
in symmetric rapidity bins of various sizes [57, 281 ] and of the two-particle correlation 
function. The data available at that time were limited to rather large rapidity interval 
sizes (Sy > 0.5). In smaller regions, however, simple cluster model fail completely (see 
Figs. 4.4b, 4.20c and 4.23a). 

Multiparticle production is described somewhat differently in ||294| , [201|| , where particle 
emission is explained by two types of sources - chaotic and coherent. 
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Some of the above approaches attempt to describe multiplicity distributions in varying 
rapidity bins in terms of the negative binomial distribution (or modifications thereof). 
In the cluster model, this is accomplished by varying the cluster parameters. In the 
clan model, the negative binomial distribution is obtained by compounding a Poisson 
distribution for the number of clans with a logarithmic distribution for their decay. In the 
statistical model with two types of sources, the negative binomial distribution naturally 
describes the chaotic sources, while the coherent sources contribute a Poisson component. 

Even though the negative binomial distribution can be phenomenologically used to 
fit experimental data in the very first approximation, there are definite distinctions from 
it in experiment. Besides, the asymptotic QCD predictions disfavour it revealing new 
features ||295|| of multiplicity distributions. In particular, the NBD cumulants are always 
positive, while perturbative QCD predicts negative values (and oscillations) of the higher 
order cumulants. The pQCD prediction is supported by experiment for the case of the 
total rapidity range. 



5.2.2 Intermittency and fractality 

In most cases, the models considered above need to have their parameters adjusted to 
be able to fit (if at all) the data on factorial moments. Power-law behaviour of factorial 
moments is most naturally obtained for cascading mechanisms and in phase transitions, 
as we shall see later. 

The concept of intermittency has been borrowed from the theory of turbulence. There, 
it represents the following property of a turbulent fluid: vortices of different size alternate 
in such a manner as to form a self-similar structure. They do not fill in the whole volume, 
but form an intermittent pattern alternating with regions of laminar flow. Mathematically, 
this property is described by a power-law dependence of the vortex distribution moments 
on the vortex size. This is the reason why the exponents <p q in the power-law dependence 
of factorial moments F q (5y) oc {5y)~^ q in ( |2.110p are called "intermittency indices". 

As mentioned in Sect. 2.4, the self-similar nature of vortices directly implies a connec- 
tion between intermittency and fractality. Fractals are self-similar objects of a non-integral 
dimension. The fractal dimension is a generalization of ordinary topological dimensional- 
ity to non-integers. 

More complicated self-similar objects exist, consisting of differently weighted fractals 
with different non-integer dimensions. They are called multifractals and are characterized 
by generalized (or Renyi) dimensions D q which depend on the rank q of the moment of 
the probability distribution over such objects. The analysis of multifractals according to 
the Levy indices goes beyond the simple definition of intermittency. 

The formal definitions are given in Sect. 2.4. For more details, we refer to the review 



papers [44, |5, 47 1 and references therein. In connection with multiparticle production, 



fractals were first mentioned in [124, 126-128 
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Somehow, the concept of fractality goes beyond a purely formal definition of inter- 
mittency, by connecting the observed dimensionality with the geometrical and thermody- 
namical properties of an object, as well as with the properties of the distributions over 
this object fll"4] , 130|| . The power-law behaviour of factorial moments reveals the fractal 
structure of rapidity distributions as decomposed in individual events. Its relation to ge- 

the 



ometrical and thermodynamical properties is discussed below. In string models [^T 1 
self-similar behaviour of fluctuations is ascribed to the fractal nature of the phase space 



available for subsequent branchings that is formulated in [296 ] as a "plumber view" of 
multiparticle processes. 



Sometimes, intermittency is ascribed [|297|| to the fluctuations of the geometrical sizes 
of emitting sources. 

Some caution is necessary when applying these concepts to multiparticle production. 
As discussed in Chapter 4, finite statistics, the rather small number of particles produced 
in an event and, especially, in a given cell, the method of bin-splitting, the rather restricted 
range of bin widths over which the power-law behaviour is observed, all influence the final 
conclusion. This is described in more detail in [|5LE58]. 



5.2.3 Random cascade models 

In turbulence, intermittency was first demonstrated in cascade models ||299|| . Modifica- 



tions of these, as applied to multiparticle processes, are rather popular nowadays |21|,[22 

In such models, one considers a series of self-similar steps in partitioning phase-space. 
Let us denote by M the number of bins obtained by breaking up the total phase space 
into A parts at each of the v iterations of the self-similar cascade. Thus M = X u (= 
for a total rapidity range AY divided into bins of width 8y). Random cascade models 
involve a probability distribution r(W) with corresponding moments 

(W q ) = JdW r(W)W q , (W) = 1 . (5.13) 

The function r(W) induces density fluctuations as the rapidity window is broken up into 
ever smaller bins. The density P m in the m-th bin is given by the product 

P - = T,fl W ^hf M -y ( 5 - 14 ) 

M n=l M (PM/ 

where the sequence of indices n defines a path leading to a given bin m with density 
P( m )- One assumes that there exists a range of scales inside of which the weights W are 
constant, i.e., they do not depend on the scale at which they operate. 

Herefrom, the intermittent character of the models follows as: 

F q = ((MP m y) = (fl Wl) = (AY/Sy) 1 ^"^ . (5.15) 

n=l 
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The intermittency indices are equal to 



<p q = \n(W q )/ln\ 



(5.16) 



i.e., random cascade models possess a multifractal spectrum [44 



The simplest type of distribution r(W) is the subclass of so-called a-models ||300|| given 
by the two-level probability distribution: 



r(W) = pS(W - W_) + (1 - p)5(W - W+ 



(5.17) 



where < W- < 1 < W + andpW / _ + (l— p)W + = 1 because of the normalization condition 
( |5.13[) . The density enhancement W + > 1 occurs with probability (1 — p) at each step of 
the cascade, while a depletion W- < 1 is present with a probability p. Combined, they 
create "spikes" and "holes" in the rapidity distribution. The intermittency indices are 
given by 



ln[pWl + {l-p)Wl)/lnX 



(5.18) 



The study of moments and multifractal analysis reveal new interesting features. (Let 
us mention that the a-model is reduced to the /3-model for W- = and describes a 
mono- fractal in that case). As the parameters p and A are changed, the model predicts 
various phase transitions | 



The moments of factorial moments are useful 
to reveal these transitions due to the fact that the distributions of factorial moments are 
extremely irregular by themselves [ 302 1. Introducing the normalized moments of moments 



and ascribing to these a power-law behaviour at small bins of the form 



(Z p q ) 



(F q (6y)) p 



oc (6y)*+*(F>(6y))<x{6yy» 



(5.19) 



one can analyse, in the framework of a-models, the dependence of the indices e Ptq on the 



parameters of the model. As discovered in [301], this dependence defines four regions in 



the parameter space, which are reminiscent of four different phases, 
as order-parameters. 



The indices e PA act 



The same conclusions have been obtained when studying [ |194|| the normalized factorial 
correlators F pq /F p F q . Another important property of these correlators is their indepen- 

This property has been confirmed by experiments (see 



21 



dence of the bin width 
Sect. 4.6). However, the a-model does not predict the correct dependence on the distance 
between bins. It predicts power-law behaviour with an exponent 



J pq 



J p+q 



(5.20) 



related to the usual intermittency indices. The experimental values do not follow a straight 
line on a double-log plot. Moreover, there are no finite intervals where they satisfy the 
above relation. When roughly approximated by a straight-line fit, the experimental values 
of (f) pq are larger than those of the a-model. 
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However, one should keep in mind that this is a toy-model, which could pretend to be 
valid for asymptotically long cascades, i.e. for extremely high energies and multiplicities. 
Otherwise, one should develop Monte-Carlo programs | |29£| | losing the beauty of analytical 
formulae. In fact, it has been clearly shown for such a well-known mathematical model 
as the Cantor set [303], that it is not an easy task to reveal its fractal dimension (known 
a priori) using factorial moments, if the number of iterations is finite. 

Nevertheless, the heuristic value of a-models in describing qualitative features of the 
process is rather high since, in particular, they suggest the possibility of phase transitions. 
The nature of the transitions and their relation to the quark-hadron transformation are not 
clear yet. An interesting observation is the existence of "non-thermal" transitions similar 
to those in spin-glass systems. They differ from the usual "order-disorder" transitions by 
producing "different order" in different regions of phase space, so that one may call them 
"clustered order-disorder" transitions. In that case, the intermittency indices increase 
with increasing rank faster than linear. Analogies to statistical-mechanics systems [1 193 
194j| provide further insight into the nature of the transitions. 



5.2.4 Field-theoretical approach and phase transitions 

Besides the scenario of a self-similar cascade, as another extreme, a higher-order quark- 
hadron phase transition has been proposed to explain strong fluctuations leading to inter- 
mittency patterns ||275| , [287j| . Evidently, the statistical mechanics description is most use- 



ful here. Hadronization of a quark-gluon plasma becomes the origin of exceptional events 
with large fluctuations observed above a strong background of conventional events. Such a 



point of view is supported by studies of a two-dimensional lattice of Ising spins ]|304j ,pi 
which shows that intermittency appears at the critical temperature. The intermittency 
indices are directly related to critical exponents of the system. Similar features have 



been found for the g-state Potts model on the Bethe lattice ||305|| at the phase transition 



but without long-distance correlation. Clear fractal structure is also observed for SU{2) 



gluodynamics near the phase-transition point in lattice calculations [ |306| . 

The scenarios of cascading and of phase transitions need not contradict each other, 
if one accepts the point of view that the role of a quark-hadron transition is to fix the 
fractal pattern formed by cascading. The fluctuations are "frozen" at the transition point 
and can be computed by just considering this point. 

A well-defined field-theoretical procedure exists to treat fluctuations and phase transi- 



tions in common media [307 ]. It requires, first of all, the definition of an order-parameter 
and its treatment as an effective fluctuation field. In case of multiparticle production, one 
could choose the rapidity density distribution of particles in individual events p( e ) (y) (or a 
function of it) as an order-parameter which fluctuates about its inclusive average p(y) at 
each rapidity value y. Its function as an order-parameter is clarified by the local parton- 
hadron duality hypothesis, which has been successful in describing the experimental data 
for electron-positron annihilation processes. This hypothesis states that the average val- 
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ues of p( e ) (y) at partonic p( p ) and hadronic p^) levels differ by a (for all rapidities) common 



numerical factor. In particular, if one defines [287] the fluctuation field as 



P(e)(y) 

p(h)(y) 



(5.21) 



then its average in the hadronic phase is equal to zero, while it differs from zero at the 
partonic level: (£(y))( p ) = P(p){y) / P(h){y) — 1 = const. Other possible choices for the 
fluctuation field exist (for example, p 1 ^ 2 (y) |[275| , [286| , |308|| or p — pqa |[273|| ). They have 
advantages and disadvantages which we shall not discuss here. 



The probability of a fluctuation is given by 

W(e) = Z _1 exp[-.F(e)] 

Z = jDe exp[-F(e)] 



(5.22) 
(5.23) 



where TJe refers to the functional differential, F(e) is the free energy and Z is the partition 
function. Adding to F(e) a term J(y)e(y) with an external current J{y), one obtains the 
irreducible Green functions: 



(e 1 ...e g ) 



5 q InZ 



I j=o 



5J{y x )....5J{y qJ 

which are related to correlation functions, so that, for example (we omit the index h) 



(5.24) 



(£l£2)= «%f) -l = K 2 (y u y 2 ) 



p{yi)p{y2> 

The factorial moments are easily obtained as 

fSy rSy 



(5.25) 



where 



roy roy 

F g (5y) = (Sy) 1 -* J Q d&... J q dC 9 -ir 9 (d, . . . , C 9 -i) 



& = yi+i-yi and r g = p q (y 1 ,...,y q )/p(y 1 )...p(y q ) 



(5.26) 



(5.27) 



At first sight, the fluctuation field theory is not directly related to the underlying 
QCD. Yet, these theories are connected through the fluctuation pattern of individual 
events pu\(y), which should be described by both of them if they pretend to be valid. 
Thus, our guesses on the fluctuation field e{y) reflect special features of cascading and 
confinement in QCD. 

For small fluctuations, one can represent F(e) by a Taylor series 



F(e) = F + Jdy 



d£ 

dy 



1 2 



+ -e 2 + ce 3 + de A + . . . 



(5.28) 
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which corresponds to the Ginzburg-Landau Hamiltonian when c = 0, d ^ and all higher 



terms are equal to zero. It has been found ||188|| that some scaling indices (but not the 
critical exponents) have universality properties in this approach. 

For free fields, i.e. c = d = . . . = 0, one gets 

(eie 2 )/ = 7exp[-|yi-Jfe|/f] , (5.29) 

7 = 7r£/&; i={b/a) 1 ' 2 . (5.30) 
This exponential form fits the two-particle correlation function qualitatively (and is often 



used, in particular, for nucleus- nucleus collisions ||217|| ), but it does not provide intermit- 
tency at small 5y. One should remember that it is related to the free-field Lagrangian but 
not to the Ginzburg-Landau potential and, therefore, describes usual short-range correla- 
tions without any phase transitions. One should also note that the approach is formulated 
in momentum space and not in configuration space. 

One is tempted to conclude that fluctuations are strong at small rapidity intervals 
and that the perturbative approach fails. The phenomenon of intermittency should be 
described by a strong coupling field theory, where perturbative methods do not work. 
In particular, the renormalization group approach and conformal theories have been 



tried |[287|, |273|| and have provided power-law behaviour of Green functions and facto- 



rial moments at small bin widths. So, it seems rather reasonable to fit the correlation 
function by an expression 

C 2 (j/i,j/ 2 ) oc 1 — |-exp[-|yi-y 2 |/C], («<1) (5-31) 

for all rapidity separations. For rapidity separations smaller than the correlation length £ 
one gets pure power-law dependence of the correlation function due to a phase transition 
phenomenon. The associated singularity should soften energy and transverse momentum 
spectra. In the simplest approximation, intermittency indices increase linearly with their 
rank. 

Let us stress an important difference of the above consideration with the previous 
treatment of phase transitions. The correlation length £ does not tend to infinity and 
the exponential law is not replaced by a power law at large Sy, as one is accustomed to. 
Instead, the power-law appears at small Sy (6y <C £), and does not influence the depen- 
dence at large rapidity. This happens, because rapidities play now the role of coordinates 
in the usual treatment, so that one has to deal with the ultraviolet (not infrared) stable 
point of the Gell-Mann-Low function. One can speculate that particles lying far apart 
on the rapidity scale reveal the dynamics of the process with a finite correlation length 
related to a particular form of the Lagrangian, while those at nearer points remind of the 
self-organising critical processes with a scaling law not tightly connected to a particular 
form of the Lagrangian (sandpile phenomenon). Thus, correlations appear as " frozen (due 
to hadronization) sounds" of cascading. 
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Similar problems have been discussed in the framework of Feynman- Wilson fluid mod- 
els ||27Ej , |286fl . One should introduce the notion of temperature, additional assumptions on 
thermal equilibrium, on Kadanoff scaling at the critical temperature, on the relative role 
of conventional and stochastic (or quark-gluon plasma) components, and so on. Imposing 
special boundary conditions | |309| | on the grand-canonical partition function, one can re- 
late Kadanoff scaling in the fluid to KNO-scaling in multiparticle processes and describe 
the fractal properties of the fluid in a wide range of scales. Formula (|5.31 ) appears to be 



valid for correlation functions of the conventional hadronic system, but for a system at 
the critical point, pure power-like behaviour with a different exponent is restored. One 
is therefore lead to consider the whole process in the framework of a two-component 
(conventional + critical) model. 

The same approach has been extended [5TH] to a multidimensional analysis of in- 



termittency using, however, the assumption that the correlation functions factorize in 
rapidity and transverse momentum. The predictions for factorial moments differ from 
QCD predictions. The factorization hypothesis allows to proceed analytically and to re- 
late intermittency to fractal properties of the system in original space-time (this problem 
has been addressed also in |[L28] , |311 , |165 |) but looks rather artificial for any field theory 



(QCD included). For instance, if one assumes that conformal symmetry is responsible for 
intermittency [[2871 , [273|j , one obtains non-factorizable Green functions and the predictions 



differ from the above-mentioned ones due to the mixing of longitudinal and transverse 
momentum components, inherent in field theories. We shall see, however, that conformal 
theory and QCD also differ. 

Numerical values of intermittency indices can be calculated in the conformal scheme 
and agree qualitatively with experimental findings. Again, the phase transition plays a 
crucial role. 

Studies of the role of phase transitions in multiparticle production are still in their 
infancy and have, until now, provided qualitative results only. Also the relation between 
hadronization and phase transitions in the simple cascade models treated in the previous 
sub-section is not yet clear. 

The relative role of parton cascading and hadronization is another matter of debate. 
The problem would be solved if extreme proposals were valid. Indeed, parton cascad- 
ing with an "infinite" number of steps would provide intermittency indices quadratically 
increasing with their rank (this corresponds to "wild" or "hard" singularities). Phase 
transitions, on the other hand, yield monofractal behaviour with a linear increase of the 
indices. In reality, the two extremes must be modified so that finite cascading would lead 
to a slower increase, while the next operator-product expansion terms for a phase transi- 
tion would induce a faster than linear rise of the intermittency exponent. Such problems 
have as yet not been treated. 
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5.2.5 The statistical-mechanics formalism 



Statistical analogy is a powerful way to analyse properties of chaotic dynamical sys- 
tems |pi2|1 , in general, and of multifractals and cascade models, in particular ||313| , |193| . |45f . 
It rests on the possibility to define a partition function Z(q) of the system in the following 

way: 

M 

Z(q) = E Vl 



m ) 



(5.32) 



m=l 



where 



Pm = -r^~ (5.33) 

is a normalized probability weight on the ensemble of bins (m) and p( m ) is a random 
(rapidity) density registered in each bin m. 

The relation to multifractal (or intermittent) properties of a system is established if 
one considers systems for which the probability inside a box m is proportional to a power 
a m of the box size and the number of degenerate boxes (with the same value of a m ) 
follows a power law as well. Then, assuming a continuous limit, one finds as in ( |2.118| ) 

a+ 



Z{q) 



M s{a) - qa da 



(5.34) 



(integration running between maximum and minimum zeros of /(«)), wherefrom one 
easily obtains the multifractal spectrum of the system f(a) (see, for example, [44, 47, 1291 ) . 

Remembering the definition of the intermittency indices 
spectrum f(a) via the relation 



one relates them to the 



-q 2 - 
dq 



+ 1 



where a is defined as 



a{q) 



dcpq 
dq 



(5.35) 



(5.36) 



The interpretation of f{a) is transparent since it weights the number of degenerate 
boxes. It, therefore, corresponds to the entropy in statistical mechanics. In a similar way, 
the rank q may be interpreted as an inverse "temperature" (3 = 1/T and the relation 
Q5.35D corresponds to the usual thermodynamical formula 

S—%, (5.37) 

where S is the entropy and F is the free energy, whose "temperature" dependence is 
provided now by 



A(ff) 



, + 1 
q 



l-F . 



(5.36 
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Two features of this analogy are particularly useful. On the one hand, application of 
the thermodynamical formalism allows for coverage of multifractals by boxes of different 
sizes, i.e. for a more precise description of individual events. This has been used in 
proposals of correlation measures with non-uniform coverage ||130| . |245| . |143|| . On the other 



hand, the multifractal treatment admits an extension of the thermodynamical formalism 
to non-equilibrium systems. It has been used to classify the phase transitions in a-models 
and to demonstrate that, in multiparticle processes, phases could exist similar to spin-glass 
states || ,§. 

It is important to note that the minimum of the X(q) (|5.38|) corresponds, according to 
( |5.35| ), to zeros of the fractal spectrum. This is a signal for a phase transition in thermo- 
dynamical systems. One should stress, however, that the similarity of the distributions is 
in itself not sufficient to justify use of statistical physics terminology in its original mean- 
ing for the quantum field systems we are interested in here. Besides, the analogy breaks 
down completely for values of q exceeding the multiplicities effectively contributing to the 
moments. Nevertheless, this analogy is used in ||314| , |315|| to derive the dependence of the 



pressure in a Feynman- Wilson liquid on its chemical potential and some peculiar features 
are found in hadron-hadron reactions (see also the review paper ||316|| ) . 



It is evident that further explorations of statistical mechanics approaches to multipar- 
ticle production are needed. In particular, analytical properties of a partition function, 
often useful in connection with phase transitions, have not been much analysed. 

The location of the (complex) roots (zeros) zv of the multiplicity generating function 
( |2.43[ ) has recently been studied in | |317| , |318| | after earlier suggestions by Biebl and Wolf 



This work is based on the analogy with the famous Lee- Yang zeros [ |319|| , whose location 



fully characterizes the thermodynamic properties of the physical system. 

For multiplicity distributions, the strength of the fluctuations of the multiplicity in an 
event is directy related to the location of the zeros in the complex z plane: the magnitude 
of the factorial cumulants, and thus the strength of the correlations, is determined by the 
roots closest to the origin. 

In the discrete version of QCD, developed in LUND, is was demonstrated that the zeros 
of (|2.43|) belong to a fractal Julia-set [|320|| with intruiging properties. Detailed studies 



of this set, and various connections with standard phenomenology, such as KNO-scaling, 
remain to be worked out. 



5.2.6 Intermittency, evolution equations and QCD 

In the previous section we considered two rather extreme possibilities, simple cascade 
models and phase transitions, as possible mechanisms leading to scale-invariance in par- 
ticle production processes. 

Further interesting results have been derived from studies of simplified kinetic branch- 
ing evolution equations for "birth-death" (or "gain- loss") processes. Many of these are 
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treated in textbooks ||322| . |323| 1 and were applied to multiparticle production ||324| f328|] . 
In general, the time-evolution of the number of "clusters" (partons, resonances, etc.) is 
described by forward or backward (retrospective) Kolmogorov equations, which relate 
the time derivative of the generating function to some combinations of that function. A 
particular example is the Smoluchowski equation treated in [ 325 ]. The terms "forward" 
and "backward" imply, that each tree graph may be viewed either as a splitting to ever 
"thinner" branches, or as a convolution to ever "thicker" branches. 

For linear evolution equations, the solutions depend directly on the initial condi- 
tions ||205|| . Non-linear equations often have solutions asymptotically independent of 
pre-history ||327| , |329[| . Stationary regimes may appear if the annihilation of clusters is 
stronger than the "birth-rate". In that case, the dispersion is proportional to the aver- 
age multiplicity and intermittent behaviour can be obtained, if the proportionality factor 
is larger than one and the mean multiplicity decreases correspondingly for smaller bins. 
This can easily be proven by means of the definition of the second-order factorial moment. 
Let us note that systems obeying non-linear evolution could exhibit quite general prop- 
erties, independent of the detailed form of the equations, such as period-doubling. Some 
ideas along this line of thought, using properties of stochastic systems and Feigenbaum 
attractors, have been formulated in [ 330 , p31 . 



More detailed analysis of intermittency in the framework of the Smoluchowski equa- 
tion [|329|1 reveals various regimes of time-evolution and cascading, depending on the pa- 
rameters of the model. The Smoluchowski equation is of the backward type. It contains 
terms linear and quadratic in the generating function G with opposite signs. Formally it 
looks like 

( ^L = G*G-G*G 1 , (5.39) 

where G(u,t) = Y^ n >iN(n,t)u n , G\ = G(l,t) with N(n,t) representing the number 
of clusters of (integer) mass n at time t and the convolution * is defined through the 
aggregation coefficients of clusters. The fractal properties of aggregates and the occurrence 
of phase transitions have been analysed in |329|| . 

Obviously, it would be desirable if an explanation of scaling phenomena in multiparticle 
production could be derived from, so to say, first principles, i.e. in the framework of 
QCD. Multiparticle production in QCD is the result of quark-gluon branching and the 
subsequent transition to hadrons. As such, the self-similar multiplicative branching (or 
cascade) process could give rise to a scaling regime. The perturbative QCD parton shower 
picture is justified for interactions with large transferred momenta (or virtualities), but in 
hadronic reactions one mostly has to deal with soft processes. Perturbative QCD is valid 
in the initial stages of high-energy cascades in electron-positron annihilation and could, 
therefore, be used as an explanation for intermittency. 

It is well known that the perturbative QCD cascade gives rise |[332 | to a mean multi- 
plicity of partons increasing rapidly with energy. Equations for higher moments of parton 
multiplicity distributions are rather complicated ||333| , |9~8|j , but reveal in any case that the 
parton number distributions are much wider than a Poissonian. The infrared limit be- 
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comes very important and one should consider infrared-safe properties. Assuming that 
the singularity is avoided in a way similar as in an electromagnetic cascade in a medium, 
one can estimate the fractal dimension of internal motion of partons in a jet and it turns 
out to be quite low for a single jet in e + e~-annihilation [ 12S|| . 



The simplest theoretical models, such as the tree diagrams of the 3 model, simplified 
QCD [|334j . |335| . |29H also based on tree diagrams, the Schwinger tunnelling transition [|336|| or 



the effective Lagrangian approach [|308|1 , indicate that the totality of all produced partons 
exhibits intermittency. 

It is not yet clear what modifications of QCD equations going beyond the tree graphs 
(so-called Double Logarithmic Approximation - DLA) would fit this region best. An 
approach to that problem has been proposed Q337H in deep inelastic processes where the 
transition from the Bjorken limit to the Regge domain proceeds through some intermediate 
region in which quadratic terms (in the fields) appear and cause some recombination of 
partons at high densities. Here, the evolution equation for the number of gluons xG(x, q 2 ) 
in a hadron with a transverse size g _1 and for small values of the Bjorken x- variable is 
taken to be 

d 2 xG(x, g 2 ) _ asxG{ ^ q2) _ C^ [xG{ ^ q2)]2 ^ ( ^ 40) 



d\n(l/x)d\nq 2 q 2 

where a s is the QCD coupling strength and C is a constant. It is inspired by QCD 
ideas and Regge phenomenology, but has not been derived rigourously. No analysis of the 
intermittency property has been attempted so far. One should note, however, that the 
general effect of such a quadratic damping is to narrow the multiplicity distribution (see, 
for example ||338|| ), which leads to decreasing factorial moments. 

One should, however, not rely on the similarity of equations (|5.39| ) and (|5.40|) , since 



this could be misleading. As is well known ||333| , |339| , |98| 1, the equations for the generating 
functionals for gluon and quark jets in QCD are non-linear, while the corresponding equa- 
tions for the structure functions in DLA are just linear GLAP equations. The generating 
functional for a parton p is given by 

G p (u,V,X,Y) = Till / dx W dx gj U ( X i) V ( X j) W ^n g ( x , x i, x ji Y ) ( 5 - 41 ) 

n q ,n g Wq.Wg. J 

where W is a differential probability to create n q quarks and n g gluons with an evolution 
parameter Y ~ lnlnQ 2 in a p jet. The equations for the generating functionals are 

dG (x Y) r 1 

q g y ] = I dx'P qq (x',x)[G q (x',F)G g (l -x',Y) -G q {x,Y)\, (5.42) 



dG g (x,Y) 
dY 



dx'[P gg (x', x)G g {x', Y)G g (l - x', Y) - G g (x, Y)) 
+n f P qg (x', x) (G q (x', Y)G q (l - x', Y) - G g (x, Y))] , (5.43) 
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where the P's are the corresponding GLAP kernels, n is the number of flavours and the 
initial conditions are such that at Y — 0, G q = u and G g = v for a single jet. 

The status of the equations for generating functionals is not completely clear up to 
now. They are able to reproduce the higher order graphs of the perturbation theory far 
beyond the tree level ||98|| . Their success in predicting the tiny features of multiplicity 



distributions in total phase space (for a review see ||295|| ) encourages speculations about 
their quite general status, with some confinement properties taken into account already 
at that stage. 

A simplified version of ( |5.43| ) with the quark term omitted (gluodynamics) has been 
studied in ||339|| . Intermittency for the factorial moments of the gluon cascade is claimed 



to be in qualitative agreement with experimental data, as well as evidence for a structural 
phase transition. However, the formula for the intermittency indices contradicts other 
QCD results. This is not surprising, since the generating function technique should here 
be applied to the subjet hitting the bin under investigation, not to the whole jet as done 
13391- 



in 



The treatment of QCD cascades has been taken further in [ 253|| within the framework 



of the dipole formalism including coherence effects. The multifractal dimension of the 
parton cascade for high order q is found to be equal to the QCD anomalous dimension 
7o = (6a s /7r) 1//2 and a first pre-asymptotic correction has been calculated. Moreover, a 
geometrical interpretation of the anomalous dimension of QCD is proposed. 

A direct solution of QCD evolution equations has been attempted for the second 



correlator in ||254|| . The behaviour of factorial moments of any rank (as well as of double 
trace moments - see below) in small phase-space windows for e + e~-collisions is treated 
both in DLA and in the next Modified Leading Logarithmic Approximation (MLLA) of 



QCD in |[256|| . Similar results for factorial moments in DLA are obtained in ||255| , [257|j . 



They are closely related to the previously derived formulae of |340|1 . 



In the approach described above, one considers three stages of the process: 

1. the initial quark emits a hard gluon, 

2. the gluon evolves into a jet consisting of several subjets, 

3. one of the subjets hits the phase space window chosen. 

Integrating over all the stages one gets the final multiplicity distribution (for details 
see |[256fl ). For comparatively large windows one can use the fixed coupling constant, 



while for smaller bins its running should be taken into account. 



Fixed coupling QCD factorial moments reveal [[254j , |256| the intermittency phenomenon 
with intermittency indices equal to 

<pQCB(q) = (q-l)(D- ^7o), (5.44) 

in DLA. This formula is valid if, for the D-dimensional analysis with M bins along each 
axis, one defines F g oc M^ 9 '. For large q the indices increase linearly. The term with 
negative sign in the second bracket is proportional to the QCD multiplicity anomalous 
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dimension 70. From ( |5.44| ), one would conclude that QCD prescribes fractal behaviour 
with codimension 



d, 



1 1 QCD 



n 4 + 1 
D 7„. 

1 



(5.45) 



The phase space term D is obviously non-fractal. The 70-term in ( |5.45|) is due to the 
energy dependence of multiplicity and gives monofractal behaviour. The gluon energy 
spectrum contribution, represented by 70/9, gives multifractal behaviour. The next-to- 
leading corrections to 70 also provide g-dependent terms. 

The calculated values should be compared to the slopes in the region 5y > 1 (which 
are rather large) since ( |5.45| ) is derived in fixed-coupling QCD (the intermittent behaviour 
in that region was discussed first in ||341|| ). In smaller bins, the QCD running coupling 
becomes important and modifies the above relations | |256| | . The factorial moments now 
behave in a semi-power-like manner so that there is no strict intermittency even though 
an approximate one can still be claimed. The second term becomes very close to 1 for the 
low-rank moments and the low-rank intermittency indices turn out to be very small for 
one-dimensional analyses, in accordance with experiment. Corrections to 70 of the order 



of 7q can be taken into account [256]. The influence of confinement seems negligible. Thus 



small as well as large <5y-intervals may be described in a unified treatment, at least at a 
qualitative level. The transition point from "large" to "small" bins depends on the rank 
of the moment in QCD, in full accordance with experimental findings. 

We should further mention that the ratio dqjd^ depends explicitly on D, contrary to 
experimental claims (see Subsect. 4.3.3). 

The "free energy" F(q) is related to X(q) given by ( |5.38| ) as 



F(q) 



X(q) « 70 



To 
q 2 



(5.46) 



In DLA, it is a steadily increasing function of q. However, with corrections to 70 taken into 
account |256|1 F(q) becomes non-monotonic. This happens for rather large values of q, 
with the result, that the distinction between factorial and usual moments becomes crucial 
and statistical analogies inapplicable. These findings show the limitations of perturbative 
QCD and provide further insight into the properties of multiplicity distributions, such as 
KNO-scaling, in full phase space ||342| . |343| . 



In particular, QCD gives rise to the prediction ||343| , |295|| of a negative value of the 
cumulant of rank 5 confirmed by experiment and to the general conclusion of a non- 
infinitely-divisible nature of total multiplicity distributions in QCD (that prohibits, e.g., 
the one-ladder multiperipheral cluster models). 

The increasing branch of the multifractal spectrum f(a) may be easily calculated 
using ( |5.35|) , ( |5.36| ) giving 

7o) 1/2 - 



/(«) 



o 1/2/ 

27o (" 



(5.47) 
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The double trace moments (DTM), redefined in analogy with factorial moments as 

j? _ 1 I ST n m( n m - l)-(w m - U + 1) V . . 



behave ||256|| in QCD as 

F v>q (A) oc A 2 *^ 11 oc A q ^-^ +q - 1 = A-Kirt+i- 1 , (5.49) 



wherefrom one finds 



4>{q, v) = <j>{qv) - qftv) = (q - 1) (l - • ^ j . (5.50) 



The second factor in ( p. 50 ) may be called "double codimension" . It is not symmetric in 
q and v and shows that increasing v one decreases effectively the anomalous dimension. 
For v = 1, as required, the double codimension becomes equal to the usual codimension. 
The scaling exponent ( |5.50| ) is not factorizable in v and q. The above redefinition of DTM 
is aimed at reducing the Poissonian noise and the role of phase-space factors, otherwise 
very important. 

In fact, it is surprising that the above expression describes qualitatively the general 
trends and even the absolute normalization of the functions K(q, v) shown in |5(|[5l|]. For 



large u, the ratio K(q, u)/(q — 1) is completely determined by the phase space factor and 
should tend to 1. This is seen in the experimental data. In the region of small q ~ 1 and 
v ~ 1 the strong compensation in fl5.50| ) prevents its use, but even there it gives quite 
reasonable values of K(q,u). This probably indicates that DTM defined as powers of 
multiplicities are not sensitive enough to dynamics, a suspicion raised earlier (Sect. 4.7). 



On the other hand, the parton cascading picture employed in [256] may be applied, 
strictly speaking, only to hard processes at extremely high energies and one should not 
rely much on the asymptotic estimates of QCD when considering experimental data. Also, 
the difference between usual and factorial moments becomes extremely important at large 
values of q or at sufficiently small bins with low multiplicities. One may suspect that 



relation ( [5.50 ) and, especially its first part, has a much wider range of applicability than 



just for e + e~-collisions and some universal relation could be valid for other reactions. If 
so, it will be important to understand whether these common features are due to common 
dynamics or to insufficient sensitivity of the proposed measures. 

The large fluctuations, e.g. those observed by the NA22 collaboration, have raised the 
suspicion that at small bins one is dealing with unusually wide distributions, which could 
have infinite moments. This has led one to consider Levy-stable probability distributions. 
The Levy indices derived from QCD factorial moment indices or DTM exponents show 
||256|| no sign of "wild" singularities. 

The DTM technique has been first applied to experimental data from hadron-hadron 
reactions (described in Subsect. 4.7.7) where direct QCD arguments are invalid. For 
e + e~-collisions, however, it may be worthwhile to analyse the data with this method. 
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Let us further note that there is a difference between predictions of QCD and those 
of variants of conformal theories considered in | 273|| , or of multiplicative models |ftlO|| . 



For example, the intermittency indices derived in a rapidity or azimuthal angle analysis 
should be equal in QCD. This is a consequence of the symmetrical form (in pseudorapidity 
rj and azimuthal angle (p) of the gluon propagator, which can be written as 

k 2 = (pr,! +PT, 2 ) 2 = 4 MT , 2 (sinh 2 V f + sin 2 f), (5.51) 

where pT,i is the transverse momentum of i-th parton. In a conformal theory, the inter- 
mittency indices for the second factorial moment are given by 

fa(5y)=2ri, 2 (M = , (5.52) 

respectively. Here, r\ is the conformal anomalous dimension estimated to be rj ~ 0.07—0.1. 
For multiplicative models, one finds 

fc(5y) = 1-D y , fa{8ip) =D V -1 (5.53) 

with < D y < 1, 1< D v < 2. 

The emergence of intermittent behaviour in solutions of non-linear equations and in 
perturbative QCD encourages further studies along these directions. At the least, they 
hold a promise of further theoretical insight. One should keep in mind, however, that a 
direct comparison of QCD-based asymptotic results with present-day experiments is not 
justified. Some effects revealed by the data seem to be of a different, as yet unsatisfactorily 
explained, origin. 



Chapter 6 
Conclusions 



Developments in physics — and in science in general — over the last decade, have brought 
exciting new discoveries and deeper insight into the dynamics of complex systems. Studies 
of classical and quantum chaos, non-equilibrium dissipative processes, random media, 
growth phenomena and many more have all contributed to reveal the pervasive importance 
of self-similarity, of power-laws and of fractals in nature. Research in these fields is still in 
full evolution and continues to uncover intriguingly simple and often surprisingly universal 
behaviour in complex, non-linear systems. 

The suggestion of Bialas and Peschanski to look for scaling in particle fluctuations 
was one of the first attempts to apply modern ideas and techniques from complex-system 
research to multihadron-production processes. In preceding pages, we have presented 
a critical overview of the impressive amount of experimental and theoretical work this 
proposal has generated since its formulation in 1986. The continuing interest in the field 
testifies of the growing conviction that new avenues need to be explored for progress in 
strong-interaction physics. 

Impressive as it may be, this work has not yet led to final answers concerning the 
fundamental issues. Approximate power-law scaling of particle density and correlation 
functions is now indeed observed, especially in two- or three-dimensional phase-space. 
However, so far it can be explained from an interplay between jet formation and more or 
less "conventional" correlations among identical particles due to quantum interference. 

Nevertheless, as often happens, the detailed scrutiny of data (and detectors) on the 
full variety of collision processes in the search for power behaviour has led to many new 
observations of interest in their own right. It has helped to recognize the importance 
of detailed studies of correlation phenomena at large and small distances in momentum 
space and new sensitive and general techniques have been developed for their analysis. 

Standard hadronization models, all too often accepted as satisfactory, have been ex- 
posed to severe and sometimes even painful tests. Intermittency analysis has revealed 
deficiencies in our understanding of the hadronization process. These defects are not easy 
to cure in a consistent manner by simple parameter-tuning and "new" physics may well be 
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needed to restore internal consistency in, e.g., fragmentation models of the LUND type. 
Present work on this subject starts to provide hints that purely probabilistic treatment of 
the break-up of colour fields has to be supplemented with effects deeply connected with 
the structure of the non-perturbative QCD vacuum. Progress in this direction would in 
itself be ample compensation for the efforts spent on attempting to establish fractality in 
multihadron production. 

Data obtained in the last years have shown the overwhelming importance of correla- 
tions among identical particles in the "intermittent" regions of phase space. This quantum 
mechanical phenomenon, discovered in particle physics in 1959, still awaits satisfactory 
incorporation into present hadronization phenomenology, if it is to be used as a reliable 
interferometric tool, e.g., in studies of quark-gluon plasma formation. 

Theoretical work has developed along a large variety of directions. Fractal properties 
have been discovered in string-fragmentation models. Within the realm of perturbative 
QCD, parton correlations and emergence of power behaviour are now studied with in- 
creasing sophistication. The relevance for present-day phenomenology remains doubtful, 
however. 

The powerful methods of statistical mechanics have been intensively exploited in stud- 
ies of random cascades as well as in equilibrium and non-equilibrium critical phenomena. 
Results of real intrinsic value have been obtained, with potentially interesting applications 
in other fields. 

In search for an explanation of "unusually large" density fluctuations, "intermittency" , 
in analogy with fluid turbulence, has progressively led to appreciate the importance and 
often spectacular manifestation of non-linear strong-coupling dynamics. 

Experiments in deep-inelastic scattering are now starting to probe hadron structure 
in a regime where the perturbative parton-cascade picture becomes blurred. Non-linear 
perturbative evolution and confinement play an increasingly important role in the very 
low Bjorken-x region now accessible in HERA. Present attempts to understand this region 
invoke QCD Reggeon-theory. 

It is intriguing to speculate that a power-law dependence of the low-x parton correla- 
tion functions could manifest itself as "gluonic" intermittency in virtual parton cascades, 
with the occasional creation of regions with very large, or very small gluon density. 
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Chapter 7 
Figure Captions 



Fig. 3.1 Contours of the two-particle correlation function, R cc (yi, y 2 ), from 205 GeV/c pp 
interactions pl |. 



Fig. 3.2 a) The charge correlation function C 2 (r]i,r)2) plotted for pp collisions at fixed 
r]i = versus rj 2 at 63, 200, 546 and 900 GeV, b) the "long-range" contribution C*l and 



c) the short-range contribution C$ ||65 |. 



Fig. 3.3 The semi-inclusive correlation function C 2 (771,772) f° r 34 < < 38 pp colli- 
sions at 900 GeV, compared to the UA5 Cluster MC, PYTHIA and FRITIOF 2.0 |§7|. 

Fig. 3.4 Normalized correlation function K 2 (yi,y 2 = 0) for (CC), ( — ), (+ +) and 
(H — ) combinations in n > 6 pp collisions at 360 GeV/c, compared to predictions from 
single chain LUND and a two-chain DPM [[71]]. 

Fig. 3.5 a,b) Correlation functions C 2 (0, y) and 6*2(0, y) for M + p reactions as compared 

with calculations in FRITIOF ( ), DPM ( -) and QGSM ( ) (non-single- 

diffractive sample), c) Correlation functions Cs(0, 2/2) for M + p reactions as compared 
with FRITIOF ( ) and QGSM ( ) WB. 



Fig. 3.6 Normalized two-particle correlation function K 2 (yi,y 2 ) for pairs of a) oppo- 
sitely and b) negatively charged hadrons produced in muon-nucleon scattering at 280 
GeV/c 0. 



Fig. 3.7 Normalized correlation functions K 2 (0,y) for the M + p non-single-diffraction 
sample and /z + p-interactions at 280 GeV/c (13 < W < 20 GeV) 



Fig. 3.8 a,b) Normalized correlation functions -^2(0, y) for the M + p non-single-diffraction 
sample and e + e~-annihilation at i/i=22 GeV. c) Normalized correlation function K 2 c (yi, y 2 ) 
a t yi = —1 4-0 in the non-single-diffraction M + p sample at 22 GeV and e + e~-annihilation 
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at 14 and 44 GeV flTS 



Fig. 3.9 Normalized correlation functions -^2(2/1, 2/2) for the entire data sample (upper 
plots) and for selected two-jet events (lower plots) in e + e~ annihilation at 35 GeV. The 
left plots show the average over the region — 1 < yi+i/2 < and the right plots correspond 
to the average over < 2/1 + 2/2 < 1- The data (open circles) are compared to the JETSET 
7.2 PS model with (solid lines) and without (dotted lines) Bose-Einstein correlations 

Fig. 3.10 OPAL data on R^i,^) as a function of (£1 — £2), for (£1 + £2) between 
5.9 and 6.1 (IV), 6.9 and 7.1 (V), 7.9 and 8.1 (VI), respectively, compared to 
lytic QCD calculations (the dashed curves indicate leading order QCD calculations for 
A=225 MeV, the three solid curves represent next-to-leading QCD calculations with, 
from top to bottom, A = 1000, 255 and 50 MeV, respectively), and b) coherent parton- 
shower Monte-Carlo models ARIADNE (solid), JETSET (dashed) and HERWIG (dotted); 
c) incoherent parton-shower Monte-Carlo models JETSET (solid) and two versions of CO- 
JETS (dashed and dotted) pL 



Fig. 3.11 Correlation functions Cs(0, y) (a) and Cs(0, y) (b) in M + p interactions at 
250 GeV/c (open circles correspond to non-single-diffractive events and are shifted to the 
right to avoid overlap) |75| . 



Fig. 3.12 Topological correlation functions 0^(0, y) in M + p reactions at 250 GeV/c 
for (H — ) pairs ]75| . 

Fig. 3.13 a) (5| n) (0,0) dependence on n and z = n/(n) for M + p interactions at 250 
GeV/c. The last ( ) point corresponds to n > 16, the last (H — ) and (++) points 



to n > 18 ]75[. b) Ct\\m ~ m\) as a function of l/(n — 1) for pp collisions at ISR 
energies |]3ql . 

Fig. 3.14 Normalized correlation functions -^2(0, y) for particles with all pr, Pt < 0.30 
GeV/c, pt > 0.30 GeV/c as compared to functions (3.2) and (3.3), solid and dot-dashed 
lines, respectively 



Fig. 3.15 Distribution in a) the rapidity gap Ay and b) the azimuthal distance Aip for 
KgKg pairs in pp collisions at 360 GeV/c. c) and d) same for K°A pairs [90|. The lines 
correspond to LUND. 

Fig. 3.16 The multiplicity-averaged angular correlation function 
((n — 1)62 (?7i, (fi] T]2i ^2)) for CC combinations in units of 10~ 3 f86|. 



Fig. 3.17 W(Aip, Ay) for inclusive non-single-diffractive 7r + p interactions at 250 GeV/c 
as compared to FRITIOF2.0 (dot-dashed), DPM (full) and QGSM (dashed) @. 

Fig. 3.18 W (Aip , Ay , pt) as compared with calculations in FRITIOF2.0 (dot-dashed), 
DPM (full) and QGSM (dashed) for Ay < 1 and pt cuts as indicated |88| . 
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Fig. 3.19 The pi dependence of the azimuthal correlation parameter B for h + h pairs 



in pp collisions at 360 GeV/c compared to LUND and DPM JrT| 



Fig. 3.20 Charmed pair azimuthal correlation: WA92 data are compared with predic- 
tions from (a) an NLO perturbative QCD calculation and (b) a model where a parton 
transverse momentum, p T , is added to the NLO perturbative QCD predictions (dotted 
line: (p|) = 1.0 (GeV/c) 2 , solid line: {p 2 T ) = 0.3 (GeV/c) 2 [g§. 

Fig. 3.21 Two-particle azimuthal correlations with respect to the sphericity axis in 
OPAL [97| compared to coherent and incoherent MC models. 

Fig. 3.22 The PPCA distribution for corrected L3 data compared to a) coherent and 



b) incoherent Monte-Carlo models ||102 



Fig. 3.23 Three-particle rapidity correlations a) -^3(0, 0, y) [the lines correspond to the 
ISR results at 31 GeV (full) and 62 GeV (dashed)] and b) K 3 (0,0,y) for M+p interac- 
tions at 250 GeV/c. [the FRITIOF (dot-dashed) prediction is indicated for the charge 
combination ( ), QGSM (dashed) for ( ) and (--+)] |7§. 



Fig. 4.1 a) The JACEE event M ; b) The NA22 event JTT9|. 



Fig. 4.2 a) logi<5 as a function of — log £77 for the JACEE event |^T],^2[ (full circles) 
compared to independent emission (small crosses); b) In F 2 and In F4 as functions of — In 5r] 
for O 16 Em at 200 A GeV (KLM) g4|, c) log F q on 2 log<5r/ and d) 2 log<fy> at 630 GeV 
(UA1) P3| . 



Fig. 4.3 Anomalous dimension d q as a function of the order q, for a) /ip and e + e 



collisions, b) NA22 and UAl, c) KLM [162 . 



Fig. 4.4 a) Slope <p q as a function of order q for NA22 [142], two versions of FRITIOF 
and a two-chain dual parton model, b) F 3 vs. 5i] from UAl [|143|| compared to GENCL, 
PYTHIA and PYTHIA + Bose-Einstein Monte Carlos. 



Fig. 4.5 a) EMC results ||140|| compared to expectations from ||108|| , b) slopes (f) q for 



EMC data as well as Webber and LUND models, c) vA data ||141|| in comparison to LUND 
model expectations. 

Fig. 4.6 Jet evolution: the self-similarity in the parton cascade derives from the simi- 
larity of each step in the evolution [|124j . 



Fig. 4.7 lnF 3 (and lnF 2 ) function of 

HRS pi, b) TASSO p4 



In by (or InM, M = AY/Sy) for a) 



c) CELLO pH and d) DELPHI [T3§ data compared to 



LUND (and Webber) parton shower models . 

Fig. 4.8 In as functions of — In 5y for JETSET 6.3 parton shower at ^=91 GeV at 
the a) parton, b) hadron level, both with cut-off Ql=l GeV 2 , c) d) with cut-off Ql=0A 
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GeV 2 |T66 l. 



Fig. 4.9 a) logF 2 as a function of — \og5f for UAl data ||143|| , b) c) slope 4> q 
function of order q for y, rj, ip or (y and ip) as variables for NA22 | 142|| . 



as a 



Fig. 4.10 The log- log plot for TASSO data using a) 2 dimensional y-(p bins, 



b) one-dimensional y bins, in comparison to a two-dimensional a-model [134, 23 1 



Fig. 4.11 Factorial moment of order q = 2 for 1, 2 and 3 dimensional analysis for a) 
DELPHI P61 , b) UAl gg and c) NA22 |T42| as a function of (log 2 M)/d and (lnM)/d, 
respectively, where M denotes the total number of boxes in a <i-dimensional analysis. 

Fig. 4.12 a) Factorial moment In F2 as a function of (In M) / d from a 3 dimensional anal- 
ysis of negative particles in pAu and central OAu and SAu collisions [|15C 
b) Factorial cumulant K 2 from the same analysis in central OAu collisions ||150|| . 



Fig. 4.13 The second scaled factorial moment F 2 as a function of the number of bins 
M 3 in the log-log scale |L77f for (a) ^p data PEH , (b) vr/Kp data P%| , (c) pAu, (d) OAu 
and (e) SAu data [ |15U| |. Solid lines represent ( fO| ) with parameter values: (a) c=0.025, 
02=0.45, c'=0; (b) c=0.02, 2 =O.45, c'=0.16; (c) c=0.02, 2 =O.45, c'=0.2; (d) and (e) 
c=0.01, 02=0.5, c'=0. 

Fig. 4.14 The factorial moments hiF q for NA22 [|142|| in various variables and dimen- 



sions indicated, as compared to ECCO [ 164 



Fig. 4.15 a) Illustration of the modified power-law behaviour. The lines indicate an 
'eye-ball' fit to the data. Only the data in bins where F 2 varies strongly are used; b) Test 
of the universal scaling law (}4.6|) for lni^ and lnF2 in e + e~, lh, hh and AA collisions as 
indicated. The straight line is adjusted to the e + e~ data and reproduced on the other 



data sets 24 



Fig. 4.16 Ratio of anomalous dimensions as a function of the order q for a) NA22 

c) and d) d q as a function 



results ||142|| and b) earlier results on all types of collisions [24 



of the order q. Continuous lines in c) and d) show the best fits using ( |2.134j ) [c) •— 
Z° decay, DELPHI, r = 1.6(8), ■- p-Ag/Br, KLM, r = 1.1(5), S-Ag/Br, KLM, 
r = 0.0(2); (d) •- O+Ag/Br, KLM, r = 0.7(2), ■- pp, UAl, r = 0.5(3) f52|. 



Fig. 4.17 X q as a function of q for a) KLM, CMS and NA22 results P5]P^ and b) 
central C-Cu collisions at 4.5 A GeV/c [|159|| . 



Fig. 4.18 Bose-Einstein correlation of order 2 to 5, as indicated. The dashed lines 
represent fits by Gaussian terms, the full lines by exponential terms. All data are corrected 
for Coulomb interaction [ 206|| . 



Fig. 4.19 The dependence of R 2 = N 2 /N2 G on Q 2 for all pairs of charged particles 
(full line), for opposite-charged pairs (dotted) and for same-charged pairs (dashed) 
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Fig. 4.20 a) hiF q as a function of — \ri5y for various pr cuts as indicated [ |142 
b) Anomalous dimensions d q as a function of the order q, for various pt cuts as indi- 
cated (lines are to guide the eye) ||142|| , c) slope 02 as a function of the average transverse 



momentum in an UAl event compared to PYTHIA 5.6 |[210[| . 

Fig. 4.21 Factorial moment F 2 and _F 3 as a function of resolution for three e + e~ data 
sets with pr cuts as indicated | |136[ |. The lines correspond to the models as indicated. 



Correction factors are given above the corresponding sub-figures. 

Fig. 4.22 F2 and F% of the first, second and third jet ordered by their energy. DELPHI 
corrected data (open symbols) are compared with JETSET 7.2 PS Monte-Carlo predic- 
tions with re-tuned settings [|136|1 . 



Fig. 4.23 a) Multiplicity dependence of the slope 03, compared to that expected from 
a number of models, the crosses correspond to a combination of independent events [ |143| 



b) slope 02 as a function of particle density for NA22 (hp at 250 GeV) and heavy-ion 



collisions as indicated | 155 |. The full line corresponds to an extrapolation from the NA22 
point to higher densities using 2 oc 1/p. The dotted lines show fits for the 16 and 28 Si 
samples and for the 32 S samples. 

Fig. 4.24 a) UAl factorial moments (crosses) decomposed into cumulant contributions: 
lower squares indicate 2-particle, upper squares 2+3-particle contributions [ |217| ], b) Upper 
row: normalized factorial moments F q (□) from NA22, together with contributions from 
2-particle correlations (A), the sum of 2- and 3- particle correlations (o), and the sum 
of 2-, 3- and 4-particle correlations (<^); lower part: the normalized cumulant moments 



K q [|142|| . c) K 3 for nuclear collisions [^T 



Fig. 4.25 hiF pq as a function of — In D for four values of Sy, as indicated 



Fig. 4.26 Dependence of lnF pq on the bin size 5y for a correlation distance D = 0.4, 
a) for NA22 data, b) for a sample of 60 000 FRITIOF Monte-Carlo events. The dashed 



lines correspond to horizontal- line fits through the points ||223 



Fig. 4.27 a) b) The increase of the slopes P? with increasing order pq compared to 
the expectation from FRITIOF, for two values of Sy, respectively; c) d) the increase of 
pg /02 with increasing order pq, compared to that expected from the a model (dashed 
line), for two values of 5y, respectively ||223|| . 



Fig. 4.28 a) Fn versus D (5y = 0.2) and F2 versus Sy for a Gaussian-shaped two- 
particle correlation function ||224|| compared to NA22 data ||223|| ; b) F12 versus D (5y = 0.2) 



as m a 



Fig. 4.29 a) (lnG q ) as a function of resolution (M = 2 M ) for UAl compared to the 
expectations from the Monte-Carlo models GENCL and PYTHIA, b) (a q ) and (r q ) as a 



function of q, c) the multifractal spectral function {f(a q )) as a function of (a q ) p33| . 
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Fig. 4.30 a) The universality function T q as a function of £ = \i — v for q = ±5 for 
UAl, b) the spectral function (f(ce q )) as a function of (a q ) [|233j| . 



Fig. 4.31 a)-c) numerical results based on the Bernoulli-trial model (|4.37|) averaged 
with a Negative Binomial distribution P(n) with (n) = 6.26 and dispersion 2.94 [ f230| . 
The line in (c) is a fit to r q = — C (q — 1). 

Fig. 4.32 Comparison of the exponents 4> q of F q and 1 — g — r^ yn (respectively 1 — q — r q ) 
of G q in a) ECCO Monte-Carlo simulations fl24T| and b) tc~ AgBr data at 350 GeV/c [p2]| . 

Fig. 4.33 a) Plot of log P vs — log 5y for n = 6 and 14. The data correspond, down- 
wards, to thresholds nth = M 7 =l,2,3,4,5; b) Codimension functions c(j) vs degree of 



singularity 7 for n = 6 and 14 [51 



Fig. 4.34 a) The integration domain fl B = £ m f2 m of Pziyi, 2/2) for the bin-averaged fac- 
torial moments, b) the corresponding integration domain Q$ for the density integral [246], 
c) illustration of a g-tuple in snake topology, d) GHP topology, e) star topology. 

Fig. 4.35 a) The fourth factorial moment of the NA22 spike event [ 119| , b) the density 
strip integrals for q = 2 — 4 | 246| . 

Fig. 4.36 Density strip integrals C q for q = 2 — 5 for a) all-charged and negatives in 
7r + p and K + p collisions at a/s=22 GeV (the integral F 2 S is also given for (H — ) combina- 
tions) ||142|| , b) all-charge and same-charge combinations in pp collisions at 630 GeV [|143|| . 

Fig. 4.37 In Ff in the NA22 data (full circles) compared to FRITIOF2.0, FRITIOF2.0 
with Dalitz decay and FRITIOF2.0 with Dalitz decay and 7-conversion (open symbols), 
for (cc), ( — ) and (H — ) combinations, as indicated [11421 . 



Fig. 4.38 The slope 02 as a function of average transverse momentum pp and multi- 
plicity n for UAl data and PYTHIA [gTOfl - 

Fig. 4.39 Comparison of density integrals for q = 2 in their differential form (in inter- 
vals Q 2 ,Q 2 + dQ 2 ) as a function of 2 log(l/Q 2 ) for e + e~ (DELPHI) and hadron-hadron 

(uai) m. 

Fig. 4.40 Density integrals F 2 S (in their differential form) as a function of Q 2 for 
like-charged pairs in UAl ||143|| and NA22 1 142|| , compared to power-law, exponential, 
double-exponential and Gaussian fits, as indicated. 

Fig. 4.41 K2(Mi nv ) for cc, H — , and ++ pairs of tracks with cm. rapidity —2 < 

y < 2, in K + /7r + p collisions at 250 GeV/c [261 . The solid line is a power-law fit (see 
text). 

Fig. 4.42 K 2 (M inv ) for a) H — , b) pairs of tracks with cm. rapidity —2 < y < 2, 

in K+/7T+P collisions at 250 GeV/c pgj| , compared to FRITIOF2. 
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Fig. 4.43 R = K2(Mi nv ) + 1 for e + e~ annihilation compared to the JETSET7.3 pre- 
diction with parameters tuned to the DELPHI data for a) unlike-sign and b) like-sign 
combinations [ |262| |. 

Fig. 4.44 In K*(Q 2 ) as a function of — \nQ 2 for all charged particles as well as for 
like-charged particles ||267 |. 



